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ON WAVES OF FINITE AMPLITUDE IN DUCTS 
PART I. WAVE FRONTS 
By R. E. MEYER 
(Department of Mathematics, The University, Manchester) 
[Received 2 May 1951] 


SUMMARY 

The shock-free flow of an inviscid, perfect gas in a duct of varying cross-section 
is investigated. It is assumed that the area of cross-section changes slowly enough 
for the variation of the velocity over any cross-section to be negligible. 

When a steady, shock-free flow is set up in such a duct, and the pressure at the 
entry, or the exit, of the duct is then lowered, or raised (in such a way that a shock 
wave is not formed immediately), then ‘wave fronts’ are formed which separate the 
region of steady flow from the region of unsteady flow. It is assumed that the rise, 


or fall, of the pressure is initiated with a discontinuity in the local rate of change 
of pressure. 

The equations of motion are integrated exactly on such wave fronts (sections 3, 4), 
and the tendencies for different wave fronts to form shocks, or decay, or to approach 
an asymptotic form are discussed. The results are summarized in section 5; they 
explain the failure to establish, in practice, steady shock-free flows decelerating 
through the sonic speed. 


1, Introduction 
Iy the theory of shock-free flow of an inviscid gas a new solution of the 
equations is obtained from any known solution if the sense of the velocity 
is reversed at every point. In particular, if a converging-diverging duct 
issuch that a steady shock-free flow is possible in which the gas is acceler- 
ated from a subsonic to a supersonic velocity, then theoretically a shock- 
free flow in the opposite direction is also possible in which the gas is 
decelerated from supersonic to subsonic speed. The former (accelerated) 
flow will be called ‘nozzle flow’, the latter (decelerated) flow ‘diffuser flow’. 
In practice only nozzle flow has been observed, and from these observa- 
tions it has been inferred that the corresponding shock-free diffuser flow 
isunstable. For the understanding of the mechanism of compressible fluid 
motion it is a matter of interest to determine whether this instability can 
be explained from the theory of inviscid fluid motion or whether boundary- 


layer effects must be taken into account in any explanation. Moreover, 
if the theory of inviscid fluid motion suffices, it is still of some interest to 
determine whether it is necessary to consider the full, two- or three- 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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dimensional flowy in such a duct, or whether it is possible to explain the 
instability by the help of the simple theory of plane waves in a duct of 
slowly changing cross-section. 

The physical argument commonly advanced in favour of the last-named 
alternative is as follows. In a compressible fluid small disturbances ar 
propagated with the local speed of sound relatively to the moving fluid. 
and any disturbance therefore gives rise to an ‘advancing’ wave and 4 
‘receding’ wave, which reach particles upstream and downstream, respec. 
tively, of the particles originally affected by the disturbance. As seen by 
an observer at rest, advancing waves always move downstream, whereas 
receding waves move downstream in supersonic flow and upstream in 
subsonic flow. In both nozzle flow and diffuser-flow, advancing waves 
originating downstream of the throat pass immediately out of the duct 
and those originating upstream travel through the throat and exit at high 
speed. Moreover, in nozzle flow, all receding waves move away from the 
throat, wherever they originate. On the other hand, in diffuser flow, all 
receding waves travel towards the sonic line, at which they accumulate 
asymptotically, since their absolute speed of propagation tends to zero as 
the local Mach number approaches unity. Thus only in diffuser flow, and 
not in nozzle flow, do the circumstances favour the accumulation of 
disturbances. 

Attempts to extend this argument to a proof of the instability of 
diffuser flows were made by Brown (1) and Kantrowitz (2). However, 
their results appear to be contradictory, and the arguments themselves 
are open to serious criticisms (see Appendix). 

In the following, attention is restricted to the fronts of the waves, where 
the disturbed region borders on the still unaffected, steady flow. It is 
assumed that the acceleration of the fluid particles changes discontinuously 
when the wave front passes over them; the case where the acceleration is 
continuous is treated in Part II, and leads to similar results. 

In addition to nozzle flows and diffuser flows, the analysis applies to all 
steady shock-free flows where the speed is subsonic throughout the duct, 
or supersonic throughout the duct, including the extreme cases where 
sonic speed is just reached in the throat. It applies also to ducts without 
throat and to ducts with several throats. 


2. The characteristic equations 
Consider the unsteady, isentropic flow of an inviscid, perfect gas in a 
duct whose cross-section changes slowly enough for the velocity to be taken 


+ A stability investigation for two-dimensional, transonic flow has been carried out by 


Kuo (5). 
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as constant over any plane normal to the axis of the duct. The motion may 
then be described in terms of two independent variables, x and t, denoting 
respectively the distance along the axis and the time. The momentum 
equation may be written ((3), p. 28) 
ou ou , lop 
+u—+4-— — 0, (1) 
ot Cx pox 
where wu is the velocity in the direction of x increasing, p the density, and 
y the pressure. The equation of continuity is 
c C 
Q- += (puQ) = 0, (2) 
ot Ox 
where Q is the area of the cross-section. Since the flow is isentropic, 
dp = a* dp, (3) 
where a is the local speed of sound. When the dependent variables are 
expressed in terms of Riemann’s variables, 





y 1y+41 (a/p) dp . mae Y.. 
~ oe y—1 2 
p 9 (4) 
ind 8 ly 1 ( (a p) dp a _u 
- a y—l 2 


0 
where y is the ratio of the specific heats, the equations of motion, (1) and 


2), are transformed into 


ér/et+-(u+-a) ér/éx = —lauH |\ (5) 
and és/et+-(u—a) és/ex = —}auH , 
where H = (1/Q) dQ/dz, 
and by (4), 
“=r—s ste ae DO “Ys 
- air (6) 
a=! (r-+-s) u—a r 8 


In the plane of the variables x and ¢ (Fig. 1) the lines given by 
dx/dt = ut+a 
and dx/dt = u—a 
are called advancing and receding Mach lines, respectively. In this plane, 
ordinary differentiation with respect to time along a line of slope 
dt/dx O(a, t) 


- « is measured downstream from the throat (Fig. 1), and u is assumed to be positive 
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Fic. 1. 1. Advancing wave front. 2. Receding wave front 
in diffuser flow. 
therefore represent ordinary differentiation with respect to time on the 
advancing and receding Mach lines, respectively, and by (5), the equations 
of motion may be written 


D,r/Dt = —tauH, (9) 
D_s/Dt = —tauH. (10) 


The derivatives D,s/Dt and D_r/Dt, on the other hand, are not deter- 
mined in terms of a, wu, x, and ¢ by the equations of motion without 
reference to particular boundary conditions, and these derivatives can be 
discontinuous. A differential equation for the variation of D,s/Dt ona 
receding Mach line is, however, obtained as follows. By differentiating 
(10) with respect to time along an advancing Mach line, taking account 
of (6) and (9), one finds that 

D, D_s 

Dt Dt 
Moreover, by (7), (8), (6), (9), and (10), 

F Bp,  @, Bil A (55 (u+a) Fe (ua) 


= }(y—1)H(sauHr-+sD, s/Dt)—fau(u+a)dH/dx. (Il) 


Di Di Dt Di} @\Dt Dt 


_ ytl(D,s auy\(D_r , Dis , 3-y H 12) 
a (Het 4) (B+ 3 a F * 
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By adding these two equations, a first-order, ordinary, non-linear differen- 
tial equation is obtained for D, s/Dt as a function of the time on a receding 
Mach line. The coefficients depend on the distribution of r, s, x, and ¢ on 
this Mach line. A similar equation governs the variation of D_r/Dt on an 
advancing Mach line. 

In the following, the unsteady flow through converging-diverging ducts 
is studied for which there exists a steady-flow solution without shock. This 
steady solution is given by Reynolds’s ‘hydraulic’ theory and the velocity, 
local speed of sound, and Riemann variables in this solution, denoted by 
U, A, R, and S respectively, may be considered as known functions of zx. 
They satisfy the equations (4) and (6) and ((3) p. 379) 


AUH = (U—A)(U+A)U'/A 


, (13) 
and A’ 4(y—1)UU'/A 


where a dash denotes ordinary differentiation with respect to x of a 
function independent of time. It will be assumed that H = (1/Q) dQ/dx 
is continuous and vanishes only at a throat, and that H’ is bounded. 


3. Receding wave fronts 

Suppose a steady shock-free diffuser-flow has been set up in a converging- 
diverging duct, and assume that a disturbance is generated at the exit only, 
in such a way that the start of the disturbance is marked by a discon- 
tinuity of the acceleration (but not of the velocity itself). The disturbance 
spreads in both directions. The downstream wave front travels away from 
the duct, but the upstream wave front—whose path in the (z, ¢)-plane is 
represented by a receding Mach line (Fig. 1)—travels into the duct and 
separates the disturbed region from the steady flow upstream. As long as 
this wave front does not meet a shock-wave, the velocity and state is 
continuous across it, and the coefficients in equations (11) and (12) may 
therefore be replaced by their local values in the steady solution. When 
this is done (use being made of (6) and (13)), and the two equations are 
added, the relation 


D_o/Dt a o(o+ F(x))+G(zx) (14) 


is obtained on the wave front, where o stands for D, s/Dt, and 


F(z) _U" _144.1)(U—A)*—(9—3y + 4U/A)(U2—A2)]. (15) 


2(y+1)A 


A particular solution of (14) is known, for on the upstream side of the 
wave front o must equal its value in the steady solution, 


= = D,S/Dt = (U+-A)S’ = —}A(1+U/A)?U" (16) 
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(by (7), (6), and (13)). By (14) the value of o on the downstream side must 
therefore satisfy the equation 


(4A/(y+1))D_(o—X)/Dt = (o—)[o+2+ F]. (17) 
With (c—X) = 1/z, this becomes 
D_2z/Dt = —((y+1)/(4A))[1+(22+ F)z], (18) 


and upon integrating, 
t 





Kz = Kx — 2 (K/A) dt 
0 (19) 
where a/, = exo { = i 


0 
both integrals are taken along the path of the receding wave front in the 
(x, t)-plane, and values taken on this wave front at time ¢ = 0 are dis- 
tinguished by a suffix 0. Let WM = U/A; then, by (15) and (16), 


22+ F = —(3+1/M)[2+(y—)MJAU'/(y+1), (20) 
which is bounded, and therefore z must vanish after a finite time when 
%, > 0, but z does not vanish when z, < 0. 

Since 
Du/Dt—DU/Dt = (é/ét+-ué/éx)(u—U) = —}(o—) = —1/(22) (21) 


on the receding wave front, by (7) to (10), (6), and (13), a zero of z implies 
an infinite acceleration and hence, a breakdown of the isentropic flow and 
the formation of a shock. A positive value of z corresponds to a sudden 
drop in the acceleration of the fluid particles as the wave front passes over 
them, by (21), and to a sudden rise of the local rate of change of pressure, 


— ép/ét = 4p(A—U)(o—3) (22) 


on the receding wave front, by (1) to (4) and (7) to (10). Therefore, if 
a gas flows steadily through a converging-diverging duct, with sonic 
speed at the throat and subsonic speed in the diverging part, then any 
disturbance at the exit that starts with a sudden rise of the local rate of change 
of pressure must lead to shock-formation after a finite time. 

To calculate the time ¢, after which the limit point appears on the wave 
front, note that K = Mii—_M)2, 


since da/dt = U—A on the receding wave front and by (19), (20), and 
(13); and hence, 


Miz/(1—M)? = [Mz/(1— My] + css (1—M)-8A-2 dx, (28) 
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where 2, denotes the position of the wave front at time ¢ = 0. If x, denotes 
the value for which the right-hand side of (23) vanishes, 


Xo 

The asymptotic behaviour of z is best inferred from (18). We may 
assume that U’ 0 in a diffuser and hence (22+ F) > 0, by (20). There- 
fore, z increases with time when z < —1/(22+-F), and decreases with time 
when 2 1/(2X+ F), and asymptotically, that is, as the receding wave 
front approaches the throat, (1/z) = o— =X approaches the asymptotic 
value of —(2X+-F) on the wave front. But F > 0, by (15), and hence, 
whenevero ~ Linitially,o + —> asymptotically, so that the local pressure 
sradient, @p/éx, just downstream of the wave front, tends to the opposite 
of the local pressure gradient in the original steady flow, just upstream 
of the wave front, by (22).+ This asymptotic solution can have physical 
significance only when the flow remains isentropic near the wave front. 
When z, < 0, i.e. when the disturbance starts with a sudden drop of 
the local rate of change of pressure, the acceleration remains finite on 
the wave front; but knowledge of the further development of the dis- 
turbance is required in order to exclude the possibility of a shock forming 
in the interior of the disturbed region and overtaking the wave front. 
The asymptotic solution shows, however, that the steady, shock-free 
diffuser-flow cannot be re-established even when the flow remains isentropic 
for all time in the neighbourhood of the wave front. 

3.1. The same results apply to the upstream wave front of any dis- 
turbance generated at the entry of a converging-diverging duct (with the 
end of the disturbance marked by a discontinuity of the acceleration) 
when the initial flow in the converging part is steady, supersonic, and 
shock-free, with sonic speed at the throat. The Mach line in the (z, t)-plane 
that separates the wave from the undisturbed region upstream is in this 
case the receding Mach line passing through the entry of the diffuser at 
the time at which the disturbance ends (Fig. 2). Again, shock formation 
is certain when the fluid particles suffer a sudden drop of acceleration as 
the wave front passes over them; by (22), this means that the local rate 
of change of pressure is negative during the ultimate stage of the dis- 
turbance. 

3.2. The development of receding wave fronts propagated into an 
entirely subsonic, or entirely supersonic steady, shock-free flow in a con- 
verging-diverging duct, where sonic speed is not reached in the throat, 
may be discussed similarly. Take first the case where the flow is subsonic, 


+ This result was first suggested by Kantrowitz (2). 
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and the disturbance is generated at the exit. Then the receding wave front 
travels through the throat into the converging part and leaves the duet 
through the entry, and this restricts the time during which the solution 
(19) has physical significance even if the flow remains isentropic in the 
neighhourhood of the wave front. 

When z) < 0, i.e. when the disturbance starts with a drop of the local 
rate of change of pressure, the acceleration remains bounded on the wave 
front, by (19) and (21). The sense of variation of o on the downstream 


t 








i gall 
> 


Fic. 2. 1. Advancing wave front. 2. Receding wave front 
in diffuser flow. 
side of the wave front is towards the local value of —(=+ F), in the diver- 
gent part of the duct, by (18) and (20); in the converging part, where the 
steady flow is accelerated, and so (25+ F) < 0, it is towards &, that is, 
the wave front loses strength. 

When 2) > 0, there is a tendency to shock-formation in the diverging 
part; if this has not occurred when the wave front enters the converging 
part, the same tendency will persist, at least for some time (by (i8) and 
(20), since U’ = 0 at the throat). But later a different influence may 
prevail; to study this it is convenient to suppose that the duct extends 
indefinitely in the upstream direction. 

Far upstream, U is small and UQ is approximately constant, and hence 
the integral in (23) remains bounded if |a-2Q| > 0 asa—> —oo. Therefore, 
even if the converging part extends indefinitely upstream, the acceleration 
will remain bounded and tend to zero ultimately (by (21) and (23)) and 
shock-formation may never occur, provided (i) the duct converges more 
rapidly than a cone, at a great distance upstream, and (ii) the initial value 


of z is sufficiently big, i.e. the initial rise in the local rate of change 0 


pressure at the exit is sufficiently small. 
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Next, take the case where the original, steady flow is entirely supersonic, 
ind the disturbance is generated at the entry and is of finite duration. 
Then the receding wave front separating the disturbed region from the 
steady flow upstream travels through the throat and finally passes out of 
the duct through the exit. The behaviour of z is similar to that in the 
previous case, but the tendency to shock-formation when z, > 0, is more 
lifficult to overcome; when A is small, QA*'Y-) is approximately constant 
by (13)), and |2-*"-)Q| must tend to infinity as x— oo, in order that the 
integral in (23) should remain bounded. For air with y = 7/5 that means 
that the area of cross-section must tend to infinity more strongly than the 
tenth power of the distance from the throat. 


4, Advancing wave fronts 

Consider now a converging-diverging duct in which a steady shock-free 
flow has been set up, and assume that a disturbance is generated at the 
entry only, with its start marked by a discontinuity of the acceleration. 
The path of the wave front that separates the disturbed region from the 
undisturbed, steady flow downstream is represented in the (x, t)-plane by 
the advancing Mach line passing through the entry at the beginning of the 
disturbance (Fig. 2). Whatever type of steady flow has been set up to 
begin with—diffuser-flow, or nozzle-flow, or entirely subsonic, or entirely 
supersonic flow—this advancing wave front travels through the throat 
ind ultimately passes out of the duct. 

On the wave front, r, s, and x are again known as functions of time, and 
in the same way as equation (14) for o, we find for the variation of 
\= D_r/Dt on the wave front the equation 


Dr) Dt = ((y+1)/(4A4))A(A+ P(x))4+-V (a). 


\ particular solution, A D_R/Dt, is again known, so y = 1/(A—A) 
satisfies the equation 


D.y/Dt ((y+-1)/(4A))[1+(2A+ P)y], (24) 
und a short calculation shows that 
2A4+P = (1/M—3)[2+(y—1) MAU (y+). (25) 
The solution of (24) is 
t 
ott f 
Li a eee L/A) dt 
Yy oo 4 | ( ye 
, , (26) 
Oe | e¢ it > 
where L/L, exp)? , | 2A f at| 
pe Joa 4 
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and both integrals are taken along the path of the advancing wave front 
in the (x, t)-plane. Therefore, (U+ A) dt = dx, and by (25) and (13), 


L= M*\(14+M)-, 
and, by (26), 


Miy/(1+M)? = [My/(1+M)°] 


Ey 
y+1 4 j | ; Sao to 
—= | M41+4+M)-A-2 dx. (21) 
Lo 
This equation shows that there is a definite tendency to shock-formation 
when y, > 0 in both diffuser-flow and nozzle-flow and also all other steady 


flows. Since : 
Du/Dt—DU/Dt = 4(A—A) = 1/(2y) (28) 


and ep/et = $p(U+ A)(A—A) (29) 


on an advancing wave front, by (1) to (4), (6), (7) to (10), and (13), y >0 
means a sudden rise in the acceleration of the fluid particle, as the wave 
front passes over it, and y) > 0 for a disturbance starting with a rise in 
pressure. Note, however, that the wave front may have passed the exit 
of the duct before the acceleration has become large on it; and even if a 
shock forms at the front of the wave while the front is still inside the duct, 
the shock travels with an absolute speed greater than the local value of 
(U-+-A) just downstream of the shock and passes the exit after a very 
short time. 

Moreover, even if the duct extends indefinitely in the downstream 
direction, the tendency to shock formation may ultimately be overcome 
if (i) the integral in (27) converges as x —> 00, and (ii) the initial rise in the 
local rate of change of pressure is sufficiently small. For a duct with 
steady subsonic flow in the divergent part, this integral converges if the 
duct diverges ultimately more rapidly than a conet (ef. section 3.2 above). 
For a duct with steady, supersonic flow in the divergent part, the integral 
converges only if |a-*"7-)Q| > o as x> oo. If the tendency to shock- 
formation is indeed overcome, (A— A) tends to zero asymptotically in all 
ducts. 

When y, < 0, on the other hand, y remains negative, and the accelera- 
tion remains finite on the wave front. Equations (24) and (25) show that, 
when M > 4, |A—A| decreases with time in a nozzle; but in a diffuser the 
sense of variation of A is towards the local value of —A—P, so long as 
M > 3. For any other duct, the tendency in the converging part is the 
same as in a nozzle or a diffuser according to the sign of U’; and similarly 
for the diverging part. 


+ Boundary layer separation may prevent this. 
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5, Conclusions 

In order to summarize the results obtained for the development of wave 
fronts travelling into a steady, shock-free flow in a duct of slowly varying 
cross-section, it is convenient to classify the wave fronts in two different 
ways, namely (i) as to whether they are connected with an expansion or 
with a compression, and (ii) as to whether they travel through the duct 
and pass out of it, or whether they travel towards a sonic throat so as to 
approach it asymptotically. 

All compression wave fronts, i.e. all wave fronts raising the local rate 
of change of pressure when passing through a section of the duct, show at 
least an initial tendency to shock formation. Ali expansion wave fronts 
show a tendency opposed to shock formation. In fact, when they travel 
into a steady accelerating flow they show a permanent tendency to decay. 
But when they travel into steady decelerating flow they show at least an 
initial tendency to establish the local pressure gradient appropriate to the 
accelerating steady flow solution, i.e. the opposite of the local pressure 
gradient of the original steady flow. 

Wave fronts are called ‘advancing’ or ‘receding’ according to whether 
they travel downstream or upstream, with respect to the moving fluid. 
Advancing wave fronts generated at the entry of a duct travel through 
the duct and pass out of the exit after a short time; advancing wave fronts 
generated at the exit do not enter the duct at all. Receding wave fronts 
generated at the entry travel into the duct if the steady flow is supersonic 
near the entry, and penetrate as far as it remains supersonic. Receding 
wave fronts generated at the exit travel into the duct if the flow is sub- 

mic, and penetrate as far as it remains so. If there is a throat where the 
steady flow is sonic, then all receding wave fronts entering the duct 
approach the (nearest) sonic throat asymptotically. 

The tendencies of development found for compression and expansion 
waves are not only initial, but permanent, for all wave fronts approaching 
asymptotically a sonic throat. For instance, if a steady shock-free diffuser 
flow could be set up, then any disturbance generated subsequently at the 
exit so as to start with a rise of the local rate of change of pressure must 
inevitably lead to the formation of a shock after a finite time (section 3). 

For wave fronts travelling right through a duct, conditions are compli- 
cated by the fact that a new tendency makes itself felt as soon as the wave 
front has passed through a throat and this may counteract the initial 
tendency. To clarify this question it is convenient to suppose that the 
duct extends indefinitely in the direction in which the wave front is 
travelling. It is then found that the initial tendency must prevail if, 
asymptotically, the duct does not diverge more strongly than a cone, in 
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the case of a steady subsonic flow, or if the cross-section does not tend to 
infinity more strongly than the tenth power of the distance from the throat (i) Br 
(for air), in the case of a steady supersonic flow (section 3.2). Otherwise, introduc 
a tendency for the wave front to decay must ultimately prevail, provided ange 
the wave front has not been overtaken by a shock. For instance, suppose 
a steady nozzle flow leading to a very low pressure has been set up ina 


a ‘ , , . (equatio 

duct the supersonic portion of which diverges more rapidly than the tenth aa 
“ ‘ 3 . oi: a op ec 

power (for air) of the distance from the throat. Then there exists a limiting Brown | 


initial strength for any compression wave front (not, in general, the same | @/2één 


for different ducts); if it is stronger, a shock must form before the tendency J reedin’ 
. . ee ‘limi 

to decay can prevail over the earlier tendency to steepen; if it is weaker, - ; 
no shock is formed in the neighbourhood of the wave front and, asymptoti- 2a 
coe 


cally, the strength tends to zero (section 3.2). 
These results are seen to explain the difference in instability between 
steady nozzle flows and steady diffuser flows. From a practical point of 


(use is | 
view, the contrast between the two types of flow is due, in the first place, aan 
not to the ‘accumulation’ of disturbances, which is a property of com- és/e§ an 
pression waves, but to the fact that in nozzle flows all wave fronts, and (ii) 
the shocks to which they may lead, are blown out of the duct immediately, steady 

, : : : . ‘ (M—1) 
In diffuser flows, on the other hand, all receding wave fronts, and for With r 
sufficiently small disturbances also their shocks (4), are caught in the duct. isintro 
The distinction is accentuated by the tendency even of expansion wave is then 
fronts not to decay when travelling into a decelerating, steady flow. ane 
From a mathematical point of view, it is of interest to note that the _ 
present stability investigation for a non-linear partial differential equation 
leads to results resembling those of stability investigations for non-linear, by (8), 
ordinary differential equations leading to ‘limit cycles’.t If attention is of (2) 
confined entirely to the equations of isentropic flow, the pressure gradient pulse. 


just downstream of any receding wave front with discontinuous accelera- 

tion that is caught in a duct is seen to tend asymptotically to the local 

pressure gradient of the steady accelerating flow. This steady solution 

shows some analogy to a ‘stable limit cycle’; the steady, decelerating flow, 

on the other hand, shows the same analogy to an ‘unstable limit cycle’. 
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ON WAVES OF FINITE AMPLITUDE IN DUCTS 
APPENDIX 


(i) Brown (1) studies wave fronts of the same type as those treated above. He 
ntroduces characteristic coordinates £, 7, where € is constant on receding Mach lines, 
ind 7 on advancing Mach lines, so that the equations (7) and (8) may be replaced by 


éx/e&—(u+a) et/e& = 0, ér/e&+4auH et/eé = 0, (A 1) 
ex/en—(u—a) ét/en 0, és/en + 4auH Ct/en 0 (A 2) 


ations 5.5 and 5.8 of (1); the definition of s given in (1) differs from the one 
\dopted here by the sign). To calculate the steepening of a receding wave front, 
Brown assumes that r, s, er/c n> C8/CNH, %, t, Ox e€, ct o€, Cx on, ot on Cn o&@y, and 
"/@één are continuous, but that @s/c&, at least, is discontinuous across a certain 
eceding Mach line & E*, 


Elimination of x from the first pair of the equations (A 1) and (A 2) gives 


ct ct ¢ ct « 
Ja — : (u a)— = (u a) 
c&en c& Cn cn c& 
, Ot he Y Hi y+1 @r ot yt+1 et es 
za — au - + -———— — ] -—4+- —_—_ -—- SCS 
a 9 9 ¢ ‘ t 5) 7 7 
cécy v4 cn = ©€”H c& 4 én O€ 


ise is made of (3), (A 1) and (A 2)), which shows that the above assumptions are 
neompatible. In fact, a discontinuity of D,s/Dt implies a discontinuity of both 
€ and ét/e& (and hence, also of ér/é& and cx c€). 

ii) Kantrowitz (2) studies the whole of a pulse superimposed on the flow of the 
steady solution in the subsonic part of a diffuser, neglecting terms of the order of 
M—1) against unity, where M is the local Mach number of the steady solution. 
With regard to such an approximation it should be noted that, if the approximation 
isintroduced into an equation like (7), or (8), or (11) and (12), and such an equation 
s then integrated, it is difficult to assess the asymptotic behaviour of the resulting 
expression without further investigation of the effect of the small terms over a very 

ng time. In fact, on a wave front of the type studied in section 3 above 


0 


{ (U—A) dt = [ dx = O(1), 


(8), if the wave front does not meet a shock wave. As it stands, the argument 
f (2) is not certain to lead to a correct description of the asymptotic state of the 


pulse. 
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PART II. WAVES OF MODERATE AMPLITUDE 


By R. E. MEYER 
(Department of Mathematics, The University, Manchester, and Institute for 
Fluid Dynamics and Applied Mathematics, University of Maryland) 


[Received 9 October 1951] 


SUMMARY 

The results of Part I are extended to a first-order theory of advancing and receding 
waves travelling into steady, shock-free flow of inviscid, perfect gas in a duct of 
slowly varying cross-section. This includes, as special cases, spherical waves and 
waves with cylindrical symmetry. 

The theory applies (i) to the neighbourhood of any wave front travelling into 
steady, shock-free flow, and (ii) to the waves caused by any disturbance of finite 
duration and extent such that the velocity and pressure differ little from their 
respective local values in the steady flow (section IV. 2). No assumption is implied 
that the local rates of change of velocity and pressure are small, and shock waves 
may form. Only a partial solution, however, is given for the interaction of an 
advancing wave with a receding wave of comparable strength (section IV. 1). 

In a first approximation, the fluid acceleration, the local rate of change of pressure, 
and the limit lines are found (sections III. 2, IIT. 4). A second approximation gives 
the velocity and the pressure (sections IV, IV. 1) and the shock path, to the first 
order in the shock strength (sections V, V. 1). 

I. Introduction 

XIEMANN (1) investigated the one-dimensional, unsteady flow of an in- 
viscid gas and found that any disturbance of finite extent and duration 
set up in an originally steady flow generates two waves, one advancing and 
the other receding with respect to the moving fluid, which are propagated 
into the gas not initially affected by the disturbance. After a certain period 
of interaction (depending on the way in which the disturbance is set up), 
these two waves will become separated from each other, and each will 
then be a ‘simple wave’. Between the two waves the flow will again be 
steady. 

Conditions are similar in the unsteady flow of inviscid gas in a duct, the 
area of cross-section of which changes slowly enough for the velocity and 
pressure to depend primarily on two independent variables only, viz. the 
position, a, and the time, ¢. If a disturbance? of finite extent and duration 
is set up in an originally steady, shock-free flow, it generates two primary 
waves, one advancing and the other receding, which become separated 

t It is assumed that a shock wave is not formed immediately. That case, however, 
can be treated by the method described below, provided the shock is weak (cf. 2). 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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ter a certain period of interaction. These waves are not simple waves. 
Owing to the variation of the area of cross-section, the waves are con- 
tinually distorted, as they travel into gas originally in steady flow, and the 
process Of distortion is accompanied by one of refraction. The receding 
wave thus generates an advancing wave by refraction, and the advancing 
wave generates a receding wave. The refracted wave interacts with the 
primary wave, gets itself distorted as it travels through the duct, and by 
the accompanying process of refraction contributes to the primary wave. 
[he primary wave passes any station in the duct within a finite time, but 
it leaves behind it the refracted wave as a ‘tail’, and the flow between the 
two primary waves is not steady, but rather a superposition on a steady 
flow of the interaction of the two tails accompanied by the effects of distor- 
tion and refraction of each of them. 

The original steady flow does not give rise to any wave, and since the 
process of refraction is a continuous one, the strength of the refracted wave 
is built up only gradually with distance irom the front of the primary wave 
travelling into steady flow. This makes it possible to treat certain types of 
transient waves by successive approximations. A detailed discussion of 
the type of wave to which the theory applies is given in section IV. 2. 

It is found convenient to take as the fundamental dependent variables 
not the velocity and pressure but certain characteristic derivatives of them 
which are not determined directly by the equations of motion (or rather 
the differences between these derivatives and their local values in the 
steady flow). Each primary wave is characterized by the distribution of 
one such derivative; the distribution of the other derivative characterizes 
the refracted wave. In what follows, the successive approximations are 
carried out in detail for the case of a receding wave generated by a dis- 
turbance at the exit of a duct. For other cases the results are quoted. Fora 


receding wave, the ‘primary’ derivative is 
Ao = [é/et+(u+a) é/éx|{(a—A)/(y—1)—}(u—U)}, 


where wu is the velocity, a the local speed of sound, and U and A are their 
respective local values in the steady flow. (See list of symbols at the end 
of section IT.) 

For a receding wave, the first approximation consists in the calculation 
of the distribution of Ac (equation 35), which results from the combined 
effect of the growth typical for a simple wave and the distortion due to 
the change in the area of cross-section. This approximation yields also the 
ocal rate of convergence of the receding Mach lines (equation 39), and 
the receding limit lines (equation 40), but the deviations of velocity and 


pressure from their local values in the steady flow are neglected. 
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Once the distribution of Ac is known to the first order, it is found that. 
as a second approximation, these deviations of velocity and pressure can be 
obtained to the first order (section IV). The value of Ao at any one point 
in the wave is found to depend only on the value of Aco in the disturbance 
at the exit at a certain instant of time, but the deviations of velocity and 
pressure depend also on all the earlier values of Ao at the exit. The second 
approximation yields also the shock path, to the first order in the shock 
strength (sections V, V. 1). Explicit formulae for the shock path are only 
given for shocks starting at, or near, the wave front, in which case certain 
effects of distortion due to the change in the area of cross-section can be 
neglected. 

The refracted wave is found to be only a second-order effect, i.e. small 
even compared with the velocity deviation (sections III. 5, IV. 2) and itis 
neglected hereafter, together with the effects of its distortion, its inter- 
action with the primary wave, and the higher-order effect of secondary 
refraction due to the distortion of the refracted wave. 

The second approximation can also be extended to give the velocity 
correction for a region of interaction of two primary waves of comparable 
strength (section IV. 1). 

The feature distinguishing the present theory from an acoustic theory of 
‘waves of small amplitude’ is that values of any magnitude of the fluid 
acceleration and of the local rate of change of pressure are admitted. The 
feature distinguishing it from a full theory of ‘waves of finite amplitude’, 
like that of Riemann (1), is that the velocity and pressure themselves are 
only permitted to differ little from their local values in the steady flow. 

The present theory is similar in many respects to the first-order theory 
of the steady, supersonic, axially symmetrical flow past slender bodies 
of revolution developed by Whitham (2). Both theories are based on 
the same physical ideas, but the mathematical approach is rather different. 
Whitham’s theory is developed from a linearized theory of a nearly uniform 
flow, whereas in the problem studied here not even the basic flow is uniform 
and the theory is developed directly from the non-linear equations of 
motion. On the other hand, a fundamental feature of Whitham’s theory 
is the correct description of the asymptotic behaviour of the primary 
wave, whereas the present theory applies only to waves passing out of the 
duct after a finite time; this is sufficient, since only steady, shock-free flows 
are stable for which all waves have this property. 

Note. The notation follows that used in Part I, and the equations 
are numbered so that numbers up to, and including, (29) refer to equations 


of Part I. For convenience, a list of symbols is given at the end of 


section IT. 
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I]. Remark on the boundary conditions 

Consider the strictly one-dimensional flow of inviscid gas. A well- 
known problem is that of the flow in a long pipe containing a piston and, 
n both sides of it, gas initially at rest and in thermodynamic equilibrium; 
the piston is set in motion with finite acceleration from the time t = 0 
nwards. Then two waves are propagated into the gas, one from each side 
f the piston, and each is a simple wave, since it travels into gas at rest. 
Ineach of the simple waves one of the two Riemann variables (cf. equation 
4)is constant, but it is not the same Riemann variable that is constant in 
both waves, one of which is ‘advancing’, while the other is ‘receding’. 
In particular, a Riemann variable is constant on each side of the piston, 
but it is not the same Riemann variable that is constant on both sides, 
though both surfaces of the piston have the same velocity. 

This simple set-up depends essentially on the assumption that the gas is 
initially at rest. If the gas is initially flowing through the pipe with uniform, 
subsonic speed, we may, for instance, imagine a disturbance set up in 
the following way. A segment of the pipe wall is made of porous material 
and suction is applied from ¢ = 0 onwards. (Suppose that a useful result 
can be obtained even when the ensuing deviations from one-dimensional 
flow are neglected—it is not the purpose of this investigation to study 
whether, and how, a disturbance could be set up that leaves the flow one- 
dimensional.) From each border cross-section of the porous segment of 
the pipe a simple wave will be propagated into the gas, but the flow inside 
the porous segment will be of a more complex type. (That the waves are 
simple waves is easily seen as follows. In the (x, t)-plane, every point outside 
the strip representing the porous segment is connected by a Mach line not 
passing through the strip to a region representing the initial, uniform flow. 
On such a Mach line, one Riemann variable is constant. In a region of 
uniform flow, both Riemann variables are constant.) 

It is not the purpose of this investigation to study the flow inside the pipe 
segment where the disturbance is set up, and it is therefore convenient to 
suppose that the porous segment is replaced by a narrow slit in the wall of 
the pipe. For mathematical convenience, let the width of the slit tend 
to zero. The slit is then represented in the (x, t)-plane by a line, or at least 
segment of a line, x = constant across which both Riemann variables are 
discontinuous. On each side of the line, one Riemann variable is constant, 
but it is not the same Riemann variable on both sides. On each side of the 
line, the distribution of the other Riemann variable can be prescribed 
arbitrarily, and it must be prescribed in order to determine the solution. 
Alternatively, the distribution of its time derivative, or of the primary 


derivative, may be prescribed. 
5092.19 
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Conditions are very similar in a duct of slowly varying cross-section jp 
which a disturbance is set up in an initially steady, isentropic, subsonic 
flow. A wave is propagated into the steady flow on each side of the position 
where the disturbance is set up. It can be shown that, as one proceeds 
along a Mach line from the region of steady flow in the (z, t)-plane into the 
region of unsteady flow, for a short distance, the rate of change of one of 

the Riemann variables (primary deriva- 
t tive) is much larger than the rates of 
V4 change of the velocity and the pressure, 
respectively, while the rate of change of 
the other Riemann variable is much 
smaller even than those of the velocity 
and pressure.t The latter Riemann vari- 
able is the one that would be constant 


7 ee on the Mach line in strictly one-dimen- 


sional flow. 








It is convenient to study separately the 
el two waves generated by the disturbance. 
In what follows, the disturbances will 

i therefore be represented by a line 

Fic. 1. x = constant in the (2, t)-plane, or a 

segment of such a line, on which the 

distribution of the primary derivative is prescribed. It is assumed to 

be bounded, but not necessarily continuous. This is found to imply that 
the velocity and the pressure are continuous on the line. 

When the initial steady flow is supersonic, the situation is complicated 
by the fact that both waves travel downstream with respect to an observer 
at rest. The region in the (x, t)-plane just downstream of the line represent- 
ing the disturbance is therefore a region of interaction of the two waves 
(Fig. 1). However, if the disturbance is of finite duration the waves separate 
eventually. The results mentioned above regarding the orders of magnitude 
of the respective rates of change hold near the wave front advancing into 
steady flow and near the wave front receding into steady flow. Similarly the 
theory described below applies to each of the waves, after their emergence 
from the region of interaction, but for the region of interaction itself only 
a partial solution is obtained here. For the formulae given below to apply 
without modification, it is necessary to prescribe, as a boundary condition, 
the distribution of the primary derivative on a line x = constant outside 
the region of interaction in the (x, t)-plane. 


+ Actually it is not these four variables themselves to which that statement applies, but 
the difference between the variables and their respective local values in the steady flow. 
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List of symbols ; 

7 — distance along duct, measured downstream. 

t — time, measured from beginning of disturbance. 

u = velocity. 

a = speed of sound. 

r= a/(y—1)+4u. 

s= aj|y 1) Su. 

p pressure. 

y = ratio of specific heats. 

U,A, P, R, S = local values of u, a, p, r, s, respectively, in steady flow 
solution. 

A = difference between local values in unsteady and in steady flow, 
e.g. Au u—U, 

M = U/A. 

p. (ut+a)/(U+A). 

p (u—a)/(U—A). 

H = (1/Q) dQ/dx, where Q denotes area of cross-section. 


v Tull. 
: ] : 
J ~ Al A. 


D,/Dt c/et+(u+a)e/ex. 
D_/Dt e/et+(u—a) é/éx. 


o = D.s/Dt. \= D_r/Dt. 

== D,S/Dt. A = D_R/Dt. 

Ao = o—*. AA = A—A. 

z= 1/Ac. y = 1/Ad. 

K = M'\(M—1)-. L = M\(M-+1)-*. 
yv+1 ff Kdzx y+l ff Ldx 

; t J A(A—U) an 4 J A(A+U) 


dr, dr’ defined in section ITT. 3. 
Suffix 0 denotes values taken at exit of duct, for wave receding into 
steady, subsonic flow. 
! denotes values taken at limit point. 
s denotes values taken at shock point. 
| denotes values taken at shock point, on upstream side of shock. 
2 denotes values taken at shock point, on downstream side of 


shock. 
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Ill. The first approximation 

Consider the receding wave generated by a disturbance at the exit of, 
duct in which a steady shock-free flow has been set up for ¢ < 0, with gyb. 
sonic flow near the exit. Assume that the disturbance is such that 

As,/Sy = O(e) (€ <1), (30)+ 
at the exit, for 0 < t < TJ, where 7, may be the time during which the 
disturbance lasts. 

If sy is given for 0 <t < J, the flow is determined, by the general 
Uniqueness Theorem (3), between the wave front receding into the steady 
flow and the receding Mach line passing the exit at time 7). Therefore, fit 
is possible to find a consistent, approximate solution on the basis of the 
additional assumption that 

Au/U = Ofe), Aa/A = O(e), AAJA = O(c) (€ <1), (31 
everywhere in the region where the flow is so determined (and this will be 
done in what follows), it may be concluded that the assumption (31) is 
justified if the assumption (30) is accepted. 

III. 1. Equations (4) and (6) of Part [remain valid in a region of unsteady 
flow when capital letters are substituted throughout, but not equations 
(9) and (10); in fact, 

D, R/Dt = —p,V, D_8/Dt = —p_V, (32) 
by (7), (8), (4), and (13). An equation for the variation of o = D, s/Dt on 
a receding Mach line can be established in a way similar to that given in 
sections 2 and 3 of Part I; after a short calculation involving use of (6), 
(7), (9), (10), (12), and (32), it is found to be 


_ tt —ytp V Loy PL 4 
Dy A?) = = (Ao—v+p_V)[Ao+22+4 F+AA+T]+ 


Py xe. H|vr—p, p_VR+(s—p_S)=|]— 








+. a ; 
-T9 em —p-)% +=" e—p, pI ) - 


—(v—p_I a uta) logH—¥— aa 7 t ral: (33) 
where F is defined by (15) and 


9 


lr = 





[(3—y)(u—U)+(y—1)A{(a—A)s+(s—S)A}]— 


V u—a—U+A u-+-3a—U—3A 
7 U+A U—A , 


~ y+ 


On a wave front receding into steady flow, (33) is seen to reduce to (17). 


+ See note at the end of section I. 


The val 
a similar 
steady flo 
Ill. 2. 


33) beco 


With z = 
D 


L 


Therefor 
The solv 


where 


and @ sl 
line thr 
taken t 
where ¢ 

On t 
first or 


and so 


which 


Fine 


to the 
rence | 
that 2 


when 
The 
by (2 
Ill 
Mach 
reced 








it of 
h sul 


17 











ON WAVES OF FINITE AMPLITUDE IN DUCTS 277 


The variation of A = D_r/Dt on an advancing Mach line is governed by 
, similar equation, which reduces to (24) on a wave front advancing into 
steady flow. 

Ill. 2. By (31), As/S, Ar/R, Av/V, p,—1, and p_—1 are all O(e) and 
33) becomes 


D 7 T l > oy i 
= C = = o){Ao+ 22+ F+C ti¢ . 
= (Ac) — (Ac){ Ao + F+-O(e)}+ O(e) 


With z = 1/Ao, this becomes 

dz y- 
2 YA ELF +O) e+ 1+(1+A0)0(6)}+A0°O(). (34) 
Therefore, if Ac is O(1), z satisfies, to the first order, equation (18) of Part I. 


[he solution is given by equation (23): 


K(x)z(x,t) = Ky%+J1 (x), (35) 
where K(x) = M*(1—M)-?, 
y+] . K dx 
I(x) =? aD 
="a | a4—0) 


nd a suffix 0 denotes values at some reference point on the receding Mach 
ine through the point (x,t). In the following the reference point will be 
taken to be the point (x,,t,) where this Mach line passes the exit and 
where z, is assumed to be prescribed (cf. section IT). 

On the other hand, if Ac is itself only O(e), then (34) becomes, to the 
first order, ov 1m, 

D_z/Dt (y+1)[2X+ F]z/(4A), 

and so Kz = constant on a receding Mach line, (35 a) 
which is the special case of (35) obtained when Kyz) > I. 

Finally, if Ao is large, (33) | 


D_(Aa)/ Dt (y+1)(Ac)?/(4a), 


yecomes 


to the first order. It will be shown in the following section that the occur- 
rence of a singularity of Ac does not imply a breakdown of the assumption 
that Aa/A O(e), and hence, to the first order, 


D_2/Dt (y+1)/(4A), 


when z is small, and the solution is, again, (35). 

The fluid acceleration and the local rate of change of pressure are given 
by (21) and (22), respectively, to the first order. 

III. 3. It is useful to consider briefly the differential geometry of the 
Mach line field. Let dr be the time interval between two neighbouring 
receding Mach lines, measured along an advancing Mach line; and dr’ the 
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time interval between two neighbouring advancing Mach lines, measure 
along a receding Mach line. Consider two points, Q and Q’, at a smal 
distance apart in the (2, t)-plane (Fig. 2; this figure has been chosen t 
represent the case of supersonic flow). Approximately, 

















dx = x(Q’)—2(Q) 
= (u—a) dr’+[u+a+{D_(u+a)/Dt} dr'l[dr+{D_(dr)/Dt} dr’) 
= (u+a) dr+[u—a+{D,(u—a)/ Dt} dr][dr’ +{D,(dr’)/Dt} dr], 
and dt = t(Q’)—t(Q) = dr+dr’+{D,(dr’)/Dt} dr 
= dr'+dr+{D_(dr)/Dbt} dr’, 
whence to the first order, 





(2a/dr)D_(dr)/Dt = D,(u—a)/Dt— D_(u+-a)/Dt. (36 
When Ao (but not AA) is large, this becomes, to the first order, 
1D(dr) yt], oa De 
a Dt ~~ a” ~te i’ 
ice. D_(Ao dr)/Dt = 0, (37 


by (6), (9), (10), and (33), since w and a are necessarily bounded. If Acis 
replaced by o, equation (37) characterizes a simple, receding wave in the 
one-dimensional, unsteady flow of an inviscid gas, and a point wher 
1/o = 0 is called a limit point.f The same definition will be adopted for 
the flow studied here, and since 1/Z + 0, (37) shows that the variation 
of o near a limit point is, to the first order, independent of the variation 
of the cross-section of the duct. Singularities of the limit type in one 
dimensional, unsteady flow, and in the analogous problem of steady 
two-dimensional, supersonic flow, have been studied by Craggs (4), Meyer 
(5, 6), and Stocker (6), and extensive use will be made in sections V and V. | 
of the qualitative results obtained by these authors. 


t Strictly speaking, a ‘receding’ limit point. A point where 1/A = 0 is an ‘advancing 
limit point, but such a limit point does not occur in a receding wave of the type studied her 
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Asingularity of o implies an infinite gradient of a in the advancing Mach 
direction (cf. list of symbols), just as in a simple receding wave. But since 
\is bounded, there is no comparably strong gradient of a in the receding 
Mach direction (in a simple receding wave, A = 0 and a is constant on 
receding Mach lines); and hence, if a is small on a receding Mach line, at a 
distance from a limit line where ca is still O(1), then a will be small of 
the same order on this Mach line, near the limit point. 

Equation (36) may be integrated, to the first order. By (6), (7), (8), 
and (31), 

D(u—a)/Dt D.(U —A)/Dt—(y+1)Ao+ O(e) 

k(y+1)Ao+[(U+A)/(U—A)|D_(U—A)/Dt+ O(e), 
and by (6), (10), (31), and (32), 
D_(u+a)/Dt—D_(U+A)/Dt = O(e), 


whence (36) becomes 


Die) = Gg et aa {p ma Di OA) GO +A + O10 
(38) 
1 (U+A D_,,, eo. 1 DM 
- salP ) i ube ot =e hy 
and by (35) and (8), 
D_(Kz)/Dt (U—A) dI/dx = —(y+1)K/(4A), 
and so from (38), to the first order, 
[(M—1)Kz]-! dr = constant on receding Mach lines. (39) 


III. 4. The direct implications of equation (35) have been discussed in 
detail in sections 3 and 3.1 of Part I. It may be worth repeating that 
a limit point (2 = 0) occurs upstream of the exit on all the receding Mach 
lines on which z is positive at the exit, and only on those. The position, 
r,, of the limit point is given by 

I(x,)+Koz = 0. (40) 
Thus, to each interval of time during which z, is positive there corresponds 
a limit line, and the limit point on it nearest to the exit occurs on the Mach 
line on which z, takes its least value, i.e. on which o takes its greatest 
value at the exit. It is also the limit point occurring at the earliest time, 
for the slope of the limit line is dt/dz = 1/(u—a) + 0 and it is a continuous 
line if z, is continuous.+ This limit point is of particular interest since 
the Rankine—Hugoniot equations for a weak shock suggest that it starts 
at the first limit point to occur in the flow (5). 

+ If z, is not continuous, edges of regression intervene (cf. 6), but the statement con- 
cerning the limit point occurring at the earliest time remains true. 
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Equation (35) applies also to receding waves travelling into steady 
supersonic flow in the convergent part of a duct. For the choice of the 
reference points in that case, cf. section II. 

This completes the proof of Part I of the instability of all steady shock. 
free flows which attain sonic velocity in a throat but are not nozzle flows 
For (31) holds for arbitrary ¢ in a sufficiently narrow strip of receding Mach 
lines adjacent to any wave front receding into steady flow, since Aw, Aq 
and AA vanish on such a wave front. Therefore, only steady flows fo 
which all wave fronts pass out of the duct after a finite time will be con. 
sidered hereafter. 

The first approximation to the flow in a receding wave where the assump. 
tions (31) hold may therefore be summed up as follows. The velocity and 
the speed of sound do not differ from their local values in the steady flow: 
AA = D_(r—)/Dt is zero, but Ao = D,(s—S)/Dt is given by (35), wher 
z = 1/Ao. A limit line may occur and its approximate position is given by 
(40), in spite of the fact that this approximation does not distinguish between 
the Mach lines of the disturbed flow and of the steady flow. In particular 
the approximate position and time of the first limit point is found, which 
may be interpreted as the position and time at which a shock wave starts, 

Similar results are obtained for an advancing wave where 


As/S = O(e), Ar/R = O(e), Ao/X& = Ole) (€ <1). (41 
In that case Ac is zero, in the first approximation, and AX = 1/y is given by 
L(x)y(x, t) = Lo yy—J (2), 


where L(x) = M*/(1+-M)-?, | 
4 4 (42 
, eyrif - L dx * 

rae J A(A+U) 


Xo J 





and a suffix 0 denotes values taken at some reference point on the advancing 
Mach line through the point (x,t); for the choice of the reference point, 
ef. section II. The fluid acceleration and the local rate of change of pressure 
are given by (28) and (29), respectively, to the first order. 

III. 5. From a mathematical point of view, the assumption made in (31) 
that AA/A is O(e), at most, is not altogether essential (cf. section IV. 1). 
It is important, however, for the physical understanding of the problem 
to note that AA is indeed small. To see this, consider first equation (35a). 
In a duct of constant cross-section, a zero of o at the exit would imply a 
receding branch line (6), which is not to be confused with the advancing 
Mach lines, all of which would be branch lines since the receding wave 
would be a simple wave. In a duct of varying cross-section, a zero of Ac 
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4 the exit implies at least a narrow strip of receding Mach lines on which 
35a) holds. Such a strip might be called an ‘acoustic’ strip, for in such a 
trip, and only in it, the present theory reproduces the predictions of 
Acoustic, or Linearized, theory. 

Now consider the advancing Mach lines. An equation for AA holds on 
them which is analogous to equation (33) for Ac, and the solution is given 
by (42), except that L(x) has to be replaced by 

y+1 f Ac 
L L(x)exp |? a Be dt\, 
\ 4 AJ 
since Ac cannot be assumed to be small throughout the receding wave; the 
ntegral is to be evaluated along an advancing Mach line, and it remains 
bounded, by (37), even if the path of integration crosses a receding limit 
ine. But AA = 0 on the wave front receding into steady flow, and so the 
solution reduces to 
L'y = constant on advancing Mach lines 
by analogy with (35a), not only on certain narrow strips, but throughout the 
wave. In contrast to the primary receding wave, the secondary advancing 
wave is therefore altogether an ‘acoustic’ wave, and in the first approxima- 
tion, where terms O(e) are neglected, AA 0. 
IV. The second approximation 
By (9), (31), and (32), 
D,(As)/Dt = Ao, D(Ar)/Dt = —v+p,V = O(e), 
whence 
As | Ao dr | 
J 
; (43) 


A? of ( € ir) | 


0 J 





where the integrals are taken along an advancing Mach line, with t measured 
from the wave front travelling into steady flow. A knowledge of Ao from 
the first approximation therefore permits us to calculate, as a second approxi- 
mation to the flow, the first-order correction to the velocity and to the 
speed of sound. 

It is convenient to choose as characteristic parameter distinguishing 
receding Mach lines the time, f), at the exit measured from the wave front. 
By (39), equation (43) may be written 

[ (M—1)K dry 


se 44 
(My l )Ko <9 dty . 
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From Fig. 2, nee ek , 


dx = (u+a) dr+(u—a) dr’ (45 
and if both Q and Q’ are situated at the exit, 
(Or/0t),.-2, = AIto/dty = (Ag—U,)/(2A,), (46) 
to the first order, whence the velocity correction As becomes 
to 
wall 5 F 
As(x,t) = OK. (M—1)KAoy dt); (47 
2Ko 


here t, = t— { dx/(U—A) is the time at which the receding Mach line 


Zo 
through (x,t) passes the exit, and (M—1)K is to be taken at the point of 
intersection of the advancing Mach line through (z,t) with the receding 
Mach line passing the exit at time ¢, i.e. at the position 2’ given by 


x 


[ A dx 


A2—U?2 = 3(to—to). 


Fo 
Since the first approximation does not distinguish between the Mach lines 
of the unsteady and of the steady flow, it is obvious that the interval is 
O(t)) over which x varies on the advancing Mach line between (2, t) and the 
wave front. If t) is small and of the same order as e, at most, (M—1)K 
may therefore be taken at (x,t), without loss of accuracy, and (47) becomes 
to 


[ Aa, dt). (47a) 


0 


As(x,t) = an |K(z) 


For an advancing wave travelling into subsonic, steady flow, 


As = O(er’), 
to 
Ar(x,t) = _ [ (r+ 1) La, dss, (48) 


lx 
rhere 3, = t— [ ii 
where . } Usa 


is the time at which the advancing Mach line through (z,t) passes the 
entry of the duct, and (/—1)L is to be taken at the position x” given by 


r 


A dv 
| A2—U? 


r 


(A suffix 0 is here used to denote values taken at the entry of the duct.) 


= 3 (9) —04). 
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IV. 1. The equation (33) can be integrated formally even if the assump- 
tion that AA may be neglected is abandoned. The integral is still given 
by (35), except that K has to be replaced throughout by 


t 

(y+! | Ay) 

(4 JA J 

where the integral is taken along the receding Mach line through the 

point (x, t). On the other hand, if AA is not neglected, 
D_(u+a)/Dit—D_(U-4 A) Dt = $(y+1)AA+ Ofe), 


and it follows that (39), and hence also (47), remain valid as they stand. 
Similarly, if the assumption Ao = O(e) is abandoned in (41), LZ has to be 


K’ K exp 


replaced by 3 


L’ = Lexp{¥+* [ > as 
| 4 J A 
throughout (42), but (48) remains valid as it stands. 

Equations (47) and (48) can therefore be employed to obtain the second 
approximation even for a region of interaction of two primary waves, where 
Ao and AA must be expected to be of the same order of magnitude. For 
instance, if a disturbance is set up at the entry of a duct, in a steady flow 
that is supersonic near the entry, Ac, and AA, may be assumed to be pre- 
scribed as functions of the time at the entry (cf. section II; the suffix 0 
is again used to denote values taken at the entry of the duct). Then 

0 
Ar(e,t) =< | (M+1)LAA, d9,+ Arty to) 


2 Ly « 
to 


(49) 


t 


Aa(e,t) = ——— | (M—1)KAog dt, +-As( to, 9) 


A 


- 0 « 


0 
at any point (a, ¢) in the region of interaction (Fig. 1), with to, 3), 2’, and 2” 
calculated as for (47) and (48). (Note also that 7, in (43), must still be 
measured from the wave front receding into steady flow, and hence, for 
the calculation of As in the receding wave, after its emergence from the 
region of interaction, the lower limit of the integrals in (47) and (47 a) must 
be replaced by the time, 7), during which the disturbance lasts.) 

It should be stressed, however, that the formulae (49) represent only a 
partial solution of the problem, because the solution given here for the first 
approximation in the region of interaction is only formal. It does not give 
the limit lines that might occur in the interaction region, nor the distribu- 
tions of Ac and Ad, respectively, on the Mach lines bordering the inter- 
action region, which would be required as boundary conditions for the 














284 R. E. MEYER 


determination of the limit lines that may occur in the primary waves after 
their emergence from the interaction. 

IV. 2. It is now easy to see for which types of waves the theory yields a 
consistent approximation. Consider first a receding wave generated at 
the exit of the duct. Equation (47) may be written 


As = —}(M—1)KF(t,))+4 [ F(t) d[(M—1)K], 
to 


where F(t)) = _ | Ao, dt, 
0 


0 « 


and the integral in the second term is taken along the advancing Mach 
line, from the wave front to the point where As is to be calculated. For 
a duct of finite length (without sonic throat) the condition 


F(t))|< Ay for 0<4<G% 


will therefore ensure the Asy/Sy < 1 for tf) < J. Note that the boundary 
conditions assumed in this investigation do not determine the primary 
wave beyond the advancing Mach line meeting the wave front at the 
entry of the duct. The time at which this Mach line passes the exit there- 
fore represents an upper limit for 7). If ¢ and e’ denote respectively the 


to 


largest values taken by |F(t,)| and | ( F(t)) dt)|, the errors in the formulae 
! 0 ' 

of the first approximation will be O(eT), at most, where T is the maximum 
time a receding Mach line spends in the duct; the errors in the formulae 
of the second approximation will be O(e'T), and AA will be of the same 
order, at most. 

For a primary wave advancing into steady subsonic flow, the corre- 
sponding condition is, by (48), 

G(d))| <A, for 0<3%,< JT, 
do 
where G(9,) = + | Ary dé%. 
Ly. 
0 

No limit on J) is implied if the flow is subsonic throughout the duct. 
For primary waves travelling into steady, supersonic flow, a consistent 
approximation will be ensured by the same two sets of conditions, by (47), 
(48), and (49). 

On the other hand, the theory will be of interest only if |Ac| (or |AA)) 
itself does not remain small at the exit or entry so that shock waves may 
occur inside the duct, which could not be predicted by Acoustic theory. 
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Therefore, unless application be restricted to the neighbourhood of wave 
fronts and to waves generated by disturbances of short duration and extent, 
the requirement is that the disturbance should be a rapidly fluctuating one. 
The mean squares of the differences of velocity and pressure from their 
respective steady-flow values at the exit should be small for the approxi- 
mation to apply. But the respective mean squares of the fluid acceleration 
and the local rate of change of pressure, at the exit, should not be small 
for the theory to be of interest. 

This suggests that the theory may be applied to the study of the pressure 
fluctuations generated by turbulence in the high-speed flow in ducts. 
These pressure fluctuations are of very small amplitude, and the shock 
waves to which they give rise are of very small strength, but their study may 
be of interest if it is desired to distinguish between the velocity fluctuations 
directly due to local vortices and those due to pressure fluctuations propa- 
gated from other parts of the flow fieid. 

It appears from the theory described above that the smallest eddies will 
generate shock waves after the shortest time. These effects of turbulence, 
however, will not be followed up here, since it appears to the writer, on 
the grounds of a rough estimate, that in wind tunnels of the types commonly 
used today even the shock waves due to the smallest eddies will form only 
downstream of the working section. More marked effects should, however, 
be expected from the receding waves in ducts with an entirely subsonic, 
oran entirely supersonic, mean flow in which a near-sonic velocity is reached 
in a throat. 


VY. The shock path 

Knowledge of Ao from the first approximation permits us also to calculate 
the shock-path in the (a, f)-plane, to the first order in the shock strength. 
The example chosen is again that of a receding wave, generated at the exit 
of a duct. 

At any point of the shock-path two receding Mach lines meet. Let a 
suffix s refer to the shock-point, and the suffixes 1, 2 refer respectively to 
the upstream and downstream receding Mach lines through the shock-point; 
in particular, let ¢, and ¢, denote the respective times at which these Mach 
lines pass the exit. 

[t is a well-known consequence of the Rankine—Hugoniot shock equations 
that, to the first order in the shock strength, the speed of propagation of a 
shock is the arithmetic mean of the speeds of propagation of signals up- 
stream and downstream of the shock. More precisely, the path of a receding 
shock bisects, at every point, the angle between the receding Mach lines 
at the point (and similarly for an advancing shock and the advancing 
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Mach lines). Hence, to the first order, AB = BC in Fig. 3, and s0, if t 
denotes the time at which a receding Mach line passes the exit, then at the 
shock (Fig. 3), 


{(@x/ét),, dtp/dt,',+-{(0x/ét), dt,/dt,.t, = 0, (50) 


where the suffix ¢, denotes that the derivative is taken at constant time. 
t = t,. 











Fic. 3. 


Another well-known property of a shock is that, even to the second order 
in the strength, there is no change of entropy ecross it. This means that 
the shock-path fits into the solution of the equations of isentropic flow in 
the following sense. The isentropic solution possesses a limit line with a 
cusp at the point where the shock-path starts. Between the two branches 
of the limit line, the mapping of the characteristic plane on to the (2, t)-plane 
is three-sheeted (5, 6). The shock-path lies in this region of the (2, t)-plane, 
and the flow passes across the shock from the first to the third sheet. 
Suppose the shock-path to be marked on both the first and the third sheet. 
Then, asa linet = constant = t, is followed on the first sheet from the shock 
to the limit line, then on the second sheet back to the other branch of the 
limit line, and finally on the third sheet from the limit line to the shock, 
the total change in x must be zero. Points on this linet = t, may be regarded 
as different if they lie on different sheets, even if they seem to coincide in 
the (x, t)-plane, and there is then a (1, 1) correspondence between the points 
and the receding Mach lines crossing this line, and thereby also between 
these points and the respective times at which receding Mach lines pass 
the exit. The total change in 2 may therefore be written 

te 
| (@a/@ty),, dty = 0. (51) 
th 

+ The possibility of the occurrence of an edge of regression is ignored for the moment, 
and that case is discussed in section V. 1 below. 
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By (39), (45), (46), and (35), 
(@x/éty)). a(M—1)Kz/(Kyz) = —a(M—1)[1+1/(Kg2,)], (52) 


where all the quantities on the right-hand side are taken on the same 
receding Mach line; those with a suffix 0 at the exit, and those without a 
suffix at time ¢,. If x, and 2, denote respectively the position of the limit 
points on the line ¢ t,, and @ the greatest value of |éx/ét,| on this line, 


between the limit points, then |2,—2,| = Of{@(t,—t,)}, which is small 


‘ 
compared with 6. To the first order, therefore, quantities denoted by 
capitals without a suffix, in (52), may be regarded as constants in the 
integrand of (51). Moreover, to the first order, Aa may be neglected in 
comparison with A. It then follows from (51) and (52) that 

te 


. Fane. (53) 


I(x,) (te —t,) Ko . 
l 


1 
i.e. the position of the shock-point where the receding Mach lines that pass 
the exit at times ¢, and ¢, respectively meet, is given directly by the mean 
value of Ac at the exit between ¢, and f,. 

Elimination of J from (52), by the help of (53), substitution in (50), and 
integration gives 


t tb j 


|} dt | Ao, dty—d(te I») | “t ) Ao, di, = 0 (54) 
fy CO Cc C 

where C is arbitrary. From that relation, one of f,, t, can be determined 
when the other has been chosen. 


The strength of the shock is given by 
re) (Po—p,)/P y(8.—8}) A+ O0(8?) 


where all quantities are taken at x = 2,, t = t,), as can be shown without 
difficulty from the shock-equations (7). Since S depends only on x, we may 
evaluate the jump of As across the shock, instead of the jump of s. The 
respective shock-points on the first and third sheet do not lie on the same 
advancing Mach line, but the rate of change of As in the receding Mach 
direction is O(e) by (10) and (31), and hence, to the first order at least, 

2 

fe) (y A) Ao(dr dt.) dty, 

th 
where all the quantities in the integrand are taken on the same receding 
Mach line; those with a suffix 0 are taken at the exit, those without a suffix 
on the advancing Mach line through either of the two shock-points. The 
quantities depending only on the value of x on this advancing Mach line 
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may again be regarded as constants during the integration, whence, }y 
(39) and (45), 


2 


M—1)K [ 
i i [ Aap dtp. (5 





mo 
| 


th 
The absolute speed of propagation of the shock is 
dx,/dt, = 4(u,—a,), +4(u,—a,)o, 
by the property of the arithmetic mean referred to at the beginning of this 
section. By (43), (44), and (47a), this may be written 


2 


th 
dx, . yti K ‘ 
=U —A4o it — 5 | +  ) deo to (56 

0 ; 








dt, 0 
0 
where the use of (47a) implies the assumption that ft, and t, are small, 
For an advancing wave, the formulae corresponding to (53) to (56) 
respectively, are 








] 
—_— = Adv dt 
I(x) t, Ly aI — 
ho oy 
| dt | Arg ned | + | } ar, dt, = 9, 
i C C c 
. y(M+1)L 
§= ——— Ady dt 
SAL, | ; 
1 1 L{f, f 
dx 
and *— U4A4%°-(mM41)2 ' 
( dt, ro 3 (M T. + Ary dt, 
0 0 
where ¢,, f, are the times at which the advancing Mach lines meeting the 
shock at x = a,, t = t, on the upstream and downstream side respectively, 


pass the line x = 2, on which Ad, is prescribed (cf. section IT). 

V. 1. It should be noted that equations (53) and (54) are formally 
identical with the equations first obtained by Whitham (2) for the shape 
of the shock in his first-order theory of the steady, supersonic, axially 
symmetrical flow past a slender body of revolution. The simple geometrical 
interpretation given by Whitham may therefore also be applied to equa- 
tions (53) and (54). This suggests that similar results may be obtainable 
also for a number of other problems. The arguments used above to establish 


(50) and (51) apply to any weak shock separating two regions where the 
equations of the flow are hyperbolic. To the first order of the shock strength, 
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the shock equations may therefore be written in the form (50), (51) in all 
such cases. Their integration in the form (53), (54) would then appear to 
depend essentially on whether the integrand of (51) may be written in the 


jorm 


,(X,, ty by (ty) +-ho(2,, ty )o(to), 


where f, is a characteristic parameter. For, to the first order, the variation 
f the d; (2 1,2) between the two branches of the limit line can always 
be neglected in the integrand of (51). If we write, for brevity, 


[ff]? for f(t.)—f(t), {f}} for f(t.)+f(4), 


and ub, dt, = F(t), 
t follows from (51) that 
$F li+¢o[Ali = 0, 


which has the same form as (53); and using this to eliminate the ¢,; from 


50), we obtain 


[Ali tdh 


:— [Ai amt = dC AR RD) —21F a = 0 


which gives (54). 

It should also be noted, however, that for the problem studied in this 
paper the advantage of writing the equations for the shock in the general 
form (53) to (56) lies in their covering a number of different cases rather 
than in the degree of accuracy to which these equations and their geo- 
metrical interpretation hold in any one case. This is best explained at the 
nstance of two examples. 

Suppose first that Ao, is a continuous function of ¢,, and analytic for 

0. As explained in section III. 4, the point where the shock-path 
starts lies on the receding Mach line passing the exit at the time, tj say, 
when Ac, reaches its highest value. Therefore, for small values of 
t* 


0 0° , 


Ao* + 1 Bt? : Lnt;? +... and B< 0, (57) 


( 6 


Ao 


0 


When substituting this into (53) to (56), we must remember that all but the 
largest term was neglected in the integrand of (51). It follows that (54) holds 
nly in so far as the largest term on the left-hand side is concerned, and 
with Ac, given by (57), this implies that 

(B/24)(t,+4)(4—4)? = 0, 


irom which we conclude (ef. Fig. 3) that 


t t+ O(t?). (58) 


5092.19 
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This shows that our approximation does not permit us to take account of 
terms O(f)? in Aoy; but if terms O(t3) are neglected, it is immediate) 
evident from Fig. 3, for reasons of symmetry, that (58) follows from the 
assumption of (57). By (57), (58), (53), (40), and (55), 


I(x,)—I(a*) = BI (x*)t?/(6Ky)+ O(t?), (59 
§ = y(l1—M)KAo*t,/(A K,)+ O(t2), (60 


where «* is the position at which the shock-path starts. By (56), (57), (58 
(35), and (40), the time at which the shock reaches position x, is, to the 
first order, 

_ th 

.) = tar)+ (4 1—1(@,)/T(a*)] [ (Aop/Ac*) dt (6 

t(x.) = t,(2*) + qn [ — I (x,)/L[ (a )] (Ao,/Ao*) dt). (61 
rs 0 

Secondly, suppose that Ac, is discontinuous at tg = 0 and, say, 


Ao, =0 for t,< 90, | 


0 


; : ; 62 
Ac, = Ac*-tot,+362+... and «<0, for t,>0/J \ 


Then the shock starts on the wave front, and one of the two branches of th: 
limit line is replaced by an edge of regression (6), at which the first and 
second sheets of the mapping meet. The arguments employed in th 
preceding section to derive equations (50) to (56) remain valid, however 
if the words ‘limit line’ and ‘limit point’ are replaced respectively by ‘edg 
of regression’ and ‘point of regression’, wherever appropriate, and this 
case is therefore also covered by (53) to (56). On substituting (62), re- 
membering that only the largest term in each expression is reliable, we find 
tAo*t, t, = 0, 
from (54), since now ¢, < 0 (Fig. 3), and 
K,/1(a*)—K,/1(x,) = 4a(t, +1), 


from (53). Since the shock travels in the direction of 2 decreasing, x, < 2* 
and hence, by (35) and since «a < 0, (t;+t)) > 0; we conclude that 


t, = OfG) (63 

at most, and I(x,)—I(a*) = al (ax*)ty/(2K,)+O(8); (64 
also, from (55), (56), and (62) to (64), we have 

8 = y(1—M)KAo*t,/(2A K,)+ O(#), (65 

tay) = ter) + [ 4+ Aot[1—104,)/T@*)B/(4a2). (66 


“quations (63) to (66) should be compared with (58) to (61). 
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SUPERSONIC THEORY OF DOWNWASH FIELDS 
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SUMMARY 

Two methods have so far been proposed for the calculation of downwash fields 
behind a supersonic wing on the linearized theory. That due to Lagerstrom (1) is a) 
extended conical flow method, while Mirels and Haefeli (2) and Robinson and 
Hunter-Tod (3) replace the wing by an equivalent vortex distribution. Sine, 
references (1) and (2) are not readily available in this country, their contents an 
briefly summarized. 

We then show that the basic formula of Robinson and Hunter-Tod (3), giving thy 
downwash in terms of the flow on the wing, can be deduced very rapidly fron 
Hadamard’s theorem. When the flow is conical, this equation can be shown to by 
identical with that of Lagerstrom. Also, this formula makes clear the nature of th 
approximation made by Mirels and Haefeli. Finally, an asymptotic expansion of th 
integral kernel occurring in the equation gives a ‘practical’ method of calculating 
the downwash behind any wing. 


1. Introduction: Lagerstrom’s method 

THE determination of the flow field behind a wing is of considerable 
practical importance, and also of some mathematical interest. If we take 
a rectangular coordinate system such that x is the direction of the free 
stream, while the wing lies in the plane z = 0, then the components v (= 4,) 
and w (= q,) of the velocity q are called the sidewash and downwash 
respectively. 

Calculations of supersonic downwash have so far been made using the 
linearized theory. This is unfortunate, since two important effects, the 
rolling-up and downward deflexion of the trailing vortex wake, are then 
missed. However, even using the linearized theory, the formulae obtained 
are rather involved, and work with more accurate equations will lead to 
even greater complications. Also, since the linearized theory gives good 
results for thin wings at small angles of incidence, it would seem that 
viscous effects should be considered before those of non-linearity. 

Throughout this paper we use the rectangular wing-tip for purposes ol 
illustration, but the arguments are obviously general. 


Within the Mach cone OAB, the velocities are known functions of the 


complex conical coordinate b, with O as origin. This coordinate is defined by 


Sue Butte ba 


x x r 


cos @—i(1—r?)!sin 6 


(Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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where 8 = (M?—1). Throughout this paper the quantities which we 
lescribe as velocities are, in fact, the ratios of the physical perturbation 
velocities to the free-stream velocity. Also, since the problem is anti- 
ammetrie in z, we shall only consider the region z > 0. In particular, 


we shall have ‘ 


u = = L(b), (1) 


Ca 


being the perturbation velocity potential (we use the symbol L since this 





\ mY 
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‘\ 
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CANCELLING /“ WING.’ a . 
/ / f f 








juantity determines the lift of the wing). In the case of rectangular wing-tip 
it incidence «, 9 
aX . 
Lib) qre(sin- IVb), 


his formula cannot be valid downstream of the Mach wave from the 
trailing edge, since it gives u + 0 on the trailing wake. Lagerstrom pointed 
out that, by the superposition of certain elementary conical flows, this 
excess lift could be cancelled, and that the actual downwash field would 
then be the sum of the flows due to the physical wing and the ‘cancelling 
wings , 


The conical flows treated in normal wing theory give constant downwash 
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over the wing, but Lagerstrom’s cancelling wing is one of constant lift. |p 
the notation of Goldstein and Ward (4) it is specified by 


re(G) = tis (t <b), 
=0 (¢>6). 


Using the methods of reference (4), it can be shown that the downwash due 
to this wing is given by 


U- . 
Ba, = re{G(b, t)}, 
21 | 7 
where G(b,t) = -{5m T\+2tan7?7' + =F! er — In 9 
r — i=v —?. BR = 1—, (1-6?) 
t b 


In the plane of the wing, t becomes a real coordinate, and 


Fins = G(b, t). 


Bis 
Lagerstrom then finds that the downwash field of the wing-tip in the plane 
z = 0 may be written 
C L 
1 db, (3 
b 


c 


1 
w = B | {G(b,)—G(b,t*) 
0 


e By/c - b 3a 


z/c—] 


where 


(A similar formula may be written down for a general station.) Lagerstron 
also shows that the flow in the Trefftz plane (a — 00) is particularly simple 
he has published an account of this work (5). 


2. Vortex-distribution methods 

Any lifting surface can be replaced by a system of bound vortices. Sine 
it also trails a vortex wake behind it, the calculation of the resultant down- 
wash is equivalent to finding the flow due to a given vortex system. Using 
this equivalence, Robinson and Hunter-Tod (3) deduce the formula whic! 
we derive in section 3 from Hadamard’s theorem. These formulae are just 
as complicated as Lagerstrom’s. 

Mirels and Haefeli (2) therefore replace the lifting surface of the wing b} 
a lifting line, and calculate the downwash due to this. Their formulae ar 
therefore approximate, but the error is not important for stations more thal 
a chord downstream of the trailing edge. The equations are much mor 
tractable than those of Lagerstrom, so that their paper represents a ver} 
definite contribution to the problem of finding downwash fields. 
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The strength of Mirels and Haefeli’s line is variable, being equal to ['(y), 
the local circulation round the wing. Since the perturbation velocity 
potential ¢ is antisymmetric in z, and is zero on /, the leading edge of the 
wing, [' is twice ¢,, t denoting the trailing edge. Logically the lifting line 
should be the locus of local centres of pressure, but this gives a non-linear 
curve, Which is nearly as difficult to handle as a lifting surface. The line is 
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therefore taken to be straight, and generally passes through the centre of 
pressure of the wing. For representing a triangular or swept-back wing, a 
lifting line with a single bend is advantageous, but for a rectangular wing-tip 
Mirels and Haefeli take the line to be straight, and normal to the direction 
of the free stream. 

Now vorticity must be conserved, so that the lifting line will trail a vortex 
wake behind it in the manner indicated in Fig. 2; the vorticity shed from 
the main line between y and y+-dy is equivalent to a line vortex of constant 
strength (el /éy)dy. Using the methods of Robinson (6), Mirels and Haefeli 
show that the downwash field due to this vortex system is 


Ye 


w(r) : K(0) : S de, (4) 
27. cy 
K(x) (x—2' )(y—y )(. R?— B?2*) 7 (5) ‘ 


Ri(x—2' P—B22Y(y—y' P+2} 
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where R? = (x—a')?—B(y—y’)?+-2%}, 


y, and y, are the coordinates of A and B, and we have taken the y-axis to 
coincide with the lifting line; eI’/éy is easily calculated for most wings of 
practical importance: for the rectangular wing-tip at incidence a 


F--2 C0) 


The computations of Mirels and Haefeli were made with formula (4), and 
were in consequence rather complicated. However, we should not expect 
the formula to hold for small x, since in obtaining it we have distorted the 
lift distribution of the wing (this will be made clear in the next section), 
This suggests that we should consider an asymptotic expansion of (4) in 
powers of x. Now 


Biy—y') B(y—y'K(y—y'P 4. Ole 6) (6) 


Sat 





9 


K(0)~ —— 
(y—y )°+2z 





y—-y 


2x* 
so that the required expansion of (4) is 


Ye Ye 

lf yy ,, - 

: a 3 

Zr J (y—y )*-+2* cy 4rrx* | 
Y1 "1 


, ol , = 
u ied — (y—y’) — dy’+ O(a-4). (7 
cy 

We shall see later that the subsequent terms of this expansion are not 
reliable. The first term represents the flow in the Trefftz plane, and was 
discovered by Mirels and Haefeli using a different method; we shall now 
show that this term does indeed give the same flow as Lagerstrom’s method. 

Xepeating the above argument, we can show that the sidewash in the 
Trefftz plane is given by 

i £ z a 
| TA (3 
2a J (y—y )*+2" oy 
V1 

In both (7) and (8) the integral can be taken from —oo to +00 if we take 
el /éy = 0 if y < y, or > ¥ (this is consistent with the definition of I’). 
Combining these two equations, we can write 

1 f F(%) de 


271 C—C 


where the complex variables ¢, ¢’ are defined by 

C= yt+2z, C’ = y’+¢2’'. 
The contour C consists of the real axis from — R to + R, together with the 
semicircle Re‘? (0 < 6 < z) as R-> ; the real axis must be cut between 
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y, and y, before the integral can be defined: the path goes along the upper 


lip of the cut. Thus, v—iw = F(t). (9) 


Now F(¢’) is an analytic function whose real part is equal to éI’/éy’ on 
y’ = 0 and which behaves suitably at infinity (this is necessary to make the 
ntegral round the semicircle converge). (9) is therefore a statement of 
Lagerstrom’s procedure (5). 

The second term is simple to calculate, but the important point is that 
it is independent of z; hence if we know the flow in the Trefftz plane and 
in the plane of the wing, it can at once be found wherever this expansion 
is valid. This holds for stations not too near the Mach wave of the trailing 


edge. 


3. Application of Hadamard’s theorem 
Hadamard (7) has shown that if the sources of a disturbance are confined 


to the plane z = 0, then the appropriat solution of the linearized equation 
wail “h Ch od 
’ Ca" cy” oz? 
Q2*°r ~ 

. o(r) = —— || d(a’.y’,0) 5, de'dy’, (10) 
we J. i Rs f 


R being defined by equation (5). The symbol */ | indicates that the finite 
part of the integral is to be taken; the region of integration is that part of 
the plane z = 0 lying within the forecone of the point r. Integrating (10) by 
parts, Ward (8, 9) has obtained the formula 


d(r) l od 2(x—2x ) dx'dy’, (1 1) 
aJ J ox’ {(y—y P+e2yR ; 

¢ ¢x is zero except on the wing, so that the region of integration is that 
part of the wing lying within the forecone of r. He suggests (11) as being 
suitable for downwash calculations, but we do not regard the formula as 
being well adapted to this purpose. It does not give w directly, and when 
lifferentiated with respect to z, it gives rise to another divergent integral, 
which must be treated by the finite part technique. 


If we integrate (11) by parts again, we find 


lf a o*h zR P P 
d(r) | . tan if . . | dy dy 
7 Cx CY \(a- 4 yy Yy )} 
d rr gf 
and hence ees : je “? . K(x’) dx'dy’, (12) 
Cz 7) J cxcy 


where A(x) is defined by formula (5). Again, the integral is over that part 
of the wing which lies within the forecone of r. This integral is finite, which 
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is a practical convenience. This formula is equivalent to Robinson and 
Hunter-Tod’s formulae (30), (33). 

We shall now show that Lagerstrom’s formula (3) is a particular case of 
(12) in which the flow on the wing is conical. If we confine ourselves to the 
plane of the wing, (12) becomes 

1 ff ed ./f(x—a’)2?—B%y—y') .,,, 
Ww a | tee: ad ; PAY Y's ax'dy’. (12a) 
Ox Cy (x—ax \(y—y 


7 


Since the flow on the wing is conical, 


 _ 1H) > Bel 


Cx ox’ oy’ x’ cb 


where b = By’/x’ (identical with Lagerstrom’s variable). If we also use 


Lagerstrom’s other conical coordinate t = By/x (12a) becomes, for a rect- 


angular wing-tip 





u B 1(b,1) 2 db (13) 
0 
cia t)2 } - 
where I(b. t) =. | vil ¢) (t—b€)*} dé 
a (1—&)(t—bé) : 
0 


(é «’/x). It is easy to show that 
I(b, t) G(b, t)— G(b, t*), 
so that (13) is identical with Lagerstrom’s formula (3). Also (12) will, in 
principle, give the downwash field of any wing exactly (subject to the 
limitations of the linearized theory) while Lagerstrom’s method is restricted 
to wings on which the flow is conical. However, on other wings, the 
integrals are likely to be rather troublesome. 
To handle (12) in general, some approximation must be made, and the 
most obvious one is to set 
h(a’, y’) b(y')H(x’), (14) 
where ¢,(y’) denotes the value of ¢ on the trailing edge and H(z’) is the unit 
function , , 
; H(x’) 0 (x 0), 
: ce > 0}. 


This replaces the smooth increase of 4 from / to ¢ by a sudden jump across 
the y-axis. From (14) 


ch 
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by’) d(x’) IT (y’) d(x’). (15) 


where ['(y’) is the local circulation round the wing, and 8(x’), the Dirac 
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delta function, is the derivative of the unit function. Hence this approxima- 
tion is equivalent to replacing the wing by a lifting line along x’ = 0, the 
precise analogue of Mirels and Haefeli’s procedure. From (15) 


o*h ] ol . ’ 
. O(x ) 
COX OY 2 ¢ y 
ieee. ae 
so that (12) becomes w K(0)— dy’, 
= 7 @ 
27 . oy 


identical with Mirels and Haefeli’s formula (4). 

It is difficult to make a general estimate of the error of the approximation 
\3). But if we consider an asymptotic expansion of the exact formula (12), 
this can be compared with the expansion (7), derived from (4). We find 


y—y' B2(y—y’) Py y')2" 


K(2’) ~ ; 


) ) > 2 


(y y')* 2° ou" x 


B%(y—y'){4x'2+ B?(y—y')2+ 98%" 


Sat 
so that the asymptotic form of (12) is 
t 
1 ¢ ( y’ B2(y—y’ ve . 
terse = § {ey | ay 4 
7 J \(y—y’)?+2? 223 j * Cx CY 
span l 
t 
i. - , no 
p~ , , , ' - 
a (y—y’) dy 2 “$ - dx’ +O(a-4). (17) 
we | Ox’ Cy 
spar l 
t 
. r “4 , l or 
Now nt da (= — on 
J ox oy cy |i 2 cy 
i 


so that the first two terms of (17) are identical with those of (7), which is 
therefore accurate, in general, to O(a-%). If the axes are chosen so that 


(y—y’) a —— as 0. (18) 
CX C Y 
asen 
then these first two terms will be accurate to O(x-*); the line which satisfies 
18) will not usually pass through the centre of pressure. Now in their paper 
Mirels and Haefeli represent a ‘narrow’ triangular wing by an unbent lifting 
line; they get best agreement with the known downwash field when the line 
is placed at the ?-chord point, while the centre of pressure is at the 3-chord 
point. The above remark accounts for this phenomenon. 

We therefore claim that the expansion (7) will give the downwash field 
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behind a rectangular wing to the same accuracy as the much more compli- 
cated calculations of Mirels and Haefeli, while the expansion (17) should 
give better accuracy, since the term in 2~* is not taken into account in 
Mirels and Haefeli’s work. Calculations to check this claim are being started. 
and will be reported elsewhere. 

In the above calculations, we have tacitly assumed that ¢¢/¢y is con- 
tinuous across the leading and trailing edges, so that they are strictly 
applicable to rectangular wings only. It is quite easy to take discontinuities 
into account, and we find that the first two terms of (17) are still accurate 
if cl’ /éy is taken to be twice the sidewash immediately behind the Mach 
wave from the trailing edge. Robinson and Hunter-Tod (3) have given 
explicitly the correction necessary for a discontinuity, so that we need not 
repeat the calculation. 
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a SUMMARY 
iven A source method is used to solve some problems of unsteady motion of a slender 
body through a perfect, inviscid, compressible fluid. The effects of variations in the 
not : . . . . . . 
forward velocity and in the angle of incidence are investigated on a linearized theory. 
The results are established for slender bodies both of revolution and of general cross- 
section. The extra aerodynamic force components produced by linear and angular 
acceleration are found to be small compared with the corresponding aerodynamic 
for force components in steady motion, provided that the changes in linear and angular 
ody velocity, in the time the body takes to travel its own length, are small compared with 
the actual linear and angular velocities respectively. A similar result is obtained for 
e effect of angular velocity. 
Introduction 
1 by A LINEARIZED potential theory is used to solve the problem of the motion 
aii of a slender body through a perfect, inviscid gas. First approximations 
per: to the aerodynamic forces are obtained for various prescribed, unsteady 
velocities of the body. In this way the effect of both linear and angular 
the acceleration on the aerodynamic force is investigated for oscillations of 
mort 
ins small amplitude. 


The work, which is developed by a source method, deals with the slender 
body of revolution in Part I, whilst the slender body of general cross- 
section is considered in Part II. In each case the effects of variable 
forward speed and incidence are investigated. 

For a slender body, the thickness ratio, e, defined as the ratio of the 
maximum diameter of the body to its length, and the angle between the 
tangent to the body surface and the direction of motion must both be small 
and of the same order. The condition that the curvature at any point on 
the contour of the body in a transverse plane of cross-section be at most 
O(1/e) excludes very small outwardly-convex radii of curvature. The nose 
of the body is taken to be pointed so that attached shocks only occur; the 
third-order entropy change across these shocks will not affect this linearized 
theory. 

In a system with unit acoustic velocity it is assumed that the axial 
velocity of the body is O(1) and that the cross-flow velocity is less than 0-4. 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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If these conditions are satisfied, then the fluid perturbation velocity wil] 
be everywhere small compared with the acoustic velocity, and a linear 
method of approach will be sufficiently accurate to first order in €, provided 
that the body accelerates rapidly through the transonic region. The fluid 
perturbation velocity is the gradient of a scalar potential function, 4, 
which satisfies the wave equation on a linearized theory. The linearized 
potential, dy, at points distant r from the axis of the body is obtained asa 
descending series in r, each term of which represents a source or multi- 
source potential satisfying Laplace’s equation in two dimensions. Thus, 
the flow is essentially two-dimensional in any plane transverse to the axis 
of the body. 

A set of cylindrical polar coordinates (r, 0,7), referred to axes moving 
with the body is taken, with 7 measured along the axis of the body from 
its base, and r and @ are polar coordinates in the plane » = constant. Let 
(i,j, k) be a triad of unit orthogonal vectors, with i along the axis of the 
body, and j in the vertical plane through 6 = 0. The corresponding fixed, 
cylindrical polar coordinates (r,,6,,8,), have their polar axis along the 
direction of mean motion, with s, = 0 at the initial position of the base 
of the body. The displacement of the base of the body along the fixed 
polar axis after a time ¢ from rest is s, f(t). The initial condition that 
the body starts from rest at t = 0 gives 


f(0) =f'(0) = 0. (1) 


The corresponding fixed triad of unit orthogonal vectors, (i,.j,, K,), has 
i, along the fixed polar axis and j, in the plane 6, = 0. 

If, in addition to the above assumptions, the fluid is also non-conducting, 
then variables p, p, and u, representing the pressure, density, and perturba- 
tion velocity respectively, satisfy the equations 


oeode V(pu) = 0, (2) 
ct 


Cd * dp — 


and -}u? = a spatial constant, (3) 


ct p 


where u = Vd. (4) 
Also for a perfect gas whose motion is isentropic, 

p = Kp’, (5) 
where K and y are known constants. In a system with unit acoustic velo- 
city, these equations reduce to the linearized equation 


ed ; 


V7¢ = —. (6) 
ot" 
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It is well known that a solution to (6) is the source potential 


q(81, t’) 
~_ 


This source, at a point (0,0,s),¢’) in a four-dimensional space-time con- 


where t’ = t—R,, and R? = (s,—s})?+73. (7) 
1 


tinuum, gives rise to a disturbance at the point (r,, 4,, 8,,¢) wheret = t’+ R,. 
For motion at a variable incidence of small magnitude, the sources are 
placed on the axis of the body, and, since s,—f(t) = y(1+O(e?)), the 


source potential is 


q(s’, t’) . 
,1+ O(e*)}, 
R 
where t t— R, to (s—s’)*+-7?, Ss n+fi(t). (8) 
_ &" (g(s',t)\ ; . (6 der; fi 
Clearly - is also a solution of (6) to the same order; for n = 1 
ox” t 


this is a doublet potential, whilst in general it represents an nth-order 
multisource potential. The general solution to (6) is obtained from a 
combination of source and multisource potentials of all orders. The 
multisource densities are obtained from the boundary condition stipulating 
that the normal velocity of the fluid and of the body must be equal on &, 
the surface of the body. 

\lthough a source method has been used in this paper, the series for the 
potential at points near the axis of the body may be obtained by an opera- 


tional method (see Ward (3)). 


Part | 
SLENDER BODY OF REVOLUTION 


A. Direct motion 
1. Solution for the velocity potential 
This section is concerned with the problem of the motion of a slender 
body of revolution along its own geometrical axis at a variable forward 
speed, f(t). The motion is represented by a line of sources placed along the 
axis of the body; these sources give the correct fluid displacements provided 
that the boundary condition on > is satisfied. The total potential is a sum 
of elementary source potentials, viz. 
$,(r, 0,8, t) ie O) as’, (9) 
F t 
G 
where the region of integration, G, includes all points on the axis of the 


body at which there is a source at a time t’ = t— R (an implicit equation). 
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Since the aerodynamic force is obtained from the normal pressure inte. 
gral over &, only values of ¢, at points on and near to = will be required 
in this work. 

Thus the only required values of ¢,(r, #,8,¢) are those for which 
(1) » = s—f(t) lies in the range, (0,/), where / is the length of the body, 

(10) 

(2) the distance of the point (7,@,s) from the axis of the body is small and O(e), 
(11) 

The region of integration includes all points, (0,0, s’,?’), for which 
(1) 7” = s’—f(t’) lies in the range (0,1), and | 1. 

(2) tt’ = R= +fr2+(s—s')3. | "4 
In the (s’, ¢’)-plane, ‘—?t/ = R > 0 represents one branch of a rectangular 
hyperbola with centre at the point (s,¢), and asymptotes (t—t’) +-(s—s’), 

Contributions to the potential, 4,(r, 6, s,¢), come from those regions of the 
body’s trajectory in the (s’, ¢’)-plane which are intercepted by the hyperbola 
t—t’ = R (see Figs. 1 and 2). For subsonic velocities only one region of 
integration occurs, but for supersonic motions there are usually two regions 
of integration, except when the point P,(7r, @,s,t), at which the potential is 
required, lies near to the body both in space and time; in this case the two 
regions coalesce (see Fig. 2). Except for transonic velocities, the region 
of integration is O(/). 

If v is the outward normal to © and r = 7(y) is the equation to the body 
surface, then from (9), the boundary condition gives 


G 
On assuming that q,(s’,/’) satisfies the Lipschitz condition, 
q,(8, t)—q,(s;, t’)| < A’ |s—s’|4+ B’ \t—t’|, (14 
where A’ and B’ are constants of O(e?), we have 
staan q,(s. t) (s—s’ e . 
Log - 1 5 | + Ole), (15) 
since A’ ris s'|F ds’ - O(c?) at most; (16) 
es 


and a similar result holds for the term containing B’. In this work @ con- 
sists of one region, only, bounded by the points, s’ = s, and s’ = 8. In 


general |s,—s| and |s—s,| are both O(/) and, from (15), 


f') 


di 


2q,(s, t){1+O(e)}. (17) 
dy 
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For some points, however, |s,—s O(e), and/or |s,—8s O(e). In these 


ses 


2-+- O(e*) 





the comparative rarity of these points allows of their neglect. If S(») is 
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Fic. 1. Subsonic motion. 
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Fic. 2. Supersonic motion. 
Thick lines in Figs. 1 and 2 indicate regions of integration. 


the area of the cross-section of the body, then the potential appropriate 
to this motion is 
. Fe id ” “7 t’ * 
| ey, (19) 
4nk 


p,(", 6,8, t) 


G 
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where n” = 8’—fi(t’). (20 
In the Appendix, ¢, is reduced to the following form, 


F()S"(m) 


mT 


$,(r, 6,8,t) = log $r + 


log(t—t’) Af (YS (n/ +t—t’)+-8’(y’+t’—1t))\ dt’, (2) 
7 Cc 
H 
where 7’ = s—f(t’), and the region of integration, H, includes all values of 


t’ for which the arguments of S’ lie in the range (0, /). 














A 
ES. F t CS, D 
/ 
F 
~TIrayectory of nose 
L - 








/ s’ 


Fic. 3. Influence regions of the body AB. 


If A,, Ay = s—f(t’)+(t—?’), then the condition that A, should lie in the 
range (0,/) is equivalent to the condition that S, (see Fig. 3) should lie in 
the domain of influence, CD, of the body AB by waves propagating in th 
direction s-increasing. A similar result holds for A,. 

2. The drag force 
The total inviscid aerodynamic force acting on the body is F, where 

F — (p—ppo)v dd, (22 


and p, is the undisturbed pressure. 
From Bernoulli’s equation the pressure is given by 


Pp e hy. 1 /(&,\* | O(e4 loge), 
Po a 2 


to 


cr 


since p po( 1+ O(e? log e)). (24) 
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The time derivative @/ét is taken with respect to axes fixed in space. The 
irag force, D,, obtained from (22) and (23), is 


l 
> 8 0 19 "9 r(O) " r "9 7(m) 
), = —F.i | f'2(t)S *(O)log| |—2f"(t) | S *(n)logh | an 
0 
: : 
| Ss’ ”n) dy log(t—t') ar| at O(e® loge), (25) 
0 H 
where 
A = f"(t’){S'(b,)+8' (bo) f(t) +f (UY S(O US" (h) +8" (he)f. (26) 
and bb, = s—f(t') +(t—t’). 


This expression (25) holds for both sub- and super-sonic motions when the 

body is pointed at each end. For a body with non-zero base area, however, 

25) holds only when the motion is supersonic (and a base drag force, 
, S(O) must be added in this case). 

If if” (t) f’*(t), and the motion is supersonic, the drag force has the 


following approximate form: 


ll 
D, Po f’= +)! | | log Y oC S"()S"(a) dndo— 
dcr ‘ | 
; Br(O 
28'(0) | log(d—n)S"(n) dy S0)log(—2”)|- 
| F Br 6f ’2(t 
f”(t)\2 | S’(q)(log FOF *( \ dy 4 
\ . 2 B 
: . 
2 | 8") dq | S"(a)log(o n) do||-+-O(elog’e), (27) 
where R2 ft) 7 


In practice the effect of acceleration on the drag force will be small 
compared with the steady velocity effect if the change in the velocity of 
the body in the time it takes to travel its own length is small compared 
with the actual velocity of the body. 


If the drag force (27) be written in the form 


D, = Ay f(t) + By f’(b), (28) 


U0 
where Ay, B, are given by equation (27), then the effect of the acceleration 
is to make the virtual mass equal to m+ By, where m is the mass of the body 
ind By = O(p,e2V), where V is the volume of the body. 
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B. Motions with oscillations 

1. Introduction 

In this section the motion of the body is taken to be a general oscillation 
of small amplitude superposed on a variable forward speed. Cylindrical 
polar coordinates, (7, @, 7), again refer to axes moving with the body, and 
D/ Dt is a time derivative taken with respect to these moving axes. The 
angle between the geometrical axis of the body and the line of mean motion, 
and its time derivatives are O(e). Initially the body is considered to possess 
an angular velocity, «’(¢)k, about an instantaneous axis, 6 = $7, » = )(t); 
but later an angular velocity, ’(t)j, about an axis, 6 0, 7 no(t), is 
superposed on the above motion to give a general spiral motion. In the 


former case the velocity of any point of the body is 


U = f’(t)i+9'(t)j, (29) 
where gq’ (t) x’ (t)(7 — no(t))— a(t) f'(é). (30) 
The boundary condition on = gives 
$,+¢4> re ae 
. (Cr+ $s)) ~ f(t) —g'(t)cos 8, (31) 
cr !. 3 dn 


where ¢, is the potential due to the oscillatory motion. This boundary 
condition is satisfied by taking ¢, to be the potential of a line distribution 
of doublets oriented in the j-direction; thus 


* d lacle’.t’)\ 
dbs : 1s’. 32) 
2 | P| R : \ 


G 
For values of r of order € this doublet potential may be reduced to the 
following form, 
cos @ ; 
do 2qo(8, t) 1 Oe), (33) 
: 
The value of q,(s,¢) is now obtained from the boundary condition (31); 
thus 
: f cos @ ‘ 
$2. = —g'(t)S(m) (34) 
Tr 


2. The aerodynamic forces 


[f R, is a position vector relative to the instantaneous axis of rotation, 
then Bernoulli’s equation gives 


) ) D © ~ 
ad - 4(V¢)? , U,.Vd 4 (R, \ Vo). x’ (t)k- O(e4 loge) (35) 
Po 


(V refers to axes fixed in space and U, 


f' (ji—a(t)f’ (Oi). 
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The drag force, D,, is obtained from a reduced form of the pressure 


distribution, (35), thus D, F .i, is given by 
D—D, ' ‘a ee rd , ae . 
£S(O){ax?(t) f’2(t)—a’2(t)n§(t)} + Vialt)a’(t)(2 f'()+ nolt))4 
Po 
(no —7,)(a'2(t) + a(t)ax”(t))+07(t) f"(H)}+O(e8 loge), (36) 
: 
where Ny V | 7S(n) dy. (37) 
0 
The effect of linear and angular acceleration on the drag force will be small 


mpared with the steady velocity effect provided that 


if’(t) < f(t), and [a"(t) < a(t)f’*(t) (38) 


respectively; also the effect of angular velocity will be small provided that 


'(t) < a(t) f(t). It is interesting to observe that larger forward velocities 
low of larger angular velocities before the latter become effectual in the 
lrag force. 
The lift force in the plane 6 = 0 is denoted by L,, where, from (22), 
34), and (35), we have 
L.=F.j, pol Via’ (t)( f(t) + no(t)) +=" (t)(nolt)— 4) + oO f"(O}+ 
Lf (t)S(O)fa(t) f’(t)+a’(t)yo(t)} | +O(e* loge). (39) 

Let '¥'(t) be the angle between the resultant velocity of the nose of the body 
ind the axis of the body, viz. 
x'(t), 


t(0) u(t) —(l— ny) ——; 
a (1) 


(40) 


n terms of ‘’(¢) the lift foree becomes 
L, = pA iVOP (OTP (OSL (O+(— aq, )o" (Oj 
+ f(t)S(O)fla’ (+O S'(O}). (41) 
If conditions (38) are satisfied, then the effect of both linear and angular 
icceleration on the lift force will be small compared with the steady velocity 
effect. 
Li, 7 = 0, then 


> 


Let M° be the pitching moment about the axis, @ 
MY = p, Vin’ Of (O+VOSL’ (O} + VOL ()- 


v" (t(j, — A) +(l— Fy )a' (Of), (42) 
l 
where n>} | 7*S(n) dy. (43) 


ig 


The analogous pitching moment about the nose is 
M! —I1L,—M?. (44) 


For a general oscillation with angular velocity a’(t)k+-f’(t)j, the above 
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method is still applicable; in this case the potential at points near to the 


body is S(n) 
¢,= +> I’ th’ (t)sin 8—g’ (t)cos 8}, (45 

27r 
where h'(t) = Bt) f'()—P'(t)(n— n2(?)). (46 


The additional lateral force, introduced by this rotation about the 
x-axis, lies in the k-direction and is of magnitude 
L, = F.k = —p V{B'(t)no(t) +B" (t)(nalt)— 4) +B’ OS (O+BOF”(O}+ 
f(t) S(O){B'(t)n(t)+ BOS (Oy). (47 
For a body moving at constant incidence « (8) to its own geometric axis 
the aerodynamic forces are given by the above equations with 


x’ (t) 0 (B’(t) 0). 


Part I] 
SLENDER BODY OF GENERAL CROSS-SECTION 


1. Introduction 

In this second part the preceding work is extended to solve identical 
problems for the slender body of general cross-section. 

The method is still to choose a solution of the wave equation which 
satisfies the boundary condition on ©: whilst source and doublet potentials 
were sufficient for the previous motion, in this problem higher-order 
multisource potentials will be required. 

If (r,@,) are a set of cylindrical polar coordinates whose polar axis 
coincides with the axis of the body, then the potential of a distribution of 
nth-order multisources oriented in the j-direction is 


r (d\*(q,(e',t')\ ;., pie Gree 
| 4. ds’, where R? s—s’)*-+-7?, (48 

. () | R J ital ’ 

G 
For r O(e), this may be reduced to the following form (see (33)), 
> ‘ ’ 
d, A. —_ Les (49 
pa-i 

where A,5 (—)"-!2(n—2)!q,,(s, t). (50 


The general expression for dy, the potential at points near to the body, 


may be reduced to give 


= A. ecos(n@+ ; 
do a, log r-— b)+ * n ( B,) (51 
y a) rl” 
n 1 


where §,, is the phase angle giving the orientation of each multisource 
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distribution; 8, may depend on 7. A complex potential w, whose real 


n 


part is 5, may be defined by 


x 
. a " 
Ww Po T ibs a, logz-+ by+ > — (52) 
n * 
where a, = A, e'?”. Thus, dy satisfies Laplace’s equation in two dimen- 


sions, so that the flow is essentially two-dimensional in any plane of cross- 
section 7 = constant, and the multisource densities a, are determined 
from the solution of a two-dimensional incompressible flow problem in 
which the velocity at infinity is zero. 

The body is taken to have an angular velocity, «’(t)k+-f’()j, with in- 
stantaneous axes as specified in Part I. 


The velocity of any point of the body is 
U = f'(t)i+g'(Oj—h' (Ok, (53) 


where g’(t) and h’(t) are given by (30) and (46) respectively. From the 


houndary condition 


U. (eo (54) 
l\avJs 
it is found that the source strength a, is still given by 
Ql, S'(n) f(t). (55) 


It is convenient to determine the aerodynamic force from a momentum 
integral rather than from the normal pressure integral directly (see (22)). 
For this we consider a right circular cylindrical control surface, S, sur- 
rounding the body and fixed relative to it. The control surface consists of 
. curved surface S,, of radius r, = O(e), and circular ends S,, S3, lying in 
the planes 7» = 1, » = 0, respectively. 

If D/ Dt is a time derivative with respect to axes moving with the body, 
ind J; is the volume of fluid lying between S and &, then the normal pressure 
integral (22) may be transformed to the following equation, 


: ° D(pu) 


F {(p—p,)/+pu(u—U)}.ndS— Di lv, (56) 


t 


« « 


S Wi 
where nis a unit vector normal to S, and J is the idemfactor or unit dyadic, 


The pressure, obtained from Bernoulli’s equation, is given by 


} } db n ( , 9 -=- 
(" f °| ms Lv)? p(y oP g (1)°? + O(etlog?e). (57) 
Po Dt 2 7 cn Cx 
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2. The drag force 

The drag force D is minus the component of F along the direction of meay 
motion, i,; from (56) and (57) we have 

) , , , , , 
D ("(1)oo_# Po! (vs)?! dS + D obo dv + 
Po J ox Dt 2 J Dt\ én 
Si—Ss V; 

[ (a 


cr 


: -g' (t)cos i) dS,+-aL, + O(¢ loge), (58 
on 


. 


where L,, is the lift force defined later by (78). 
The several terms of (58) will be reduced separately. Let C, be the con- 
tour of the body in the plane » = 0, then from Gauss’s theorem and the 


equation for the potential, (52), we have 


1° _.. If, ad | 
2 | (V$o)° dS = 5) > $o ae dr {Ay bo} n ee G(r), (59 
Sr S3 Co 
where G(r,) is a function of r, and the a,,’s (see (61)). The divergence theorem 
gives , ; 
—— | °dS— | —%(v.i) dz. (60) 
re Dt ‘: Dt Dt ( ) 
Vi Si-Ss3 z 


The integral over S, may be determined explicitly from the equation for 

¢d,; thus 
Odbo{Obo _ 
cn | cr 


q'(t)cos a dS, 


l 


> 


al 
2n dy" dy—G(rs)—na,9'(O-| 7 [ «, x(t) dn. (61 


ho 


: "lie 
‘ ] ‘ 
0 0 
Let W be the area in the plane 7 = constant bounded by the circle, r = ',, 


and the body contour, 7 = 7(7, 6), then let 


I(m) | [bo La obo| dw, (62 
J lea oy | 
Ww 
and so | Cho dS = refI(l)—1(0)}. (63 
ox 


On defining C to be the contour of the body in the plane 7 = constant, and 
C, to be the circle 7 = r, in the same plane, we may apply the two-dimen- 
sional form of Stokes’s theorem to give 


I(n) = —i{ $ dz— $ a] fj v(—1)}; (64 
C1 C 
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fom the complex potential defined by (52), we have 


I(n) 7,-+-0 { w dz $ ub, dz. (65) 


( c 
Since w is analytic in the z-plane cut from the origin to oe’ and all the 
singularities of w are enclosed by C, then the series for w converges in the 
neighbourhood of some point z z, on C. If C. is another closed contour 
having 29 in common with C, such that the series for w converges every- 


where on C’., then the second term in (65) is equal to 


a } w dz 279 2 27;. (66) 


CU, 


On integrating by parts, we have 


4 oy 3 
® wb, dz [zbpde t 210 dr 22Be Ze | 2 obo dr, (67) 
, ; 4 OT b Ov 
( Cc 
a: 
ind so I(7) 70, — > z 0 dr. (68) 
J & 
s 
Ifz = z,() is the centre of area of the cross-section S(7), then 
a : ee ee 
» f'(t)(v.i) dr f(t) — {S(n)z,(n)}- (69) 
. ‘ dy q 


Also, from the divergence theorem, 
b zv.(g'(t)j+h'(t)k) dr ig’ (t)+th’(t)}S(m). (70) 


On simplifying (68) by means of (54), (69), and (70), we have 


— 2 P oe ‘ 
I(n) 7,+f (t) ; M()zZ,(n)s—(g (A th’(t))S(m). (71) 


If A(t) = 0, the drag force, obtained from (58) et seq., is given by 


D Lf Ch * Dd, (ed d> 
D bo—— dr roy 10 __ g'(t)(j.v)} ——+ 
| > a iain oo oe 
1 
\ Cy . da P : : 
7 | (a, v (7) 2b 1 2a(t) De dy- re (t)S() 27a,— 


f'(t)- (S(m)24(0)}| x’ (t)no(t))—7a9 bo OFS) -- 


c 
WVQ+ O(e* loge), (72) 


where 
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The drag force for an oscillating slender body of revolution is obtained for 

ma, = —g'(t)S(n) (74 
(cf. (36)). For direct motion without oscillation, a(é) 0, and the drag 
force is D, where 


l 
ca 
05 
7 


Dy i $ 3 Ody a [ Ddo dbo dX _ oF °dn — 
Do $7 °° & | Dt ev f'(t) én 
Co p>» 0 
— (A bo)y-9 + Ole8 log*e). (75) 


The third term of (75) has the alternative form: 
l t 


"ay oP 4 
m &. S'(n) dn | log(t—t’) : 
aT C 


0 0 


9 
7 log 3(45)» 0 


[f(t ){S"(n' + (t-#’))4 


+ 8" (m’—(t—t’))}] dt’. (76) 


Contributions to (76) come from regions of ¢’ quite near to ¢; thus if t’ lies 


in the range (t—&, t), then for f(t)— f(t’) = V(t’), we shall have = = J 
oo 
when V, > 1; whilst, for ¥, < 1, 
bit [= 
2. > a . “oF k aS (77) 
yH+1 I—l 


3. The lateral force 
The lift force in the j, direction is L,, where, from (56) and (57), 


L, = F.i, 


Z 


po f ‘(t) Cho dS | Po | (2 Pecos A a $$ (20 _ aitvoos a\\ dS . 
) ee ee aa \ er 1 
* D [ed — . 
Po il set) de + Ole log’ (73) 
Vy 
A complex lateral force, L, may be defined by L = L,+7L,; in the above 
notation, 
l 
2 D) , 
L ! ¢°(t) — —\(1(n) —za,) dn + O(€* log?e). 79) 
—As oe | oF 
0 


On substituting the value of /(7) from (71) into this equation, we have 
£ 


Po 


f”(t)S(0)z,(0)—f ?(t){8’(0)z,(0) + S(0)z,(0)}— 


, 
g 


| O(e® loge). (80) 
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The angular velocity, 6’(t)j, introduces extra terms, 7L,, into (80), where 
L, is given by (47). 
A complex pitching moment, ./°, is defined as the moment of the lateral 


force about the base; (71) and (79) give 


l 
) 0 ° 0 ° P 7 0 D 3 
M' M°®+iM? ipo | n ans (0< — Gln) Fay) (81) 
L : 
If n?S(n) dn nV, and Z4(n)S(m) dn = Z,V, (82) 
0 0 
then 
l 
; , hs Da,)\ z ; An Soe 
ye po] 277 \J (t)a,;+7 Di | dyn+f (t)S(0)z,(0)-+ f (t)z, J .s 


« 
0 


+ V| v(t) f ‘(t)no(t) ; u(t) f*(t)— x" (t) A + 
+ Hy{x" (t)jo(t)+-a(t)f”(t)+-a'(t)no(t)}]]. (83) 


Thus although the drag force depends on all the multisource densities, 
1,, the lift foree depends only on the first one, ay. 


APPENDIX 
Reduction of the Source Potential 


In this appendix the source potential (19) is reduced to a more convenient expres- 
m, (21). 


The source potential is 


bir, 0 Ll per) (a.1 
b.(r,0,8,t (t’)S’(7 re a. 1) 
1 ro ‘  _ 
here t t—R, 7” s’—f(t’) and R (s—s’)?-+-r*. (a. 2) 
r points (7, 0,8, t) near to the body both in space and time, @ consists of one region 
nly bounded by the limits s’ s, and s’ S20 
On integration, we have 
“ds (R 8 8’) ” 7 
R [oe r I: _ 
Pines mav be integrated by parts to give 
(R S s’\ 1°? r (R s 8’) c aa W\ gry 4? , 
dard, f’(t’)S’(7 log) 5 i]. — | log) 5) jz )f (t’)} ds’. 
83 
(a. 4) 
: [Rie—s [R+s’—s)_ |) (4 ; 
Bi i “ log, = —log| 5 | —log| =). (a. 5) 
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and so 


89 





R+s—s’ 
33 
(pF vs. er : ie dhdiciadil 
: | log{- Pas {S8'(n") f"(U)} | ds’ +2 log(4r)[S’(”) f(t’). (a. 6) 
. . \ - les : 7 


If the first term in the right-hand side of (a. 6) is denoted by Q, and 


{R+ |s—s’|) 








E = f"(t’)S"(y")log) 1, (a.7 
then (a. 4) and (a. 6) give 
4d, Q+[E]y =e; +[E]——2.—2S'() f(t)log $r. (a.8 
At both limits, s’ = 83, s’ = 8,, we have 7” = either 0 or J, (see Figs. 1 and 2), 


Consider the expression E: 

(a) for bodies pointed at each end, FE vanishes at both limits since S’(0) = 0 
and S’(1l) 0: 

(6) in the supersonic case one region of integration, only, occurs when the point 
(r,@,s8,t) lies on or near to the body; in this case both limits correspond to 
n” l and, since S’(/) 0, EF vanishes at both limits even when S(0) + 0, 


It can be shown that 


, 2 > (R fs Nicer 9’ | C = sapere . : 
log; ——_——' ] — {S8’(n") f (t’)} | ds’ = O(e*), (a. 9) 
\ z les . 
83 8 P 
r ¢ c vs ” , , « 
and so Q = ( : }[to R—{8"(") f(t | 1 O(e*), (a. 10) 
cs 
For |s—s’ O(l), we have (t—t’) R |s—s’|+ O(e*), and 
n” s—f(t’)+s’—s s—f(t’)+(t—t’)+ O(e?). (a. 11) 
But values of s’ for which |s—s’ O(e) only contribute an amount of O(e* loge) 
to Q@ and so 
‘ rn, © lial ” , , 7 , , , > 
Q log(t—t’) PAL EMS (+ t— 0) + Sy —t +E) dt’+O(e*? loge), (a. 12) 
C 


H 
where H extends over all values of ¢’ for which the arguments of S’ lie in the range (0,1). 
From (a. 8) and (a. 12) we have 
ra , , c , , vs , , 
4rd, 28’(n) f (t)log ir log(t—t’) pid (t’)(S"(n’ +t—t’)-4 
c 

H 
+ S"(n’—t+0’))} dt’ O(e loge). (a. 13) 
The author wishes to express his sincere thanks to Mr. G. N. Ward for 
his help and encouragement and to Professor Goldstein and Professor 


Lighthill for their advice and criticism during the preparation of this paper. 
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THE ENERGY DISTRIBUTION BEHIND A DECAYING 
TWO-DIMENSIONAL SHOCK 
By J. E. PHYTHIAN 





(Department of Mathematics, The University, Manchester) 
[Received 10 July 1951] 


SUMMARY 

In this paper the work of Lighthill (1) for one-dimensional unsteady motion is 
extended to the problem of steady, supersonic flow past a two-dimensional slender 
body. Goldsworthy (2) and others have shown that, in the limit when the Mace! 
number, M, of the undisturbed flow tends to infinity, the two motions become 
identical. 

As in (1), an accuracy hypothesis is formulated concerning the degree of accuracy 
of the Friedrichs theory (3) for nearly plane, two-dimensional, steady shocks, “A 
comprehensive check of this theory on the lines of the accuracy hypothesis is mad 
for a fluid with general thermodynamic properties: in the process of justifying this 
hypothesis the picture of the flow as a whole is improved. 

A relationship is found between the increased total energy flux across a plan 
transverse to the flow in the residual wave, and the entropy flux across the region of 
the shock upstream from this plane. The relationship is checked by comparing the 
value of the pressure fluctuations in the residual wave (due to entropy change at the 
shock) as determined from a consideration of conditions in the residual wave, wit! 
the value obtained from an investigation of the local reflection of the simple wave it 
the shock. 

The work is extended to the problem of the actual flow past a two-dimensional aero- 
foil when there are two shocks forming an N-wave, and in this case, conditions very 


far downstream are investigated. 


1. Introduction 

Tue work of Lighthill (1) for one-dimensional piston motions is extended 
to the problem of the two-dimensional, steady, supersonic flow past a thin 
symmetrical body. Goldsworthy (2) has shown that these two motions 
become identical as the Mach number, //, of the undisturbed (two-dimen- 
sional) flow tends to infinity. 

For both these motions, Friedrichs (3) has put forward theories which 
neglect the third-order entropy change at the shocks. By assuming an 
accuracy hypothesis for the Friedrichs theory for one-dimensional motion, 
Lighthill has been able to sketch into this earlier theory the effect of the 
third-order entropy change, and then to check the hypothesis. 

In this paper a similar accuracy hypothesis is formulated for the Fried- 
richs theory for steady, two-dimensional flow, and results similar to those 
derived by Lighthill are obtained. At first the decay of a single nose-shock 
attached to the body is investigated, but later this work is extended to 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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the problem of the NV-wave. For this two shocks occur, one attached to 
the nose and one attached to the tail of the body: the region lying between 
these two shocks is a simple wave. 

On the Friedrichs theory the flow behind a single shock is a simple wave 
followed by an undisturbed region—the residual wave. Errors in this 
theory are chiefly due to: (1) neglecting the third-order entropy change 
it the shock, (2) applying a second-order boundary condition at the shock. 


These errors suggest the following: 


Accuracy hypothesis 

If the simple wave is altered in the following manner, then the changes 
produced are of the same order as the errors arising in the Friedrichs theory: 
|) the specific entropy of the whole wave is altered to its known maximum, 
2) the slope of the body surface is altered continuously by an amount of 

order the cube of the maximum shock strength. 

If this hypothesis is valid the Friedrichs theory will be of considerable 
value at any distance downstream. The theory is checked by comparing 
the value of the pressure fluctuations in the residual wave, as determined 
from a consideration of conditions in this region, with the value obtained 
from an investigation of the local reflection of the simple wave in the shock. 
The work is modified for the .V-wave, and in this case conditions well down- 
stream of the body are investigated. 

A brief résumé only of the work is given, since the analysis is similar to 
the work of Lighthill (1). 


2. The entropy-energy flux relation 


The notation for the physical quantities used in this paper is as follows: 


Quantity Und sturbed flow Disturbed flow 
Pressur¢ r P+p 
Specific volume J V—v 
Specihe entropy 2 0 8s 
Fluid velocity U, 0) (U+u, w) 

w : 
Flow angle A 0 6 = 0(68) 
l u 
. (ep 
Sonic velocity . 7 VaPy a (V—v) | ) 
N \cv/s 
Temperature , ‘ r 


In the undisturbed flow, derivatives will be indicated with capital letters 
8 suffixes; thus P,. is the derivative of the pressure with respect to minus 


he volume at constant entropy. The usual partial derivatives will be 
sed for the disturbed flow, i.e. the afore-mentioned derivative is (ép/év),. 
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The steady, two-dimensional, isentropic equations of fluid motion may 
be solved to give values for the fluid perturbation velocity components in 
the region of the simple wave; these are 


w = UfO—6?(M?2—1)-!+ 0(68)! | 


4, P 
u = U{—6(M?—1)-!+ 6>(M?—1) nly MVP, y | 


. 1 O(03) {° 
| 4(M?—1)P,-} \ | 


The flow behind the shock is a simple wave in which the flow quantities are 
constant on straight characteristics inclined at an angle (tan~'Z_,) with the 
axis of the body. If €(@) is the intercept of a ¢, characteristic on the y-axis 
(see Figs. 1 and 2), then these characteristics have equation y &(0)+al,, 
where ¢., is obtained to second order in @ from simple wave theory. 
























ch | yA 
> J 
| Shock 
a Simple 
| ; wave u - + 
| ris ie 
| i) ‘ 
52> Residual - aco 
/ wave | 
aa, we ; | | wave 
—_— Piston curve body surface 
t Eat” r 
Fic. 1. One-dimensional un- Fic. 2. Two-dimensional 
steady motion, steady flow. 


The specific entropy gain at the shock is determined from the shock con- 
ditions. Also, the rate at which the body is doing work in moving through 
the fluid, and the flux of increased total energy across that part of the plane 
x constant which lies in the simple wave (see Fig. 2) are both determined, 
to second order in @, from the simple wave theory. The difference between 
these two expressions (a second-order quantity) must represent increased 
total energy flux across the residual wave; the latter is denoted by R(z). 
The form of R(x) suggests that it depends only on shock discontinuities 
upstream from the plane x = constant; if S(a) is the total entropy flux 
across this region of the shock, then it is found that 
T,-P PP,y M? \. 


. (2) 
2TP, 8P%(M?—1)} 


R(x) TS(x)) _ 


; : U(x 3M? } “i 
for a perfect gas ad — (3.M"—4)y iJ (3) 


T S(x) Sy( M?—1) 


where y is the ratio of the specific heats. 
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Equation (2) is identical with the result found by Lighthill (correspond- 
ing to the case M 00) except for the factor M?/(M?—1). 


3, The pressure pulse check 

On the accuracy hypothesis the velocities and thermodynamic variables 
ye of third order in the residual wave. The pressure fluctuations in this 
region satisfy the wave equation to a first approximation; and so a solution 
for p under the boundary condition, [dp/dy 





. 0, is 
p Fi{x+y,(M*?—1)}+ F{x—y,/(M*2—1)}. (4) 

Moreover, the explicit value of p may be obtained from relationships 
existing between the third-order velocities and thermodynamic variables 
in the residual wave; thus 

VTM (47',, P,ypM* | 8 

a ot XK 

16a9,(M?—1)\ T P,(M?—1) © M?} 





id . 2 ‘td 2 ' 5 
S’{4(a+-y,/(M?—1))}+ S’{4(x—y,/(M?—1))}]. (5) 
The pressure variations in the residual wave are due to a pressure pulse 
propagating from the shock along ¢-characteristics and of strength 


ét y , dS(2) 
{(3y—5)M4— 4(y—3)M2—8!} 
16U(M? pi? 1 \Y ) dx 


; (6) 
fora perfect gas. The first term in the square brackets in (5) represents a 
pressure pulse which may be considered as a reflection of the simple wave 
n the shock; this pulse is propagated along characteristics with slope 
lydx = €_ (< 0). The second term in (5) represents the reflection of this 
reflected wave, at the body surface, and is propagated along ¢, character- 
istics. 

Consider the reflection of the simple wave at the shock. The reflected 
pressure wave, determined from a study of local conditions just behind the 
shock, is propagated unchanged (to a first approximation) through the 
simple wave and into the residual wave. The theory is checked by com- 
paring this reflected pressure wave with its value given by (5). 

The pressure gradients along the two families of characteristics are 
btained from a knowledge of the shock discontinuities and the exact 


haracteristic equations; these are 


dé (W2—a?)t(V —v) dy , = 
t , on — = ) 
dp aw? oie dx 3 (7) 
where W2 — (U+-u)*?+w? (see 4). (8) 


5092.19 


Y 
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A reflection coefficient, 7, is defined as the ratio of the pressure gradients 
across the simple wave and its reflection at the shock ; thus we find 
14 p2 ” 2177p 9 
M*P+,, {TJ M?V Pypy Pe. \,3 


128V P3(M2—1)2\ 7 P,GP—1)* We" 


where v is the specific volume perturbation. 


(9 


From this value of r and the value of |dp/dx| along ¢, characteristics, 
obtained from simple wave theory, the equation for the pressure pulse 
along ¢_ characteristics is obtained: it is, in fact (6). 

The pressure pulses reflected at the body surface gradually slow down the 
characteristics in the residual wave. 


4. The N-wave 

The previous work for a single shock may be extended to the case when 
the body is finite in length, so that two shocks occur, one at the nose and 
one at the tail of the body. Both shocks decay in strength as they recede 
from the body whilst the flow between them is taken to be a simple wave. 

On the Friedrichs theory, the fluid behind the rear shock is undisturbed 
to second order, and the rear shock conditions are satisfied to this order, 
Since conditions behind the rear shock are known with an error of order 
the cube of the maximum shock strength, this error may also be incurred 
in the flow quantities in the simple wave. Thus, at large distances down- 
stream, the motion of the rear shock may be given quite inaccurately by 
simple wave theory. 

The work is developed in a manner similar to that outlined above for a 
single shock. The total increased energy flux across that part of the plane 
x = constant which lies in the residual wave, and the entropy flux across the 
regions of both shocks lying upstream from this plane, are both determined 
explicitly, and relation (2) is found to hold between them. This work is 
checked by a pressure pulse method. The pressure pulse is still given by 
(5), although, because of the new value of the entropy flux, the pressure 
pulse becomes 

M®P,,V? (47,V M*VPpy , 8 \ 
192(M2—1)52| T =P, M?—1)' My 


where 6,, #, are the values of @ at the front and rear shocks respectively. 


(63 — 63), (10) 


In the latter equation the term in 6, represents the reflection of the simple 
wave at the shock, whilst the transmission of this reflected wave through 
the rear shock gives rise to the term in @,. 


At large distances downstream, the equation to the rear shock is un- 


reliable, for the reflected pressure waves will gradually overtake the rear 


shock and the equation given by the Friedrichs theory will be modified. 
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Verv far downstream the energv flux is divided between the following 
odes: 
a) energy transferred to infinity in the N-wave, 
energy flux across the residual wave associated with entropy flux 
across the shocks, 
energy flux across the residual wave due to the pressure wave behind 
the rear shock. 
Far downstream the ratios of the energy fluxes in these three modes 
ecome 
T, P M?PP,, \-. (1 T, P rt 
Eras SP?.(M2—1 ’ TP, wvR) 
PT, PPyy M? L P (11) 
2P,T s8P}(M*—1)° M?*VP, 
The total energy flux is 7’S(co) in this case. 
This is not the final state, however, for ultimately the energy associated 
ith the pressure wave, (c), is gradually transferred to the N-wave, (a), 
nd all the energy which has not gone into increasing the entropy of the 
iid without pressure change is in progress to infinity in the N-wave. 
[his final energy transfer occurs at a distance downstream of O(2;/6,), 
shere a, is the length of the body and @, is the value of 6 at the nose. When 
the body is not symmetrical, the energy flux in the N-wave dominates the 
nergy fluxes in the other modes. 
The author wishes to thank Mr. G. N. Ward for his help and encourage- 
ent and Professor Lighthill, who suggested the problem, for his criticism 


nd advice. 
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THE DIFFRACTION OF SOUND PULSES BY AN 
OSCILLATING, INFINITELY LONG STRIP 


By F. J. BERRY 


(Department of Mathematics, The University, Manchester) 


[Received 26 June 1951] 


SUMMARY 


In this paper the solution is obtained for the problem of a plane pressure puls 
normally incident on an infinitely long strip of finite width, capable of motion a 
a spring-supported rigid body. 

Numerical calculations have been performed for the case when the incident puls 
is sharp-fronted and of constant unit pressure, and graphs drawn to illustrate tl 
way in which the mobility of the strip affects the pressure distribution on the bac} 
of it. 


1. Introduction 


THE problem considered is the two-dimensional one of an infinitely long 
strip of finite width subjected to a plane-fronted sound pulse at norma 
incidence. The strip may be regarded as the obstacle produced when a wal 
perpendicular to the ground is ‘reflected’ in the plane of the ground. 

For an incident sinusoidal wave-train there are the classical solutions 
of Rayleigh (1) for the extreme case of wave-length large compared wit! 
the dimensions of the obstacle. Other more recent investigators in this 
field include Copson (2), Sieger (3), and Morse and Rubinstein (4). In 194 
E. N. Fox (5) obtained a solution of the stationary strip problem with a1 
incident pulse of a general two-dimensional form, and carried out compu 
tations in the particular case of a sharp-fronted pulse of constant uni 
pressure. 

In this paper, based on this solution by Fox, the two-dimensiona 
problem is investigated when the strip is capable of motion as a spring 
supported, rigid body. By a general principle propounded in Fox’s pape! 
to find the pressure at any point it is sufficient to obtain the pressure 01 
the back of the strip, since these two quantities are directly related. 

Adapting results obtained by Fox, the exact solution for the pressurt 


on the back of the strip has been obtained for all moments up to a time 4! 


after the onset of the pulse, where 2/ is the width of the strip and a is the 
velocity of sound. Calculations have been carried out on the pressure 0! 
the back of the strip up to a time 2//a, and on the velocity of the strij 
up to time 4//a. 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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No calculations have been carried out on the pressure at points off the 
strip, but, as Fox points out, it is doubtful whether such calculation would 
rield any new results of appreciable physical interest. In effect, points on 
the back of the strip are in the deepest ‘shadow’, and thus the results 
\btained indicate the most pronounced effects of diffraction to be expected 
in the problem. 


2. Strip moving as a rigid body 

Let us now consider the problem in which an infinite strip of finite 
width is struck at normal incidence by a plane-fronted pressure pulse. 
For simplicity in writing formulae we make the system non-dimensional by 
taking the width, 2/, of the strip as the unit of length and the time, 2l/a, 
taken by a sound wave in travelling this distance as the unit of time. The 
wave velocity is then unity and we may write the incident pressure pulse 
1s pt+-y), where time is measured from the instant at which the pulse 
strikes the strip. 

In the (x, y)-plane the strip is represented by the straight line y = 0, 
}<x < 1 when at rest. The strip is considered as being a rigid body, able 
to move in one-dimensional motion parallel to itself, as if supported by a 
spring with damping. 

Let the total pressure set up be p(x, y,t) and write 

Prl&, y, t) Polf+y)+ p(x, y, t) (1) 
where p(x, y,t) is the diffracted pressure. 

Taking the corresponding velocity potential as ¢, and the velocity of the 
strip at any moment as U(#), then for a point on the strip we have 

og / Cy U(t); 


I, a We are assuming small oscillations, 





¢ ,; : 
. d U(t) approximately. (2) 
cy 
« y=0 
Hence, eliminating ¢ between (2) and the equation p, = —p é¢/ét, 
where p is the density, and substituting from (1), we find 
Cp dU dpy ( ) ) 
_ : za 
cy p dt dt 





y=0 

Now let us compare this with the postulate of Fox’s stationary strip 
problem, when the incident pressure is [pU(t+-y)+po(t+-y)]. Denoting the 
total and diffracted pressures thus obtained by p(x, y,t) and p(x, y, t) 
respectively, then the appropriate boundary condition is that 


a 0. as Ea ae dU dp, 
4 y=0 





| ey cy Prat dt’ : 


y=0 
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Thus we see that the present problem and that of Fox are mathematical} 
identical, and we may write 
p(x, y,t) = p(x, y, t), 
or, pra, y,t) = pP(x, y, t)—pU(t+y). (3 
2.1. First interval: 0 <t <1 
Fox showed that in the interval 0 < ¢t < 1 the contribution to pP(a, —0,1 
from the edge x = 0 is 


t-—z 
or pret pel) || x }) dz H(t—z). 4 
W 3 t—z t—xr—z 
0 
where H(t—z7)=1 ift>z2 (4a 
20 ift<-z. 
The contribution from the edge x = 1 is 


pp(x,t) = p,|(1—2), t]. 5 
Thus, in the interval 0 <¢ < 1, 
Pr(x, —0,t)—py(x, +0,t) = 2[p,+p.,—pU(t)—p,(t)]. (6 
We can now formulate the equation for pU(t), since, for the motion o/ 
the strip, we have the equation 
dU : - : 
m ate +k | U(z) dz = | | pr(x, —0, t)—pp(a, +0, t)] dz, 7 
0 0 
where m, c, and k are constants connected with the mass, damping, and 
spring stiffness involved. 
On substituting for the integral on the right-hand side of (7) from (6), it 
is found that the integration can be carried out, after changing the order 
of integration, and thence we find 


, t t 

dl ; be : 

m Tt +(c+2p)U+(k—2p) | U(z) dz = 2 | Po) dz pot) 8 
0 0 

This is the equation of a forced damped oscillation, the precise natur 
of which will be determined by the relations between p and the constants 
m, c, and k. 

Assuming that U(0) = 0 and that the roots g = a and q = B of thi 
equation ma? \ (c+2p)q Lf: 2p 0 
are distinct, equation (8) can be solved to give 

t 
; 2 f 2 1—B 
U(t) ~ | Pol(t—z) mgt ee P eB: dz. \v 
m J x—B x—B 
0 
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natical The pressure on the back of the strip is found by inserting this expression 
for U(t) in (4) and (5) and applying (6). 


2.2. Second interval: 1 < t 2. 


Fox showed that in the interval 1 < ¢t < 2 the contribution from the 
edge x 0 to p(x, —O, t) is 
(xv, —(0) Pg(%, b) p(x, t)— F(x, t), (10) 
; Ll f p,(l+y,t—x—2y) [x\} 
where F(a, t) fe y ”( } dy, 
4 7 uy 7] 
(4 F(x. t) 
1 2 7 2y 1 ; 
r Ll {x\} i Zz) z 1+-2 4 ‘ 
_ ( Py | # (2) + Pol) Y__\? dz H(T—2), 
Tr c+y\y 4 |-+-7'—zx -y¥—zZ T —xz- 2y—z 
0 0 (11) 
ind T' (¢—1) 
The contribution from the edge x | to p(x, —0, t) is 
) p4(x,t) = pal (1—z), t]. (12) 
notion 1 
In the course of the solution we shall require | p, dz; let us first evaluate 
0 
x. - p, dx. If we write 
w(z) pU(z)+p,(z), (13) 


ping, a0} then, after performing a change of order of integration and evaluating the 


inner integral so produced, we find 






om (6), it ' ‘ 7 
- il : if ic c T—2)! 
the ord p(x, t) da w(z) dz | w(z)| tan 1(7 —2z)? + ( — : dz. 
2 wT 14-7 —z 
0 0 0 (14) 
1 
Dealing in a similar manner with | F(x,t) dx we find 
0 
se natul 1 7 
' _— wi = T'—z)* e 
constants F(a, t) da w(z)|tan—1(7 z+ - | de. (15) 
! 7. 1+ 7'—z 
0 0 
p ot tl Hence. 
| Prl2, 0, 7) Prl*, 1() f) da 
T 1 
a 4 sil , 12)! 
2 w(z) dz—wi(t) w(z)| tan-1(7’'—z)! + Li .. dz\, 
} F wT T'4+-l—z } 
0 0 
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and this can now be substituted in the equation of motion. If at the same 
time we put 


t 
: : 2 I—a ., 1-8 2 
U,(t) = U(t)— [ rote—=)| exr_ 1 cl dz, 
ak r—B a—B 


we find that the equation satisfied by U;(t) is 


; t 
dU, | sii -_ 
m di + (c+2p)U,+(k—2p) | U,(z) dz 
0 
jag : 
ate t—1—z)¥ 
=-— | [ pl (2)+p9(2)]| tan 1(t—1—2)$ +! — ) dz, (16) 
i 0 ‘ 
where the right-hand side is known, since it involves U(z) for 0 < z < | 
only. 


This can be solved by Heaviside operational methods, and we find 


t 
» fe man _— 
U(t) — | palt—2)|. Sn” 30 | ee 
m J 


r—f r—B 
0 
t-1 t—y-1 
8 - y? e - ” 
_ “ __ 4 tan—lyt| da | 1,(t—y— 1—z)| we™ — Be] dz— 
mr( x —f) | F Ly d | zt Polly / ; P | 
0 0 
t-1 t-—y-1 
16p rT y! : 
= ~ __ittan-lyt| da | ».(t—y— 1—z) x 
m®n(a—B)3 | ; 7 | yo Polt—y ) 
0 0 


J g 1—a«a)ze®-L B(1 — wpb: 1. 2>B—a—B e%2__ pz dz. (17 
N | ( )z | B( B)z x—B [ | 


3. Numerical calculations 

Calculations were carried out for the case of an incident pressure pulse 
H(t+-y), where H(t) is defined as in (4a). Also c and k have been taken 
as zero, to simplify calculations, justification for this lying in the fact that, 
from mechanical considerations, the inertia effect will be predominant in 
the early stages of such an oscillation. 

During the period 0 < ¢t < 1 the velocity and the pressure on the back 
of the strip for x = 0-1, 0-3, and 0-5 have been calculated at intervals oi 
0-1 in time. During the period 1 < t < 2 the velocity, only, has been 
calculated at intervals of 0-2 in time. The expression for the pressure on 
the back of the strip was complicated and any calculation on it would 
hardly have been justified. 
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The interval O <t< 1 
Taking c = k = 0 and po(z) = H(z) we find 


: l 2 l—« 1—B . 
U(t) —-—+ (6 ct, (18 
Pp m x—B) B ) 


x 
where , ptis’+Semr p= —p—(p*?+2pm)! 
m , pr ’ 
and p 8/3 » (density of air), 


m = 41? x (thickness of wall) x (density of wall). 


Also 
H(t—z) 2pf1— |] — 
D, no | J 118)], (19) 
T1(xa—pb) m x B 
v(t—x) . 
where I(x) Ze: t-a)yh | ll - du, (20) 
u*+ax 
0 


this being in the form computed. 


Calculations were carried out for four values of p/m as follows: 





( Case 2 Case 3 Case 4 

1/40 1/4 9/8 81/20 

X : ; 1/5 1/2 | 3/4 g/10 
1/4 I 3 9 


Case 1 corresponds to a brick wall 9 inches thick and approximately 
l4 feet high. The other cases are perhaps more likely to occur in submarine 
blast problems. 

The interval 1 <t < 2 

It is found, on substituting for p,(z) and using known relations between 
. and £8, that 








» —_ 
pU(t) 14 =P * w_} B ape _— 
m(a—B) x B 
3 
ee \ 2 | 
: 1 tan-ly? l—a)t—1—y) 4... ——__ J e**-1-9). 
m*n(a—B) . \l+y y }\ ( oH y)4 x— B 


t —B)(t—1—y)—- : 


]em--0) dy. (21) 


The values for pU(t) corresponding to case 1 were found to be so small 
in the first interval that they were not computed in the second interval, as 
they would not have exhibited the fluctuations due to the additional term 
sufficiently well. The following is a table of values of (—1)pU(t) for the 
interval 0 < ¢ < 2, and displays the oscillatory nature of the motion. As 
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there are no spring or damping forces involved, the strip would continue 
to oscillate slowly back to the rest position. 





t Case I Case 2 Case 3 Case 4 | t Case I Case 2 Case 3 Case 4 

° 0000 0000 0°000 0°000 od 0°024 0°203 o°56 o'812 
orl 0005 0°046 o*192 o*516 o'9 0°024 0°206 0°553 0°795 
o'2 0°009 0086 0°329 o°74I 1°O 0°025 0°205 0°5358 0°76 
o°3 O'ol3 o'1Ig 0°424 0°320 I°2 O°197 0°506 0°772 
o*4 o'016 o'146 0°459 0845 I°4 o19g2 0°503 0°797 
o's o’oIg 0168 0°530 0847 1°6 o'1g2 0°520 0817 
o°6 0'021 o'184 0°554 0°840 1°38 0*200 0550 0°836 
o°7 0°022 0°196 0°564 0°828 2°0 0°216 07589 0°857 








Figs. 1 and 2 are graphs of the pressure on the back of the strip plotted 
for the interval 0 < ¢ < 1, and contain four sets of curves labelled 1 to 4, 
corresponding to cases | to 4 respectively. Each curve labelled A is for 
the point x = 0-1, B for x = 0-3, and C for 2 = 0-5. The discontinuities 
in slope displayed by each curve represent the arrival of diffracted wave- 
fronts from one of the edges. 

The set of graphs taken together cover the transition from the case of 
a very heavy (i.e. nearly stationary) strip, for which p/m is small, to that 
of a very light one. Asis seen from the graphs, the pressure, [ p, + p.— pU (t)], 
on the back of the strip is greater for the lighter strips, which means that 
these would be less effective as shields. However, the magnitude, 

2[1—py—Po+ pU(b)], 
of the pressure discontinuity across the strip is less for the lighter strip, 
which means that the total force acting on such a strip would be smaller. 

Whereas in Fox’s problem the pressure on the back of the strip did not 
become equal to that of the incident pulse, i.e. unity, until time ¢ = 1, we 
see that, in all the present cases, such ‘equalization’ takes place before this 
time. (At a time when equalization takes place, the sign of the pressure 
difference across the strip changes.) The following table shows the approxi- 
mate time at which equalization first takes place for the various curves: 





Curve C . e 0°98 0°34 


Case 1 Case 2 Case 3 Case 4 

Curve A . ‘ 0998 0°95 0°93 o18 
Curve B . ‘ 0°98 0°87 0°75 0°35 
0°59 o's 


Sets 3 and 4 exhibit definite maxima, and set 4, in particular, shows 
some rather striking features. Firstly, equalization occurs very soon after 
the arrival of the first diffracted wave. Then, after attaining a maximum 
of between 1-0 and 1-2, the pressure on the back falls off until it is below 
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unity, rising again above unity soon after the arrival of the second 


diffracted wave. Thus the strip is being subjected to a force of frequently 


Pressure 








0 Of 02 03 04 O05 06 OF 08 09 10 
Time 
Fie. 1. 
A 





varying sign, and it seems likely that this phenomenon would persevere 
throughout the motion. 

Thus, summing up, the heavier strips are subjected to forces of greater 
magnitude and of mainly constant direction, while the lighter strips are 
subjected to smaller forces but of more frequently varying direction. 
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The method which has been used could, of course, be extended to cover 
a longer period of time, but the additional labour involved would be greater 
than the result would warrant. 

In conclusion, the author wishes to express his gratitude to Dr. F. G, 
Friedlander and Mr. D. S. Jones for all their help and encouragement in 


connexion with this paper. 
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SUMMARY 

A treatment is given in this paper of the two-dimensional diffraction of a plane- 
fronted sound pulse, normally incident on a semi-infinite plane screen, each particle 
of which is able to yield under the action of the incident pressure. Two methods of 
finding the pressure change across the screen are given, both depending on the method 
of successive approximations. Graphs are given illustrating the results found by the 
first method for several types of screen material. 

1. Introduction 

Tue problem considered in this paper is the two-dimensional one of the 
semi-infinite plane screen subjected to a plane-fronted sound pulse at 
normal incidence. The solution for the stationary screen is, of course, well 
known and has been discussed by many authors. The stationary screen is, 
however, by no means an accurate representation of that employed in 
practice, as such a screen is neither immobile nor rigid. Although it is not 
possible to treat the practical problem as it stands, some characteristics 
can be reproduced by allowing the screen to move in ways which are 
approximately like those observed. 

In a previous paper (1) the problem considered was that of diffraction 
by an infinitely long strip of finite width, capable of motion as a rigid, 
spring supported body. Such an obstacle, although rigid, still reproduces 
one feature experienced in practice, in that a disturbance at any point 
immediately affects the entire screen. In the present paper another aspect 
of the screen used in practice is considered, namely that of non-rigidity. 
Now, each particle of the screen is supposed capable of yielding indepen- 
dently under the action of the pulse. 

This may be considered as giving some indication of what occurs when 
i sound pulse strikes an infinitely long strip, clamped along the middle, 
ind able to move according to the laws of flexural vibration. Firstly, 
neglect of the term due to bending in the equation of motion of the strip 
will be expected to be a good approximation near the edges of the strip, 
as there the bending moment will be nearly zero, only becoming large near 
the middle of the strip. Also, in the interval of time between the incidence 
of the pulse and the first arrival at one edge of a diffracted wave from the 
other, the effect can be treated as a superposition of waves arising from 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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semi-infinite plane screens. Thus, in such an interval the solution obtained 
in this paper may be expected to be a fairly good guide to the state of 
affairs near the edge of the strip. 

An integral equation is set up for the pressure saltus across the screen, 
This is then solved by two methods: firstly by a direct application of the 
method of successive approximations, and secondly by applying the same 
method to a modified form of the equation. For the purposes of computa- 
tion the first method is found to give the more satisfactory result, since. 
although the series so obtained is not so rapidly convergent as that 
generated by the second method, the terms are more easily calculated, and 
there is the additional advantage that one set of computations suffices for 
a fairly wide range of screen materials, whereas the second series has to be 
recalculated for each separate type. 

The first three terms of the series have been computed when the incident 
pressure pulse is sharp-fronted and of constant unit pressure. Graphs of 
the pressure saltus are given for four types of screen material showing the 
modifications introduced by allowing the screen to yield. 

2. Formulation and solution of the problem 

Let us consider a semi-infinite plane screen represented in the (wv, y)-plane 

when at rest by y=0 (0<2x< 0). 


Let the total pressure set up when the screen is struck at normal incidence 
by a pulse p,(t+y) be p(x, y,t) such that 
p(x, y, t) Polt+-y)+py(@, y, bt), (1) 
p,(x, y,t) being the ‘diffracted’ pressure. The unit of time has been chosen 
so that the velocity of sound is unity, and the pressure, p, satisfies the 
equation ap ep ap ; 
att yi = BE @) 
Let v(x,t) denote the velocity of the particle of the screen which is a 
distance x from the edge at time ¢, displacements being assumed small. 
Then the boundary condition with which we shall be concerned is 
m nto = [p(x, y, t)}" pp = p(x, —0,t)— p(x, +0, t), (3) 
ot y ) 
where m is a constant related to the density of the screen. 
Also, by equating the acceleration of a particle of the screen and that 
of the sound wave, we have the relation 


Vv L{¢ ] - 
ot Pp oy y=0 Pp P 
where q,(z,t) = E plod) ; (5) 
cy y=0 
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nd p is the density. Using (1) and then substituting in (3) for év/ét we 
ybtain y= +0 , . 

pl Pp m| qy(x, t) | po(t)|. (6) 
Now according to Friedlander (2) 


son y 


qy(x',t’) 
~ | [ (t t’)? <* 


Py(X, y, t) 


, dx'dt’, (7) 


‘\9 


a’)?—y?| 








‘a tha 
ra i 
| A " / 
WZ ——— 
| P4 a \ 





\4 








I xr 
Fic. l. Fic. 2. 
Letting y 0 in this formula, we obtain 
( ee i q,(2’, t’ een 
cae | | WNT) da'dt’, (8) 
T (t{—t )*—(x—av )?}! 


vhere the area of integration is now defined by 


0 —_— y= t—|a—z’|, 


for f > ax, and will be while for ¢ < 2 it is the area 


t’ t rv H 


denoted by A+B, 
, which will be denoted by C (see Figs. 1 and 2). 

For the remainder of this section, we shall only be concerned with 
v being obtained otherwise in section 4. 


times f¢ r. the solution for f 


Eliminating g,(x,t) from equations (5) and (8), and writing A for 2p/m, 


ve find that 
0 9 a , , 
© dx'dt’ A {| pol! ) dx'dt’. (9) 


1 0, A er L Pil, 
Pil, Ss ee | | zs sa y 


A+B A+B 


where y ({—t’ )}?@—(x—2’)*. 
Here it may be noted that A 0 is the rigid screen case, and then (9) is 


equivalent to the equation set up by Fox (3). 
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This equation can now be treated by the method of successive approxi- 
mations, i.e. we put 


ae A es a, . 
[pil = $o—-—$1 +592 T --+(—3} $,+...; (10) 


7 


where, for all 7, ¢, does not depend on A. 
On inserting this in (9) we obtain 


9 fal ad , “4 
do = | | Pol ) dxdt’, (11) 
ajJJ y 
A+B 
, , bla’ .t’ ne 
and for r >0 br44= [ | P,(2 dx'dt'. (12) 
J, 


By performing the necessary integration, we find 


—_ 
by 2pa(t)—= | mA, = 7)" dt’ H(t—2), (13) 
0 
where H(t) is the Heaviside unit function defined by 
H(t) lL 2et>0 (13a) 
2) we<— 0. 


Thus, when p,(t) = H(t) 
4 t—a\} 
be = 2H(t)—“tan-1| | A(t—2). (14) 
Cs \eae 
This is the term which arises in the theory of diffraction by a stationary 
strip during the same interval as we are now considering, being, in fact, 
a combination of solutions for the stationary semi-infinite screen (see Fox’s 
paper (3)). 
In this case it is found that we can perform one of the integrations 
involved in finding ¢, and that we can write 


do, = 2ntH(t)—2n(t—a)H(t—2) 


1 
-\ 1 )», ) ‘- 
4H (t 2)fttan ( set x) | (1—u2)! tan (< du— 
x 7 Uu 
0 
1 
9 ¥ 4/ D__9a,2 l 
— —(t—z) | (e+2—v?)! tan-! nilirccueens mis du| (15) 
7 7 u\e+2—u? 
0 


where « = 22/(t—2z). 
When these integrations were carried out (by numerical methods) it 
was found that 


d, = n| 2tH(t)+ (a—2)(t—x)H(t—2)], 
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proximately, where a is roughly —1-6 and varies slightly with x. This 
was used to obtain an approximate value ®, for ¢,. Thus we find that 


17-| 2(5t—2x)ah(t a)'+-4(a—2)(t—a)(2x)*(t—ax)!|H(t—2). (16) 


_ { Y 1 - 
©, = 2nt?|}— H(t)—tan H(" “ )!| Hea] + 


[It can be shown that the series is convergent, but the proof is omitted. 

If, for the moment, we consider the strip problem, recalling that the 
lution just obtained will be a good approximation near the edges during 
the interval in which the effect can be considered as a superposition of 
two semi-infinite screen solutions, we may thereby obtain some idea of the 
rror involved as follows. 

In the (x, y)-plane let the strip extend from 2 0 to x :. 

The equation of motion to replace equation (3) would be 


m LB r" v(a,t’) dt’ [ p(x. y,t)}) ies (17) 


vhere m and B are constants, and, of course, further restrictions on v and 
ts derivatives are imposed by the support. 

From (8), (10), (11), and (12) it can be seen that the expansion (10) can 
e rephrased in terms of q,, and is then equivalent to setting 


, A 
q(x, t) pald)+0| } 


_ 
7 


Using this as a basis for estimating the bending term on the left-hand side 

17) (in a purely formal manner), we see that it gives rise to a term of 
der O(AB/pz), which is a second-order effect. Thus, to the first order at 
east, our approximation should be valid. 


3. Calculations 

The quantities dy, 6,, and ®, as defined by equations (14), (15), and (16) 
were computed for values of 2 from 0-2 to 1-0 at intervals of 0-2, and of t 
irom 0 to 2-0 at similar intervals. 

Taking this scale for x and ¢ necessitates care in fixing the value of A to 
nake the computation conform to a physical system. The actual value 
required is 

A = 21 (density of air) x (density of screen material)-! > 

< (thickness of screen)-!, 
where / is the equivalent of the unit of length in centimetres (e.g. if we were 
lealing with the strip, 2/ would be the actual width of it). 


5092.19 
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Hence the expression db, , 2, 
do— —- 

7 7 


was calculated for A = 1/40, 1/20, 1/10, and 1/4. It was found that fy 
A = 1/40, the third term was making a maximum contribution of 0-00) 
while for A = 1/4 it was making a contribution of 0-100. 

Figs. 3 and 4 are graphs illustrating the results obtained, and contaiy 
four sets of curves, labelled 1 to 4 to correspond to the cases A = 1/40 
1/20, 1/10, and 1/4, respectively. Each curve labelled A is for the point 
x 0-2, B for « = 0-6, and C for x = 1-0, curves for intermediate points 
being omitted for the sake of clarity. 

The effect of the motion of the screen is already quite marked in the 
set of curves 2. Curves 1 show very little deviation from the corresponding 
curves for the stationary screen, as they represent a fairly weighty wall- 
the value of A taken being as for a 9-in. thick brick wall about 15 ft. high, 

A qualitative comparison with the corresponding results in the previous 
paper (1) shows that the pressure discontinuity here is greater than that 
previously calculated, except during the period of time which elapses 
between the arrival of the pulse and an instant very soon after the arrival 
of the diffracted wave, when the reverse is true. This means that, after 
this initial interval, the new type of screen is more effective as a shield 
although the force on it is correspondingly larger. 

A similar comparison amongst the sets of curves given here reveals a 
more complicated state of affairs. The set of graphs taken together repre- 
sents a survey of the transition from the case of a very heavy screen to 


a much lighter one. Now, near the time ¢ = 2, the pressure discontinuity 
is larger for the lighter screens, except near the point x = 1-0 (curves C)— 


which, for practical purposes, is the least important, being too far from 
the edge of the screen to be a good representation of the corresponding 
point on the strip (the mid-point). However, near the beginning of the 
motion, and for some time after the arrival of the diffracted wave, the 
pressure discontinuity is smaller for the lighter screens. The instant at 
which one condition gives way to the other varies appreciably as we make 
the transition from the heavier to the lighter screen. Thus we can only say 
that the heavier screen would be the more effective as a shield in the 
earlier part of the motion, while later the lighter screen would be the 
more effective. 

No computation has been attempted for points off the screen, but it is 
doubtful whether any new results of appreciable physical interest would 
be yielded, as it is at: points on the screen that the strongest effects are to 
be expected. 
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4. An alternative approach 
An alternative method of dealing with equation (9) suggests itself when 
this equation is written out for? < x. It is then 
0 


y=4 * Pe 
[pil}- +2 ff Pale 0 dx'dt’ = = [| Pol! ) dxdt’. (13 
Id 


# Y vw. i 
C 


We note that the area of integration now involved is C, and this is precisely 
the area of integration involved, whatever the time ¢, for the case of norma 
reflection by a yielding infinite plane screen. 

The physical interpretation of this is that the effect which diffractio: 
makes upon the screen does not travel along the screen quicker than 
through the air. In other words, a point on the screen does not ‘know 
of the existence of the edge until disturbances set up there have had time 
to propagate to the point in question through the surrounding air. 

In consequence of these remarks, let us consider for a moment the 
problem of a plane screen of infinite extent, represented by the plane y = 0, 
at rest, when the incident pressure field is P,(¢+-y) giving rise to a total 
pressure P(x, y,t) such that 

P(x, y,t) = P(t+y)+P (a, y, b), (19) 
where P, is the ‘reflected’ pressure. 

Suppose that we had a boundary condition 


oP 
. =AXAP 
= y=0 | ly on 


ie ALPiy +071 AP,(t)— Po(t). (20) 


y=0 


° f P, 
i.e. 
3 CY © 
Again using Friedlander’s result, 
9 PPTAaAD/|A 
epee 2 OP, /ey 
[Ply —0 Pil, = | a 
mJ. y! 
A 
Eliminating [éP,/éy]|,,-9 from (21) and (20) we have 
gar Py=*° a eine Kae 
[A], et : [| | uy, Sd’ dt! = == = )— Pot 
7 aye y! 


WT 
C ( (22) 





Y=0 dy'dt’. (21) 





)| da'dt', 


Jv @ / wa 


and this is precisely the same equation as (18) if we put 
APj(t)— Po(t) = —po(t). (23) 
Now, since this problem is one of pure reflection, we can assume that P 
will be of the form F(t—y) for y > 0, and then applying (20) and (23) we 
find that 


t 
F(t) polt)—A | po(t—z)e™ dz. (24) 
0 
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Now returning to equation (9) and using the value just obtained for 


py) : when ¢ < x, we may write 
4 a) ld Y +0 ; 9 an ’ t’ ; , 2 F’ t’ , . 
" ir 1 : ly 9 dxdt | Pol! dz'dt’ — — ) dx'dt', 
i TJ. r Je 7 Y 
1 A B 
nd hence, carrying out the integral over B, 
Aff [py sili ne & H(t iF 
7JJ y aJJ 7! 
A 
t 2 (t—zx)/2 
2 ¢ Fit’ t—x 2 Ft x : 
an ( ,) at’ — ts \ dt’. (25) 
7 J t t AJ t HX 2t 2 t—t t—x—t 
t—x)/2 0 


The double integral on the right-hand side has not been evaluated for 
reason which will become apparent later on when a special value is taken 


for po(t). 


We now make a change of variable by putting 


t+x 20 t’ +2’ 2a’) . ' 
| throughout and | in the double integrals. 
t—2zx yA w) H o 
ae ‘ > . / 1 +0 > 
We also write Y'(a,B) for [ p,fa(a«, B), t(x, B)} I) aif (26) 
Then we have 
* ¥( (x’, B’) some 
‘Y'(a, B)- = | he 7 any, da’dB o(a, B), (27) 
(p p )|! 
0 B 
vhere 
(a’-+") 
J | 
o(2, B) Wr Pol: Fi da'dp 
(B—p’) 


B 
2 Fit) B \} 2 F(t) [«a—B\} 
ex o It. 28 
E 3 A ) = | edge) (28) 


For the sake of simplicity, let us consider the special case of p(t) = H(t), 


this being the one usually used for computation. Then the double integral 


er the area A on the right-hand side of equation (28) vanishes since 
t') = &(t’) 0 unless ¢’ 0, and this lies outside the area A; and yy 
ecomes 


>F eM (1 B\L,, 2f e* [aB\s 
Po(a, B) : | }" dt —— | : (29) 
77 r+ B—tla—t, a) atp—t 2p—t 


RQ 0 
p 
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We now put 


H(a,B) = Yo(a,B)— HP) sg (CAB) 


r 


(30 
7 7 
where, unlike the ¢,, the ys, are now functions of A as well. On substituting 
in (27) we obtain 


ez - 
w44(a, 8B) = [ [ Pirlo B’) — da'dp’, (31 
JJ [(a—a’)(B—B’)}! 
0B 
for r > 0. Of course, the same expansion still applies when yf, is given by 
(28), but then the actual integrals which appear always involve one more 
integration than in those appearing when yf, is given by (29). 
The series (30) can be seen to be rapidly convergent as follows. It can 
be shown by iterating equation (27) that 


Be 
s,+9(a, B) | | os.(a’, B’)k(a, B, x’, B’) da’d’, (32) 


where 


1 r , 
KasBua’s pt) = 4sin-{ 2) sin” yp) for B <a! <p 
_— . 


x B—Pp’ (33) 
7 for B < a’ : J 
Then, if we assume that || < M (some constant) for some range of « 
and f, application of (33) shows that 
— Ap" 8 ,; 
> = = lebo,| u| Iy{2Ay (aB)} _t ~ ot 2A, (aB)} | (34) 


n=0 


where /,, J, are Bessel functions. 
Also, taking || < 4Ma [it is actually < 4MB(a—B)!] it is found that 





A2n } 1M 2 ; 
S 2n41 th. "te (; } | [,{2A\/(«B)} wale T,{2d,/(oB)} (35) 
oa a \B < 


The factor e- in equation (29) makes the series all the more rapidly 
convergent. Unfortunately, this is also the chief disadvantage of the 
method, as, for purposes of computation, the integrals in (29), for example, 
would have to be recalculated for each separate value of A, whereas, in the 
method adopted in section 2, one computation is sufficient for a wide range 
of values of A. For this reason, no computation has been attempted on 
this series. 





In 


Mr. 


i) 








riven by 


he more 


It car 


ge ol 


nd that 


rapidly 
of the 
‘ample 
, in the 
e range 


ted or 














DIFFRACTION OF A SOUND PULSE 343 


In conclusion, the author wishes to thank Dr. F. G. Friedlander and 


Mr. D. 8. Jones for all their assistance in connexion with this paper. 
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SUMMARY 


In this investigation we employ a combination of complex potential and Fourier 


OF THE ELASTIC 


integral methods to solve problems of generalized plane stress of semi-infinite plates 
with one straight boundary. In the first part of this paper relations between th 
complex potentials (Q(z), w(z) are obtained suitable for giving any desired conditions 
of stress or displacement along the real axis, thus reducing problems of the half 
plane y > 0 to the determination of the single potential function Q(z). A simpk 
method is given of obtaining (Q(z) from specified conditions along the real axis and 
at infinity. The method is applied to the solution of problems of the half-plan 
y > 0 in equilibrium with specified stresses or displacements along the straight 
boundary and balancing stresses at infinity and to the problem of the half-plan 
subjected to an interior force and with the straight boundary free from displacement. 
Introduction 
For the majority of two-dimensional elastic problems, solutions in terms 
of complex potentials are simpler and more informative than Airy stress- 
function solutions. However, a difficulty in the use of complex potentials 
lies in the lack of general methods of obtaining them. For the Airy stress- 
function treatment the results of years of research into the solution of 
differential equations is available. Whilst the tentative method has proved 
powerful (see Stevenson (1), (2)) it is important to establish general methods 
of solution in terms of complex potentials. For problems involving material 
bounded by a simple closed contour, an elegant method for the determina- 
tion of the potentials in terms of given boundary stresses or displacements 
has been given by Muschelisvili (3, 4), who has also given a treatment of the 
elastic half-plane. The publications (5, 6) containing the latter investiga- 
tions are not easily obtained. However, Muschelisvili’s basic assumption 
that the straight boundary is stress-free except for a finite number of finite 
segments shows that the treatment is not exhaustive of all possibilities. 
Sneddon (7, 8) and Girkmann (9) have considered the problem of the half- 
plane with specified stresses over the straight boundary in terms of the Airy 
stress function. The treatments are purely formal and introduce unnecessary 
analytical difficulties. 
Fundamental equations 

The following investigations are confined to elastic material in equili- 
brium, subjected neither to body forces nor body couples, and in a state 0! 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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yeneralized plane stress. The stress and displacement components may be 
expressed in terms of complex potentials Q(z), w(z) as follows (Steven- 


son (1, 2)) 20 = 2(%%+7y) = O'(z) +0’), 


(1) 
20 = —2(zx—Yy—2ixy) = 20"(z)+w"(z), (2) 
8uD Su(us wv) KQ(z) —202'(Z)—a' (2), (3) 
where « is an elastic constant and Q(z), w(z) have continuous derivatives of 
the second order at all points in the interior of the material. From (1); (2) we 
see that 
2(EQ— MD) 4( yy +-ixy) Q’(z)+-0'(2Z) +20" (z)+w"(z). (4) 
Notation. Throughout this work considerable advantage is gained by 
choosing the straight boundary of the material as the real axis and we shall 
write 


PIy-9 = pq etc. (5) 


Relations between the complex potentials suitable for giving any 
desired conditions of stress or displacement along the real axis 


i) Let w/(z) 2Q(z)+2 [ Q(z) dz. (6) 
Equation (4) becomes 

4 (Hy + ixy) = Q'(z)+0'(Z)—2t1yQ"(z). (7) 

Thus 27y° = re[Q’(z)],-9 = re[Q’(zx)], xy = 0. (8) 


We asume, in general, that lim yQ”(z) = 0 at all points of the real axis. 
We may relax this condition to include those cases in which this limit exists 
but is not zero or unique at a finite number of points of the real axis. From 
8) we see that this combination gives zero shear over the real axis and if we 


require wy = f(x) (9) 
we must choose Q(z) so that 

re[Q’(x)] = 2f (x). (10) 

ii) Let w(z) = —2Q(z). (11) 

Equation (4) becomes 

4(yy+izy) = Q'(Z)—Q'(z)—2iyQ”"(z). (12) 

Thus yy = 0, 27y° —im|Q’(x)]. (13) 

If we require ry? = (x) (14) 
we must choose Q(z) so that 

im[Q’(x)] = 24(2). (15) 


(iii) Let w(z) —20(z)+(1—k) [ Q(z) dz. (16) 


. 
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Equation (3) becomes 

8uD = n{Q(z)+-O(zZ)|—2iyOQ'(Z). (17) 

Thus 4uu° = cre{[Q(x)], y 0. (18) 


We have assumed that lim yQ’(z) = 0 at all points of the real axis but this 
may be relaxed as above. 


[f we require u® = (2), (19) 
we must choose 02(z) so that 

re{Q(x)| = * yb(x). (20) 

(iv) Let w(z) -2Q(z)+(1+«) [ Q(z) dz. (21) 

Equation (3) becomes , 

8uD = f{Q(z)—O(2)|—2iyO' (2). (22 

Thus wf” = 6, 4uv® = xim[Q(a)]. (23) 

[f we require v = x(x) (24) 
we must choose (Q(z) so that 

im[Q(x)| -t x(x). (25) 


Equations (10), (15), (20), (25) reduce problems of the half-plane y > 0 with 
specified stresses or displacements along y = 0 to the determination of 
functions which are analytic in the half-plane y > 0, of suitable orders of 
magnitude at infinity, and which have specified real or imaginary parts on 
the real axis. We will consider, in general, problems in which the potentials 
at infinity are to be of orders such that 

O'(s) = Of2-*), w"(z) = O(z-). (26) 
Thus the stresses at infinity are O(z-!). These are the lowest possible orders 


if the stresses applied along the real axis have a non-zero resultant. 


Determination of the function F(z) which is analytic in the region 
y > 0 and has a specified real or imaginary part along the real 
axis 
Clearly we need consider only the case for which the real part is given 

on the real axis. 

Let f (2) satisfy conditions sufficient for it to be expressible as a Fourier 
integral (10), e.g. let f(x) be of bounded variation and at each point 


f(x) = Sf (w#—0)+f (a +0)] (27) 
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(‘onsider the function 


x 


, :, ae 
F(z) G(x,y)+iH (x,y) = Sr(uje®" du, (29) 
7 
where Tru) ( f (tje-™ dt. (30) 


It can easily be shown that, subject to (27), (28), F(z) has the following 


properties: 


(a) G(x,0) =f (x) (see (5)). (31) 
(6) lim G(x, y) f (a). (32) 
(c) F(z) is analytic in y > 0. (33) 


(d) F(z) 


; t)¢ 
&. ) at in y > 0. (34) 


If, in addition to (28), we assume that 
( vf (a)| dx < o. (35) 


(e) F(z) = O(z-) at infinity in y > 0. (36) 
Thus, in general, F(z) is the required function which is analytic in the 
upper half-plane and has real part f (2) on the real axis. Thus equations 


10), (15), (20), (25) are satisfied respectively by choosing 


9 Pf 
Q2'(z) I(z) — | fy(uje®" du, (37) 
“4 
21 , 
Q'(z) = J(z) | buyer” du, (38) 
9 
4u [ : ‘ 
Q2(z) | pp (uje™" du, (39) 
K7T 
0 
tui = 
Q(z) | X(upei™ du. (40) 
K7T 
0 


Applications of the above theory 
The following examples illustrate the simplicity of the above theory. 


( : — i, lx| < 6, 


(i) Given yy = f(x) 
0, |x| > b, 


ay = 0. (41) 
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From (37) 

oo 
I(z) = - [ 2z log z—(z—b)log(z—b) —(z+b)log(z+5)]. (42) 


70 


Elastic material 





External forces 





The integration may be performed most simply by means of (34), whence 
miQ(z) = (z—b)* log(z—b)+ (z +5)? log(z +b) — 22? log z. (43) 
From (6) 
37biw(z) = 8b?z+ 2z5 log z—(z-+ 2b)(z—b)* log(z—b)— 
(z— 2b)(z+b)? log(z+-b). (44) 


At infinity Q(z) = O(log z), w(z) = O(zlogz), so that (26) is satisfied. The 
functions Q(z), Q’(z), Q’(z), w 


rt 
stresses given by (1), (2) are finite and uniform in y > 0 and the displace- 


(z), w’(z) are continuous in y > 0 so that the 
ments given by (3) are finite and uniform in the finite part of the upper 
half-plane but are necessarily logarithmically infinite at infinity. 


| (7-1), lb<a<b 


| 


(11) Given 


to — 2. — 1 << ae < 2}, yy? 0 (45) 
xry® (2) b . vy 


From (38), (34) 


J(z) = : [ (z+6)log(z+6)—(z—b)log(z—b)+ 


oath 
+ (22—b)log(z— 4b) —(2z+-b)log(z+ 4))]. (46) 


The potentials Q(z), w(z) may now be evaluated from (38), (11). In this case 


J (z) O(z-2) so that Q(z) = O(z-4), w(z) = O(1) at infinity. Since the 
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stresses applied over y = 0 are in equilibrium, it might be expected that 
the stresses at infinity would be of a lower order. Uniqueness of the solution 
11) shows that this cannot be the case. 


QS Elastic material 








(iii) Given w® = (x) ——— , r= ®, (47) 
(a—b)?+a? 


where a and 6 are real and positive. From (39) 
4yi 


(2(z) —-, 
an(z—b-+ia) 


(48) 


The second potential w(z) may now be obtained from (16). 


- x—b 
iv) Given Th us(x) bs y? (0. (49) 
(a—b)?+a? 
In this case (2) does not satisfy (28) but may still be expressed as a Fourier 
integral (see (10), p. 16). From (39) 


4 
Q(z) —_ (50) 
«(z—b+ia) 
ind w(z) may be obtained from (16). 
(v) Given u® = &(2) log| (~—b)?+-a?], vy? — 0. (51) 


The function (2) cannot be expressed as a Fourier integral. We assume 
tentatively that 


4 
Q(z) = —log[z—b+ia] (52) 
= 


and w(z) is as given in (16). We see immediately that (20) is satisfied, so 
that (52) is the required potent ial. 
(v1) Given Vind x(x) " u® = 0). (53) 


(2 b)?+ a* 


From (40) Q(z) alia - (54) 
xa|z—b+ia] 


and w(z) is given by (21). 
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(vii) Given v° = y(z) = tis - uw — 0. (55 
, [(~—b)?+a?]’ bi 

. 4yr 
From (40) Q(z) = —_- (56) 


and w(z) is given by (21). 
(viii) Given 
pe x(x) log| (a—b)? +a], yu 0. (57 
We again revert to a semi-tentative method. Let 
Sur js 
Q(z) = loglz—b-+ia]. (58 


If w(z) is given by (21) we see that (25) is satisfied so that (58) gives the 
required potential. 

Solutions (iii)—(viii) will be applied to the solution of the problem of the 
half-plane subjected to an isolated interior force and with the straight 
boundary free from displacement. 

The potentials giving an isolated force F = X+7Y at z = b+ ia in an 
infinite plate are (see (1)) 


Q(z) —vF log|z—b—ia], (59 
w(z) ve F (z—b ia)log(z—b—ia)+vF(b—ia)log(z—b—ia), (60 
2 
where v= ; (61 
m(x+1) 
Using equation (3) and setting y = 0 we obtain 
2va2X 2vaY (x—l ¥ : 
Suu? — sue nee li —vKX log| (a—b)?- a*| +-yX (1—k), 


(cx—b)*+a? (a—b)?+a? 
(62) 
Suv 


> 


2va2Y 2vaX (x— b) 
(x b)? += (x- b)? | gq? 


—veY log| (~— b)?+-a?]- vY(1—x). (63) 
We require to annul these displacements by the introduction of potentials 
which have no singularities in the upper half-plane. From the results given 
in (ili)—(viii) and obvious rigid body displacement terms we obtain for the 
required additional potentials 


I37K 
Q(z) vF log(z—b-+ia) + nl , (64) 
K(z—b+7a) 
9 it 7 ail oq 
w(z) 2Q)(z) + anal +v(F—x«F)(z—b+ia)|log(z—b-+ia). (65) 
* 


It is clear from the above solution that we may solve problems which 
involve terms in the boundary stresses or displacements which are not 
expressible as Fourier integrals. We have only to introduce potentials 
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which give these terms at infinity. The latter potentials may give any 
additional terms which are expressible as Fourier integrals so that the choice 
of suitable potentials should not be difficult. The remaining terms may be 
dealt with by the general method given above. 

The above treatment of the elastic half-plane has the additional advan- 
tage that it leads on to a comprehensive treatment of the infinite elastic 
strip in terms of complex potentials. 

A brief treatment in terms of complex potentials of the half-plane with 
specified stresses over the straight boundary has been given by Milne- 
Thomson (12). In terms of the notation used above Milne-Thomson’s state- 
ments imply that for any given function 7° = f(x), the corresponding 
function J(z) is obvious. In the only application given he chooses f(x) to 
be constant so that the need for determining J(z) does not arise. The ques- 
tion of uniqueness of /(z) is not considered. Further, he introduces a single 
function of z to take care of conditions at infinity without showing that this 
s, in general, sufficient to give specified stresses at infinity. 
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SOLUTION OF TWO-DIMENSIONAL ELASTIC PROBLEMS 
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SUMMARY 

In this paper we give a general method of solution of problems of generalized plane 
stress of plate-like material bounded by a single unclosed contour. The theory is 
restricted to material of such a form that conformal mapping on to a half-plane is 
possible. That part of the boundary which is mapped on to the straight boundary 
of the half-plane is subjected to specified stresses or displacements, equilibrium being 
maintained by evanescent stresses at infinity. The method is illustrated by considera- 
tion of material with a parabolic boundary subjected to simple boundary loading. 


Fundamental equations 
THE notation is that of Stevenson (1, 2). We use complex variables 
z= #+1y, € = €+in related by the equation z = z(), where z'(¢) # 0 in 











Fic. 1. 
7 > 0; accents denote derivatives. Elastic material is chosen to occupy the 
region 7 > 0 of the z-plane and the theory follows closely that given by the 
present writer for treatment of the corresponding problems of the elastic 
half-plane. In terms of complex potentials Q(2), w(Z) the stresses are given 
by the equations 
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nw~ 
~~ 


20 = 26+77H) = = + 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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2 2(€& — yn — 21En) 


39 that 


~ (C Q)'(¢ = 0)’ (Z MC = 
10-8) Mie-1-éea) al ¥ i (S) | 22) « (5) e w O| [20 
~ \ © AG y Zz 


Let the displacement D, = u)+iv, referred to rectangular Cartesian axes 
(x,y) and have components wu, v in the directions of € and 7 increasing 
respectively. We write D ut+iv. Then (1), (2) 





, Q'(f) a'(Z) 
8uD, Su(Uy+ tv) KQ(C)—2(C) — = ——>. (4 
af) 28) 
=F 1 
And D Ea . Dy, (5) 
2 (C) 
so that SuD = [2'(0)2"(f) fn Q(0)2(Q—2(0)Q' (QZ)—a'(Q)].. (6) 


Relations between the complex potentials suitable for obtaining 
any desired conditions of stress or displacement along the 
boundary 7» = 0 





ae Or an i 
Let \(Z) uC), > (f) = | A(C) dé. (7) 
:) Let w(C) | 2’(C)P(f) dZ (3) 
d y' (CZ) y 
so that ne ae —A(€) (9) 
IZ 2'(C) 
ind (3) becomes 
> Qf) ae) , #0Z)—20) d Q'S 
$(777 +7&n) «(6 } (S): = (S) a AS) (10) 
C) =(6) 2'(¢) df 2(C) 
. - Te : ()'(€) 
hus (17 &7 +-v& 2re : 
Thu 1/77 +7En),,-o 4(7 +-1€y) re (8) (11) 
From (11) we see that the combination gives zero shear over 7 0 and if 
we require — 0 f (€), (12) 
we must choose ()(Z) so that 
()'(€) ; 
re : Py (é). (13) 
(c) 
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and (3) becomes 


~ O'(7 (CL z 3(7)—2(¢ Q’(z 
4554+-iFs) = 20) _* Olea 1 + 1)— als) d O'(C) (16) 
2 (¢) (CV L 26 z(¢) = df 2'(Z) 
~ Oy 
Thus 4(4n+7én)° = —2iim a (S) (17) 
g) 
From (17) we see that the combination gives zero normal pressure over 
n = 0 and if we require rae = 46), 
we must choose 2(Z) so that 
a . 
im — — 24(€). (18) 
Pe 
(iii) Let w(l) = — [ [eQ(0)2(Q+2(YQ'(Q] dé. (19) 


Then (6) becomes 

8uD = [2'(£)2'(2)] — KQ(C)2’ (2) + f2(Z)—2(C)jQ'(Z)]. (20) 
Thus 8uD° = 2relz’()z’(E)]-4[ KQ(E)2’(E)]. (21) 
From (21) we see that the combination gives zero normal displacements 
over 7 = 0 and if we require 


u® = #(€), (22) 

we must choose Q(Z) so that 
re[ 2’ (€)2’(€) |-*[«Q(E)2’(E)] = 4b(E), (23) 
i.e. re[ KQ(€)Z oa = [2’(€)2’(E) |*4pp(€). (24) 
(iv) Let w(Z) | [«Q(f)2’(f) —2(2)Q’(f)] a. (25) 


Then (6) becomes 
8uD = [2'(f)2'(Z)] [02 0)2(f) —«Q(Z)z’ (2) + f2(FJ—2(QHQ(D]. (26) 
Thus 8uD° = 21 im|[2’(E)2’(E)}-3[ «Q(E)2’(E)]. (27) 


From (27) we see that the combination gives zero tangential displacements 
over 7 = 0 and if we require 


v® = x(§), (28) 

we must choose 22(Z) so that 
im[2’(E)2’(E)}-*[«Q(E)2(E)] = 4yx(€), (29) 
i.e. im| «Q(E)2’(E)] = [2(€)2’(E) |! 4uy(E). (30) 


From equations (13), (18), (24), (30) we see that problems involving elastic 
material occupying the region 7 > 0 of the z-plane and subjected to specified 


stresses or displacements over n — 0 reduce to the determination of func- 
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tions of ¢ which are analytic in 7 > 0, of suitable orders of magnitude at 
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infinity, and which have specified real or imaginary parts when » = 0. As 
in the case of the half-plane we may show that if 

F(Z) e's"f_.(u) du, where f,(u) = f (tje-™ dt, (31) 
then F(f) is analytic in 7 > 0 and re F(€) = f(€). Also imi F(€) = f(E). 


Thus, we may satisfy (13) by 


x 


Q2'(¢ or as. 
; I(¢) = | ets“ f.(u) du; (32) 
(¢ 
0 
18) is satisfied by 
2'(¢) 25 fa 
EI J(¢) aia | e'S"h(u) du. (33) 
z(C) wT | 
0 


We may satisfy (24), (30) by equations of similar forms. 


Orders of magnitude at infinity 


If | tf (t)| dt < a, [ \td(t)| dt < x, 


then 1(¢), J(f) = O(¢-*) at infinity in yn > 0. (34) 
Thus, if z(Z) O(C*) at infinity (s ~ 1), it follows thatQ(Z) = O(f*-1) = o(z) 
and w(f) = O(¢*8-!) = o(z?) so that the stresses at infinity are evanescent. 


Ifs = 1, then Q(Z) = O(log Z) and w(f) = O(log Z) and the above conclu- 
sion still holds. 


Application of the above theory to material with a parabolic 
boundary 


Consider the transformation 


({-+ia)? = 2iz. (35) 
The boundary 7 = 0 is given by the Cartesian equation 
a? 2a*y as (36) 
and the locus é -b is given by 
x2— 2b7y 54, (37) 


The curves given in (36), (37) intersect at the points (ab, (a?—b?)/2) 
denoted by A and B in Fig. 2. Also 


2'(C) = a—2¢. (38) 
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Thus z'(f) is zero at € = —ia but does not vanish in the region ” 
occupied by the elastic material. 
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We now apply the above theory to the determination of stresses and dis- 


placements in elastic material which occupies the region 7 > 0 and is sub- 
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19 os ae (—F§", if b, ~ ; 
(ii) Given nn” = fo(€) €x° = 0 (46) 
\ 0, \€| > 6, 
. 27 i om C—b = 
From (32) I,(f) 2bf+ Clog -—. (47) 
(C) i C+b 
/ ee oe eS ee 
nd 03(2) “| 2abe 2ibl?+-al? log > —1@ log =], 
sisal rr : C+5 C+b 
whence 
672,(C) 21) tabl? + 31bC(b? "hap | + | 3(¢4—b4) + 41a(f3 —b’) llog(f—b) — 
tq [3(¢4—b*) + 40a(C3+-b) |log(f+-b). (48) 
From (7) 
mA,(C) 460? + Siabl+ 2a7b— 268 + 
b4-2iab? h4—2iab3 ri ” rh 
ones I — | 2¢?— 41al?— 2a? ]log in (49) 
c= b ¢ b C+5b 
nd 
ae - et eee eo ee 
671(C) 6603 + 16iabl? — 6b3C + 3C4log -—— — 8ia@log =—— — 
oar ; : i C—b ¢—0 
} (+b — wil ° 1 
6a2C2 log = 4. (3b44 6a2b?)log = 4+. 4rab? log(l?—b?), (50) 
. ¢ b —- 
whence 
nd dis ; 
is sub- ats G) 
3| — 31b05— 6abl4— (3a? + 5b?)ibl3 + 6ab(3b?— 2a?) l? + 31b3(4b? + 3a?) 0] +- 
L[i(Z6—b8) 4 2a( 05-5) + ia2( 244) + 403( £3 +63) 4 
(39 | f3(37h 4a)(C? b?)— 2ab3(3b+ dia)(C - b) |log(f+ b)— 
1a(C' b6) 2a(C°—b°)- ia>(C4 b+) 4a3(C3 — 53) - 
(40) b3(3ib+ 4a)(C2?—b*) — 2ab3(3b—4ia)(f—b) log(f—b). (51) 
toe ; ( €2—b2, |é| <b, * ; 
(4] il) Given n° €&n° = 0. (52) 
ot - b, 
This is the problem of hydrostatic loading up to the level A B in Fig. 1, the 
solution being interpreted in terms of plane strain. The complex potentials 
(49 ire given by a linear combination of the two preceding pairs of potentials, 
- (2.(C) 67Q,(C) —Q,(C) w(C) b*w,(C)—w(C). (53) 
(43 ; 
From (1), (2) 
S7 r. qe ™ 
(44 20) — (& 7? —b?)(0,—0,) 4 “(b+ €log 3) (see Fig. 1), (54) 
— 8y €*@+an+a*+71én | r : 
20 | | 2b+ & log —— (0, —8,) | + 
7 (a-+-7)*-+6* | Ny v 
n | - 
+7] log ++ &(6,—6,) \ (55) 
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whence 


nén[(a+n)?+€] 


= 2n| 26+ n(2€2-+-a?+an)log 1 + £(0,—0,)(22—72+a24an)|, (56 


lg 
and 
mipal (a+ n)P+€] = 2by(y?—€—a*) + 2Eq2(a-+ 2y)log 1 + 
I 
+ (0, —6,)[ (€2—b2)(€2+-a2+ 2ay+ 42) — 4?(q2-+-02—302)], (53) 
and &é ©Q— 77. (58) 


The stresses are finite and uniform in n > 0 and are evanescent at infinity, 
The displacements may be evaluated from (4). We write 

(¢) 
( ( 


From (42), (48) we see that D® is uniform and continuous in » > 0. From 
(40), (44), (47), (50) 
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(60) 
Thus Df is uniform and continuous in 7 > 0. Equations similar in form 
to (56)-(58) were given by Conrad (3). His solution to the problem is, how- 
ever, not valid as the integrals chosen for complex potentials are divergent 
and the mathematical operations on such integrals are not permissible. The 
integrands may be modified as done by Filon (4) when dealing with the 
infinite strip. There is, however, no need to work with potentials in the form 
of integrals as the problem requires only elementary functions. We have 
contented ourselves with obtaining solutions in which the stresses at infinity 
are evanescent without specifying orders of magnitude at infinity. We now 
show that the solution is not unique. From Stevenson (1) the condition that 
a boundary is stress-free is that 
Q(z) +20'(z)+w'(z) = constant (61) 
along the boundary. Thus, taking the constant to be zero, in the present 
case the boundary 7 = 0 is stress-free if 


20(f)(€-+7a)—Q’ (C)(C—ia)*? + 2iw'(f) = 0. (62) 
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Equation (62) is satisfied by the potentials 


us 


Q,(f) =f,  wo(L) = (22+ 152). (63) 

6 
These potentials give evanescent stresses at infinity and may be super- 
imposed on the above solutions without violating the boundary conditions. 
If we put Q(C) ¢ in (8) the term in w(Z) of highest degree is iZ3/6. Hence 
63) may be employed to remove the term in Q(€) of order Z at infinity and 
that in w(C) of order 2° from all solutions obtained from (8). The stresses 
given by (63) are 


2¢€[(n +a)? er (n+a)(2€?+ 2a?+ 2an+ 7”), (64) 


poms 


277| (yn +a)*+ &}? n(n+a)(n+ 2a), (65) 


2€n[(n+a)?+é]? = én(y+2a). (66) 
In all the solutions obtained above the stresses at infinity are O({-*), i.e. 
()(2-') (see equation (34)). We may use the potentials given in (63) to obtain 
solutions to the problems defined in (39), (46), (52) with stresses at infinity 
which are O(z-!). For large |f| (42), (45) become 

Q,(¢) 450 /7+ O(log €), w,(C) —2ibf3/32+ O(log lf). (67) 


The required modified potentials Q, ,,,(¢), 4,,() are 


““1,m\'>/> 


4} 4b 


' : ) 7 ‘ es y _ 
G2 mS) (2,(¢) ge {2,(¢), 4 (4) w,(C) + . wo(C). (68) 
Similarly, for large |¢| (48), (51) become 
Q,(Z) 1h30/307+O(logf),  — w(£) = —2ib323/97+ O(log Z). (69) 
The required modified potentials in this case are 
" ’ 4h3 . ae pa 4h3 y - 
(ds in(S) (2,(¢) + ——22,(4), We m(S) Wy(C) + ; Wo(C). (70) 
OTT OTT 
Thus the modified form of (53) is 
) ; 8b3 0 ; Rh 7 _ 
G23 m(S) 023(¢)-4 = £25(); W3 »»(S) w(C) Tr 3m Wo(C). (71) 
OTT v7 


The modifications are equivalent to superimposing multiples of the stresses 
given in (64)-(66) on to the stress systems given by the general method of 
solution. The additional stresses have a profound effect on the stresses in 
the finite part of the material even though they do not affect the boundary 
conditions over 7 0 and give evanescent stresses at infinity. The solution 
emphasizes the importance of specifying stress components at infinity other 
than those which are O(z-1). In the above example the orders of magnitude 
of the modified potentials are the lowest possible since the stresses applied 
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along 7 = 0 have a non-zero force resultant. Their uniqueness (5) may be 


shown by a method similar to that given for the half-plane (6). 
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1. An integral occurring in Probability Theory is 


} oe t(t sin 6)"-? i eee a 
Z,,(z, a) : -exp{— }(t?+-a*— 2at cos 6)} dédt. 
J J [(4n—$)2!"- v7 . 
0 0 is x 
This corresponds to the non-central Chi-square Distribution (a + 0) and 
represents the probability (1) that an observation (2,, 29,..., 2,,) falls within 
, given distance oz of any point P(a,, dy...., a,), the origin being at the 


point O. where 


OF oa (> a?)} 


\ uys* 

An interesting application of this in Ballistics is to the problem of esti- 
mating the relative risks at different distances from an aiming point when 
each projectile has a zone of risk. 

In any specific problem, the value of Z,,(z,a) for small values of z and a 
will be required. The object of this note is to give a few relevant formulae. 

[am grateful to Dr. D. 8. Kothari, Scientific Adviser to the Ministry of 
Defence, Government of India, for encouragement and for permission to 
publish this note. 

2. Since ears? _ ST H(z) (2.1) 

er r. 


where H,(x) is the Hermite polynomial defined by 


r 


=A J - Dy 
Aa) = S$ SCPE gots, (2.2) 
aa 2)(1,9) 
The function exp{— }(t?+-a2— 2at cos 6)} can be expanded in powers of ¢ and 


uso in powers of a. 
If we consider the first alternative we get 


r 


, t(tsin 6)"-2e-!@ wot 
Z,(z,a) | | iss . S H(acos 6) dédt. 
T(dn—14)2!"'-17 Gr! 7 
- = r=0 


+ My thanks are due to the referee for advising me to rewrite the introductory section 1 
f this paper. 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 














362 ON THE PROBABILITY FUNCTION 


Term-by-term integration with respect to t is obviously permissible, and 


hence 
Z,,(z,a) = —— el aon a - | sin”-2@ H.(a cos @) dé 
‘ D(4n—4)2*"°-Wa Gy r! (r+n) | ”" 
= r=0 0 
aw eS et ns 
Qin-la/a &, (2r)!(2r +n) Z 2/(1, 9) (r—j+4n) 
F ) ' : - 
since . 


| sin’”’6 cos”6 d@ 0 when n is odd 


0 


> 


C(4m+4)0(4n-+-4 ' 

: a)? (an +3) otherwise. 
[(4m+4n+1) 

3. The other alternative, contemplated in section 2, gives 


, " f t(tsin 8)" Gar 
0 0 il “ r=0 * 
Since Z,(z, 0) 3 : t2-lp—3l dt 
b 4n\~> - Din 1T'(4n) ; 
alia’ 
we get the expansion 
x j= . 
, Se i(—2r,2 — : : 
Z,(z,a)= = S a oo D or 21(r—j+43)Zop_25+n(2; 9). 
mony ual j=0 (159) 


4. An expansion similar to the one in section 3 is obtained by expanding 
exp(at cos @) and integrating term by term. We then have 
F we oF . 
Z,,(z,a@) = exp(—4a?) 5 : — Zn+or(2; 9), 
5 y 2°] (r- 1) ae 
r=0 


where we have used the duplication formula for the gamma function. 


5. An interesting reciprocal formula will now be obtained for the function 
Z,,(z,0), the probability corresponding to the central Chi-square Distribu- 


tion. Taking |h| to be small, we have 


a Jin-3T/1 1 ° = gb 
S 0Z, (2,0) — = ets) e-th dg — YS _* H,,(h) 
a j &, Pn+2)" 


true for small values of z. 
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SUMMARY 
Conditions are given which are sufficient to ensure that the current density normal 
an edge is zero at the edge and that there is no line distribution of charge on the 
lige. An extra condition is given which makes the components of the field parallel 


the edge finite. The solution is then shown to be unique. Simpler conditions are 
given for two-dimensional fields. The agreement of various known solutions with 
ir assumptions is discussed. Finally it is shown that certain simple current densities 


ead to a unique solution for the diffraction by corners of flat surfaces. 


Introduction 

BouwKAMP (1) has indicated how one may obtain any number of solutions 
tothe problem of the two-dimensional diffraction of plane harmonic electro- 
magnetic waves by a perfectly conducting semi-infinite plane by taking 
derivatives of any known solution. Two different solutions for the three- 
dimensional diffraction by a semi-infinite plane have been given by Wiegrefe 
2) and by the author (3). Various solutions have been given for the 
diffraction by a circular disk; two different ones are those of Méglich (4) 
and Meixner and Andrejewski (5). Asa consequence it has become necessary 
to impose conditions that ensure a unique solution. 

The question has been discussed by Bouwkamp (6), Maue (7), Meixner (8), 
ind the author (3). The results of all authors are in agreement although 
their viewpoints are different. In all cases the edge has been taken to be a 
curve with continuously turning tangent and apart from Maue’s work the 
edge has been on an (essentially plane) infinitely thin conductor. Bouwkamp 
has argued that the field and currents must behave near the edge as they do 
in the Sommerfeld two-dimensional solution for the semi-infinite plane. 
Meixner uses the fact that near the edge the fields are quasi-stationary and 
the assumption that the integral of the electromagnetic energy over any 
region must be finite to deduce the forms of the fields. It is then possible to 
deduce the currents and charges; they agree with those found by Maue using 
an integral equation formulation. The author attempted to show that the 
issumption of finite current and charge on the conductor gives a unique field 
and that the current perpendicular to the edge vanishes, but the proof is not 
completely rigorous. 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3, (1952)] 
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It follows from the work of Meixner and the author that, whether or not 
the bounding curve of an infinitely thin surface is plane, a sufficient condi- 
tion for a unique solution (apart from radiation conditions at infinity) is that 
the integral of the Poynting vector over a small cylinder surrounding the 
edge should tend to zero as the radius 6 of the cylinder tends to zero i.e, 


lim { (EA B*+E*,B).dS = 0, 
te) 


wie, 


where S, is the cylindrical surface and the asterisk denotes a complex 
conjugate. 

For the actual solution of problems it is necessary to reduce this condition 
to simpler forms of various kinds. If the problem is formulated in terms 
of the field it is sufficient that the components of E and B parallel to 
the edge should both be O(1) and that the remaining components 
should be 0(6-!). In two-dimensional diffraction when the fields depend 
only upon £, or B, the integration is round a small circle and sufficient 
conditions are that E,, = O(1) or B, = O(1) in the two cases, as has been 
indicated by the author (9). Frequently, however, the problem is formu- 
lated in terms of integral equations involving the current and charge densi- 
ties. In this case one wishes to know firstly whether line integrals along 
the edge are necessary and, secondly, conditions on the current and 
charge densities to ensure that the solution to the diffraction problem is 
unique. 

In this paper we are concerned with the integral equation approach and 
the analysis is kept fairly general; the curvature of the surface is not neg- 
lected except in section 9. The work is complicated by the fact that we have 
to allow the charge and current densities to be infinite at the edge since all 
known solutions have singularities at the edge. In section 1 there is a brief 
discussion of the validity of the integral expressions. Section 2 shows that 
the singularities (if any) of the field are derived from the intensities and 
potentials of certain static charge distributions. Section 3 is devoted to 
enumerating the conditions imposed on the current and charge densities and 
the reasons for their choice. As a consequence of these conditions we show 
in section 4 that the current density normal to the edge is zero. In section 
we give sufficient conditions for the field components parallel to the edge to 
be finite. It is shown in section 6 that as a consequence of these restrictions 
the field is unique. Section 7 deals with the modifications for two-dimensional 
fields and section 8 indicates how some known solutions agree with our 
restrictions. Finally, we discuss the uniqueness of certain current distribu- 
tions in the corners of flat surfaces. 
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|, The validity of the integral expressions for the field 

Weshall use the definitions of regular surfaces and regular curves that are 
siven by Kellogg (see 10, ch. 4) with some restrictions. Thus if z = f(x,y), 
where (x, y) lies in a regular region of the (x, y)-plane, is the representation 
of the set of points forming a regular surface element we shall assume that 
‘is continuous and has continuous second partial derivatives bounded in 
ibsolute value by some constant M. The (x, y)-plane will usually be the 
tangent plane. If y = g,(2),z = g(a) is the representation of a regular arc, 
we shall assume that g, and g, are continuous and have continuous second 
lerivatives bounded in an absolute value by M. Apart from this, in this 
section we shall assume that the regular surface has a continuously turning 
normal and that the regular curve has a continuously turning tangent except 
when specifically stated otherwise. 

Let S be an unclosed perfectly conducting two-sided regular surface sur- 
rounded by free space. Let the field whose electric intensity is Ey) and whose 
magnetic flux density is B, be incident upon S. We use m.k.s. units and 
suppress the time factor e'*’. The incident field excites a surface current 
density I and a surface charge density o on S which satisfy the continuity 
equation div I+iwo = 0. There may also be present a line charge of linear 
density o, on C, the bounding curve of S. Then a solution of Maxwell’s 


equations (under certain conditions) is 


° . | ” 
~ . tw ° ° . tw 
E, = E, ie {k2I +-iwo grad wb} dS — igh ‘) o, grad ys ds 
dork Th 
; > 


B, — B, + : | (TAgrady} ds 


where E,, B, are the electric intensity and magnetic flux density at the point 
of observation (x,, ¥,,2,), k 27/(wave-length), 4 = e-*"/r and 
y2 (x, —2)?+(y, y)?4 zy z)2, 
x,y,z) being a point of S or C. The current density I must lie in S and 
therefore satisfies I.n 0. where n is a unit vector along the normal to S, 
being chosen so that it is always directed from the same side of S. 

In order that the integral expressions may have a meaning when the point 
of observation is in free space we assume that I and o are summable 
Lebesgue) on S. Neither of these assumptions implies the other. The assump- 
tion implies that I and o are finite almost everywhere on S. We shall assume 
that the only place where finiteness and continuity may fail is on C. If 
part of the surface goes to infinity a difficulty arises because a summable 
quantity is also absolutely summable and this may not be true for the 


current and charge densities occurring in practice. The difficulty can usually 
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be overcome by taking the surrounding medium to be slightly conducting 

until after the integrals have been evaluated. Another device which is some. 
‘ Ww 

times successful is to define ( as lim {. For simplicity we shall assume in 
Segoe 

the following that S is finite in area and C is finite in length. Also, for the 

line integral in (1) to have a meaning, we require o, to be summable on ( 

and therefore finite at almost all points of C. We assume that the finiteness 

of o, can fail only at points of C where the direction of C changes dis. 

continuously. 

The method of verifying that E,, B, satisfy Maxwell’s equations at points 
of free space is well known (see, for example, Stratton (11)). The verification 
requires firstly that it shall be permissible to take derivatives under the 
integral sign, secondly that Stokes’s theorem can be applied to > o,¢ di 

e C 
and > o,grad¢ ds, and thirdly that o, = nAI.t on C, where t is a unit 
é 
vector tangent to C, its direction being determined from n by the right-hand 
screw law. 

There is no difficulty about the taking of derivatives since the integrand 
is the sum of products of a summable function and a continuous function 
when the point of observation is not on S or C so that the integrals exist. 

As regards the third step, the remark already made about o, implies that 
n ( I.tis finite at all points of C where the tangent is continuously turning, 
i.e. that the current normal to C is finite. 

Stokes’s theorem for a two-sided regular surface can be proved whenever 
Green’s formula in two dimensions holds for a regular region of the plane 
which is bounded by the projection of the boundary of the surface, i.e. when 


| & ~. | dandy = | (M dy+N dr), 
J \cea oy , 

i J 
where J is a closed regular curve and J its interior. 

In most proofs (see, for example, Verblunsky (12)) it is required that J/ 
and N shall be continuous on J+J. This restriction is too heavy for our 
purposes since we do not wish to assume that all components of the current 
density are finite on C, although we may take them to be finite and con- 
tinuous on S. We can, however, easily extend the theorem to cover our case. 
If 2M /éx and éN /éy are summable, if M is the indefinite integral of ¢M cx 
with respect to a, if V is the indefinite integral of @N /éy, and M and N are 
continuous on J, then Green’s formula holds for J,,, where J, is a closed 
regular curve in J. In particular, this is true when J,, is a member of the 


m 
sequence {/J,,} of simple closed regular curves x = 2,,(t), ¥ = y,,(t) on I J 
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(t) = y(t) uniformly for 0 < t < 1, where 


i 
m— « m 


m 


ach that lim x,,(¢) = x(t), lim y 


a(t), y = y(t) (0 <t <1) is J. Hestenes (13) has shown that this 
N,, be the 
,) Where 


sequence exists when J is a closed rectifiable curve. Let M,, 
values of M and N on J,, and define My’+ Na’ as lim(M, 


, Ty 
m¥mT N m* 


m 
m= 


primes indicate derivatives with respect to t. Now the measure of J—T,, 
tends to zero as m —> oo and hence 


, * (eh N 
lim | a : dady = 


Ox c Yy 





j (= aS on ) dady. 


Ox oy 


m—> =x 


ad 
Thus Green’s formula with M dy+N dz replaced by (My’+ Nz’) dt holds 
whenever we can assert that 


m—> & 


lim [ ( M., Ym T Nin . ) dt wees [ (My’ i Nz’) dt. 
J J 


, 


< K for all m 
ind almost all ¢, where K is independent of m and t, or (ii) {1 Yin +N 2m} 
monotone and converging to My’ + N2’ with both My’+ Na’ and M,y,+™, 2; 


Some conditions which are sufficient are (i) |M,, y,,+NnXm 


summable, or (iii) {J,, y’,,+-N,, 2} converges uniformly to My’+ Nz’ at 
almost all points of C. In most physical situations the assumption (ii) will 
be satisfied; it requires that the current density perpendicular to the edge 
is finite almost everywhere on the edge and monotone in the neighbourhood, 
and we have already assumed the first part of this in connexion with the 
thirdstep. In the particular case when J has a continuously turning tangent 
we may take J,, to be a(t)—m-1y’(t), y(t)+-m-1a'(t). For large enough m the 
equations 2,,(f) = 2,,(t). Y(t) = Ym(to) will have no solution other than 
t,. Then, if JZ, o(m), N o(m) as m > , 


m 


My'+ Na’ lim (JZ 


m 


mY +N, 2’) 

so that if My’+ Na’ is finite on J, (ii) is satisfied when the sequence 
smonotone. We shall assume that one of the conditions given can always 
he complied with. 

Accepting (1) as the expression for the field we obtain three integral 
equations for the determination of I and o from the boundary conditions 
that the tangential components of E, and the normal component of B, 
vanish on S. Only two of these integral equations are independent. 


2. The reduction of the fields due to the surface distribution 
Write (1) as E, = E,+ E+E’, B, = B,+B where 


sy (ta9/4crk?) { 07 grad ub ds. 
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Now b—r-| = |r—(e-*r—1)| = O(1), 
grad —grad r-! ikr-l(e-tr_ ]1)+-r-2(e-tr 4 jkr—1)||\gradr| = O(1 
and hence 
ia i > ° Y 
E= — sa | {k?r—1I +-iwo grad r-!! dS+ O(1) 
4nk? | 
S 9 
a 
B= | TA grad r—} dS+O(1) 
In particular 
E iw [ f[2p-1J 15 2 dS+oi 
“Utah? | eer *L, + twa(y—yy)r} dS+-O(1), 
§ 


where the suffixes x, y, z denote components parallel to the coordinate axes 
Thus the part of Z, which is possibly non-finite is the sum of the potentia 
due to a surface charge distribution whose density is proportional to J, and 
the y-component of the intensity of a surface charge density proportiona 
to o. In fact it is clear from (2) that E and B can be expressed in terms 
of the intensities and potentials of certain charge distributions. Therefore 
if o and I satisfy appropriate conditions, the theorems obtained by the 
author (14) can be applied. The conditions and the reasons for their choice 
are given in the next section. 


3. The conditions imposed on the charge and current densities 
We restate conditions B given in (14). They are that 
(i) ois summable on S. 

(ii) |o—o,| < A(e)r{, where o, is the density at P of S, o the density at 
a point distant 7, from P, « the shortest distance of P from C, and 
A(e) and a are constants and 0 < « < 1. The condition is satisfied 
uniformly on S, except near C’, where the condition holds in 7, < a 
0 <a < 1,and A(e) is constant in 7, < ae but such that 


A(e) = o(e—!-*) ase > 0. 


(ili) |o o(e!)ase> 0. 

(iv) The line integral of o along any line on S, whose points in comm 

with C (if any) form a set of (linear) measure zero, is finite. 

It will be assumed that I also satisfies the same conditions with o, A.» 
replaced by I, A’, «’, where 0 < a’ < land A’ = 0(e-!-*). These conditions 
will be known as B’. Since I satisfies B’ the separate components of I als 
satisfy B’. 
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The reason for the choice of the modified Hélder continuity conditions 
B(ii) and B’ (ii) is that the continuity alone of the charge density is in- 
sufficient to ensure the existence of the tangential components of the 
ntensity on S. This is shown by a counter-example given by Kellogg (10). 
In practice, of course, I and o will be differentiable on S and we can take 

L, a’ 1. The modification of allowing A and A’ to become infinite 
s ( is approached is introduced in order to allow I and o to become infinite 
n(. Ifthis were not permitted all known solutions to diffraction problems 
would be excluded. We may note, however, that the component of I normal 
‘ since we 
have already assumed that component to be finite on C. The use of Stokes’s 


to ( can be regarded as satisfying a Hélder condition with finite A 


theorem in section 1 is also justified. 

Conditions B (iv) and B’ (iv) result from a consideration of the case when 
Sis plane and C is a straight line. There we can expect variations parallel 
to('to be of no importance and that the charge per unit length parallel to C 
is finite, i.e. the integration of o perpendicular to C gives a finite result. 

Conditions 6 (iii) and #’ (iii) are virtually implied by conditions B (iv) 
ind B’ (iv) in any physical problem. Conditions B (i) and B’ (i) express the 
fact that the total current and charge on S are finite. Thus the conditions 
mposed, although appearing somewhat artificial at first sight, do in fact 
havea physical basis. It may be assumed in the actual solution of problems 
that conditions B (iv) and B’ (iv) will be dominating ones. In a later section 
we see how some known solutions agree with our conditions. 

It may be noted that it is not possible to deduce immediately restrictions 
mn the field from those imposed on the current and charge densities since 
these densities are related only to the discontinuities of certain field compo- 


nents across WN. 


4. The consequences of the conditions imposed 

Let O be a point of C in a neighbourhood of which C has a continuously 
tuning tangent. Let Oryz be asystem of axes with Oxy as the tangent plane 
toS at O and Ox as the tangent to C at O. Let PXYZ be a system of axes 
with PXY as the tangent plane to S at P and such that PXYZ becomes 
Yryz when P coincides with O. Then, since o and the components of I 
satisfy conditions B and B’, Theorem 1 of (14) shows that, at P, 

E, o(e7! log e'). E o(e-1+ |loge'), ‘!B 0 (e—!) 
SéE> VU, W here € PO. Therefore 
E, o(e-*). (3) 

Also since the component of I normal to C satisfies a Hélder condition 

with finite A’, so does oj, i.e. , satisfies conditions B,, namely 


i) o, is summable on C and 


5099.19 
<1 


Bb 
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(ii) |o,—o49| < A,r}, where A, and y are constants, 0 < y < 1, oy is the 
density at O and 7, is the distance between O and the point of ( 
where og, is the density. 

One of our boundary conditions is that E,;-+ £y+ Ey} = 0onS. Since the 
relation for Ey in (3) is derived as if o were present but not I, and since ¢ 
and o, satisfy B and B, respectively, we can apply Theorem 4 of (14), and 
thus o, = 0. Moreover it is clear that Theorems 2a, 4a, and 46 also apply, 
Thus we have proved that if I and o satisfy B’ and B then there is no line 
distribution of charge on C and no isolated charge at a vertex of C. This 
may be restated thus: if I and o satisfy B’ and B respectively, if S hasa 
continuously turning normal, and if the component of I normal to C, where 
C has a continuously turning tangent, is finite, then that component is zero, 
Also, if C is the continuously turning edge on a surface where the direction 
of the normal changes discontinuously, then the component of I normal to 
C is continuous across C. 

It also follows from Theorems 2 and 3 that at Q, where Q is not necessarily 


on 8, E,| = 0(e~), B,| = 0(e-) (4 


as « > 0 where e is the shortest distance of Q@ from C. 


5. The behaviour of E;, and B,, 


Copson (15) has conjectured that the component of E, along C vanishes 
at C. This is known to be true in a number of solutions, but solutions are 
also known in which it is not true, e.g. Wiegrefe (2). We shall give sufficient 
conditions for the component to be finite and also give a counter-example 
to show that the conditions already employed are insufficient. It will be 
assumed that the results of the preceding section hold. 

It is convenient to make a change in notation so that the (2, y,z) coordi- 
nates become (2), 2, x,) and the point of observation Q becomes (X,, X9, X3). 
The equation of Sisx, = f(x,,x,) and the parametric form of Cisx, = G;,(t), 
= G,(t), x, = f (x, x2) where G,(0) = 0, G,(0) = 0, G{(0) = 1, GA(0) = 9, 
the primes indicating derivatives with respect to the argument. Let &,, &, £5 
be curvilinear coordinates so that S is &, = 0, €, > 0, C is €, = 0, & = 9, 
and O is £, = 0, & = 0, &, = 0. In particular we may take 2, = G,(§,), 


t, = G2(E1) +o, %3 = Egtf (21,2; 


> 
“ 


OX}. CX}. 
h;; = h;; = >; k : La 
k 0&; 0&; 


). Let h,;; be defined by 


Let Q be (,, 2,53) in the curvilinear coordinates and H,; be the value of 


h;; at Q. 
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Suppose that o(f,,&,) satisfies conditions B, of (14), ie. in |€&| < &, 
t, < &, then the following expansion holds: 


o(€,,&) = o(€)+&, 0’ (&) + €o"(E,, &5) 


ghere o” is summable on S and o(€,), o’() are both O[£ ‘{log(1/£_)}-?-F] 
where 8 > 0. (It may be noted that we can replace €, by e.) Then we may 
pply Theorem 5 of (14) to that part of E,, which is due to o and hence that 
part of E,, is finite. There remains the part due to I to be considered. Apart 
from multiplicative constants 


E; > g Om, rl. (hy, hog—hi,) dé, dé,+ O0(1), (5) 
. : 


where 7? = > (X,—2;,)?. 
k 


Let I... I:, be the components of the current density along the curves 
constant, £, = constant respectively. Then J, Jz,, the components 
long and normal to €, constant, are given by 


) 


, ‘ x 
le, (1/shy,) S "7 = Te, +{hye/y (hy hoe) }le,, 
—_ S13 


I’, = a/(1—h2q/hy, hog) Ie,. 


> 
< $2 


1 
Since we have proved in section 4 that Jz, > 0 as €, > 0 and since 
y (1—hijg/Ryy hoe) F O 
it follows that J; > 0 as €,> 0. Now 
, % OX, Ga, @X,\| 
a fe > tal t_*n\| < MRI 
é ye alee, &]| 


D 
where R? = > (=,—§&,)?. 
k 


Since r? > K, R® (from equation (21) of (14)) and I is summable on S we 
may repls ace > [_., 0X;,/EZ, in (5) by Ig, vhy,. Thus the only part of Ez, which 


may be non- finite i is the potential of a charge distribution Jz. Hence if J, 
satisfies B, we may apply Theorem 5 of (14) and hence this part is O(1). 
Thus we have proved that if o satisfies B and By, and if I satisfies B’ and 
I; satisfies B, then Ez, = O(1) in the neighbourhood of C. 

Let us now consider B;, which is given by 


& & 
l f . eX, . 
B; ay iy (X.—z.)\ 
$1 4nvH,, | | | . , A six ez, I,,,( k vx) 


(hy, ho—h2,) dé, dé,+O(1), (6) 
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where A;;, = 0 if i, 7, k are not all different and otherwise 4,,, = +| 


according as 7, j, k are in cyclic or anti-cyclic order. Now 


ex, ex, eX,eX; eX; eX;\| , 
— I.—— I, - A tant —=* — —* "|| < MRYT | + | 7, }. 
OB, * OB, ~° WViheg “*\05, 08, 05, 05,)| ou ” 


Therefore, since J;, satisfies B, and J;, is a continuous bounded function 
, , Poe ] ox, 
we may replace J,, in the integrand of (6) by —— I;,——? and the equation 


Nico ** OB 


still holds. Further we may replace X,;.—., by 


since X;,.—2,.— > = (€;—=,)| < MR?. 


Then, since r? > K, R? and &, = 0, 5, = 0 


er ge 
$1 So 
. 


_ 2 


Be, < Gv. | | (\Te,|/ Vag) R-8y) (hy higg—h3,) dé, dé, + O(1) 
g 0 


where J = 6(X,, X,, X3)/é(E,, ,, 3) and R? = + (F,—é,)?+53. Now], 
is bounded and zero at £, = 0. Therefore by defining J;, as zero off S we see 
that the expression on the right is, apart from bounded multiplicative 
factors, the normal derivative of the potential of a continuous charge 
distribution which is known to have a finite limit when =, > 0. Hence 
B; O(1) in the neighbourhood of C. 

[f C instead of being the boundary of S is the continuously turning edge 
where the direction of the normal to S changes discontinuously, we may 
apply the above analysis in turn to the portions separated by C. The only 
difference is that J;, is non-zero but continuous across C’. Hence in this case 


E. 


i 


O(1), Bz, = O(1). Clearly the same is true of the appropriate com- 
ponents of E, and B,, since o, = 0. 

As a counter-example we may consider the case when S is plane and Ca 
straight line so that €,, €,, €, are Cartesian coordinates. Take 


2)}]. 


I, = 0, Ie, iw€,/{log(1/&5)}, o 


wr 


1 [E,{log( I 


Uv 


Then the only non-finite contribution to Z,, is from o and it follows from 
the counter-example in section 8 of (14) that this becomes infinite like 
loglog(1/=,). Thus conditions B and B’ are insufficient to ensure the finite- 
ness of £;.. 
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6, Uniqueness for the edge with continuously turning tangent 

We have now obtained sufficient conditions for the solution (1) to be 
unique. Let S, be a small cylinder of radius 6 surrounding the edge. On S,, 
B, \Bt.n = Eye, BY,— Biz, Ey,, where E,, and Bj, are linear functions with 
pounded coefficients of the components of E, and B* respectively. From 
sections 4 and 5 it follows that, if o satisfies B and B, and if I satisfies B’ 
ind I, satisfies By, then Ey¢, BY, — By¢, Ey, = 0(8-1). Since the radius of 8, 


sswehave | E, \ B*.dS = o0(1)as8 + Oand hence the field is determined 


Si 

miquely by (1), assuming that the appropriate radiation conditions at 
nfinity are satisfied. We are referring here to the case when S and C are 
finite; when they are infinite restrictions on the behaviour of the charge and 
wrent densities at infinity have to be introduced to ensure a unique solu- 
tion. In view of the last paragraph but one of the preceding section it is 
mmaterial whether C is an edge or the boundary of S. 

One consequence of this result is that there is only one solution differen- 
tiable on S such that Co O(e™ 7/5 O(e*:), Ie, = O(1) where Oy Po —l, 
, > —1, since such a solution c 8 ar ly satisfies the conditions imposed. This 
is the type of solution in which one would probably be interested when 





ving a practical problem. 


7. Two-dimensional fields 

The behaviour of two-dimensional fields is similar to that in three dimen- 
sions, but since the conditions are simpler it seems worth while to state them 
separately. Let C, be a regular curve with end-points F, and F,, and consider 
the behaviour near F,. Suppose that the fields are independent of the z- 
oordinate when the notation of sections 1-4 is used. Then 


tw 


. ‘. 
E E, (—k*?.bI—iwo grad %) ds —-—.. o, grad 
; 4ork? | ) 4nk2 *° ¢ 
-_ (7) 
B, B, r (I grad us) ds 
where miH\(kr), 7? = (y,—y)?+(z,—z)*, HY is the Hankel function 


‘the second kind and in the o, term 7 is measured from F,. Formulae (7) 
satisfy Maxwell’s equations if (i) it is permissible to take derivatives under 
he integral sign, (ii) it is permissible to’integrate [ op ds, | o grad ds by 
Cy C1 

parts, o being the quantity integrated, and (iii) o, = I, at F. These con- 
ditions can obviously be satisfied. We can replace % in (7) by —2logr 
without making the expressions less finite. Then E, and B, are obtained 
irom the potentials and intensities of certain charge distributions. 
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Let conditions B (i)’, (ii)’, B’ (i)’, (ii)’ be the same as B (i), (ii), B’ (i), (jj 
except that S is replaced by C, and C by F,. Then, if o and the components 
of I satisfy B (i)’, (ii)’, (iii) and B’ (i)’, (ii)’, (iii) respectively, Theorem 6 of 


(14) shows that E, = 0(e-), B, = 0(e"). (8 


Then Theorem 7 of (14) shows that co, = 0 at /. Similar results hold at £ 
and at a point where the direction of C, changes discontinuously. Hence 
we have shown that if o and I satisfy B (i)’, (ii)’, (iii) and B’ (i)’, (ii)’, (iij 
and if the current density perpendicular to the a-axis is finite it is zero at an 
end-point of C, and continuous across a point where the direction of ( 
changes discontinuously. . 
The discussion corresponding to section 5 is considerably simpler. The 
possible non-finite part of E, is the potential of the distribution J,. If 
[, = Ofe{log(1/e)}-®-*] this part is finite and Z, = O(1). It is easy to show 
that B, is mainly the normal derivative of the potential of the current 
density component perpendicular to the x-axis and hence B, = O(1). It 
now follows by a use of (8) that, with suitable conditions at infinity, there 
is only one solution of the type (7) if o satisfies B (i)’, (ii)’, (iii), if I satisfies 
B’ (i)’, (ii)’, (iii), and if, in addition, J, satisfies the condition just given above, 


8. Some known solutions 

In this section we briefly indicate how various known solutions agree with 
our conditions. Before doing so it is convenient to remark that B (ii) is 
satisfied with « = 1 if max(|@o/éx|, |éo/éy|) = o(e-?) inr, < ae. 

In the Sommerfeld solution for the two-dimensional diffraction by a semi- 
infinite plane we have o = O(y-!) taking the plane as z = 0, y > 0. Since 

éa/éy = O(y-*) = o(y-?), 
B (ii)’ is satisfied. The other conditions are clearly satisfied. Similarly for 
the current densities. 

The Wiegrefe (2) solution for the three-dimensional diffraction by a semi- 
infinite plane does not satisfy B’ (iv) or B, since the current density parallel 
to the edge is O(y-!). On the other hand, the solution given by the author (3 
does, since the current density parallel to the edge is y-te-!*" ©? + O(y), 
The formulae given by the author agree with those obtained for a special 
case by Copson (16). 

In the solution of the diffraction by a circular disk given by Meixner and 
Andrejewski (5) we have, taking polar coordinates p’, 4’ on the disk 
[, = (1 —p'*)tcos¢’, Ty, = (l—p”)-? sind’, o & p'(1—p”)-* cos¢’. Since 
the derivative parallel to the edge is @/¢¢’ there is no difficulty in seeing that 
the conditions are satisfied. Méglich’s (4) solution fails to satisfy B (iv) 0! 
B, since the charge density is Of{(1—p’*)-!} near the edge p’ = 1. 
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BY AN EDGE AND BY A CORNER 


9, Behaviour near a corner 


By a corner we mean the neighbourhood of a vertex O where the direction 
of C changes discontinuously but where S has a continuously turning 
normal. This definition excludes the type of corner that occurs at the 
end of a rectangular wave-guide radiating into space, but we shall extend 
our results to cover this case. 

We have already shown in section 4 that, when C has a vertex at O, the 
field is o (e-!) and that there is neither a line distribution on C nor an isolated 
charge at O. However, the uniqueness theorem of section 6 cannot be 
deduced at once since the analysis of section 5 breaks down in the neighbour- 
hood of a vertex. Here we prove uniqueness in the simple case when the 
neighbourhood of O is plane and bounded by two straight lines and when 
the current and charge densities have certain simple forms. 

We suppose that @ is near one of the straight lines and that, with this 
line as polar axis, the spherical polar coordinates of Q are (p,6,¢), S being 
¢=0. Let the angle between the straight lines which contains S be 6’ 
which, for the moment, will be taken to be less than 7. The fields from 
parts of S which are a distance from 6 = 0 greater than d, where d is fixed 
but small, will be finite. Consider now the potential of a charge density p\ © 
in0< p,<d,0 < 6, < @, where QO is a function of 6, only and © behaves 


like 6} as 0, > 0. A similar behaviour is assumed near 0, = 6’. In order to 
comply with B (iv), i.e. that the charge on a line (not 6, = 0 or 0, = 6’) shall 
be finite, we require that A > —1 and vy > —1. The potential is, apart 


from multiplicative constants, 

dip 0 

p* . | pitl@r-1 dp, dé, 

0 OU 
where r? = p?—2p,(cos @ cos 6, -+sin 6 sin 6, cos ¢)+ 1. Let T' be a large num- 
ber but such that 7’ < d/p. For T < p, < d/p, r- may be expanded in 
inverse powers of p, to give, when A is not an integer, U + Vp'+1 where U is 
a series of non-negative powers of p, cos@, and sin@cos¢ and V is a 
series of non-negative powers of cos 6 and sin 6 cos ¢. For 0 < p, < T when p,; 
isnot near unity and @, not near @ the integrand is finite and the’contribution 
of the integral is U, (say). There remains the range l1—», < py < 1+, 
(n, > 0), 0 < 6, < 6” where », and 6” are small. Expand p4+1 in a power 
series about p, 1; the first term is 

lim © . 
p" i | @r-! dp, dé, = pr} () log( R, Rs) dé, 


1 "1 0 0 


where R,, R, are the values of p,— cos 6 cos 6,—sin 6 sin 6, cos¢+r when 
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Pi 1+7,, l—7n, respectively. Since 7, is fixed and @, 6” small the only 
part of the integrand which is singular is that involving R;. From that part, 
putting O 6, we have, as in reference (3), 
6 
p+ [ 6% log(@2-+03— 200, cos p) d0,+- Wp! 
0 
_ 2arp\*10"*1 cos{(v+-1)(p—7)} Wp 
sin(v+1)z 
when v is not an integer. Here W is a series of non-negative powers of 6 
and cos¢. The remaining terms in the expansion of p+! give terms which 
can be included in W and terms involving 6”*4,..... Hence the potential of 
the charge distribution is 
U+(U,+V-4 W) p14 or+1yar+3 ! 2p ‘tort cos{(v+ 1)(¢ my (9) 
sin(v-+1)z 
Using this expression we now deduce the field due to the charge density 
p}© and the current density whose radial and transverse components are 
(const) p*O,, (const)p**0, where @, = 6%, O, = 6% as 0, > 0. In order that 
the current density normal to 6, = 0 should vanish we require v, > 0. In 
order that the total current on a line should be finite we require A, > —1, 
A, > —1, v, > —1. Since the charge density can be derived from the 
current density by the equation of continuity, A, and A, can be less than 
zero only if Ay = As, vy = vo—1 and a certain relation holds between the 
multiplying constants. We consider the case when A, A,, Ag, v, v, all lie in 
(—1,0) so that the fields may not be finite. It is easy to verify that the use 
of Stokes’s theorem required in section 1 is permissible for these charge and 
current densities. 
To obtain the field derivatives of the expression (9) are taken and here 
we note that 
eU 1aU 12aU 1a 1aW 12, 
Gp’ pb ch’ p 20’ 0ad’ 0 &db’ 6 &d 
are finite for all p, 6 near the origin. After that remark a tedious, but 
straightforward, calculation shows that 
E,, = O(p’), Eg = O(p*0"), E,4 = O(p*6") | ao 
B,, = O(p*+p*), By = O(p*O"+-p*), Bg = O(p*O"+ p*) | 
Consider the integral of the Poynting vector over a sphere of radius p, 
with centre O. Let 6, (> 0) besmalland fixed. Since E, = 0(e-!), B, = 0 (e“) 
from (4) we have for 0’—@, >0>6, 0>0+6, that E, = 0(95'), 
B, = 0(p)*) as py > 0, so that the integral of the Poynting vector over this 
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portion is 0(1) as py 9. For 6, > @ we see by using (10) that this portion 
of the sphere contributes Ofp\+26}*+2(p210j:+- paz)} = 0(1) as po > 0 for fixed 
g,. Similarly for the neighbourhood of the other edge. 

Now take the Poynting vector over the surface of the cone 6 = 8/d’ (< 45), 
oy < p <a’ where d’ (< d) is fixed. It is, using (10), 

Of pA+Ar+2 (S414 Gyr+d) 4 pd+Ae+2(§r+1+ §)4.8(8"4+8"+41)} = o(1) (11) 
as p, and 6 tend to zero. Similarly for the other edge. The edge and corner 
may now be surrounded by the sphere p = po, the cone @ = 8/d’, pp < p<’, 
asimilar cone round the other edge, and the cylinder of radius 6 and centre 
the edge for p > d’. From the above it follows that as py > 0 and 6 > 0 the 
integral of the Poynting vector over this surface tends to zero; hence if 
appropriate conditions are satisfied at infinity the solution is unique. 

There is clearly no difficulty in obtaining the same result when the angle 
of the corner is greater than z but less than 27. The result also extends to 
current and charge densities which are the sums of products of the type p*0. 
We may also remove the restriction that A and v are not integers. For if v 
is an integer s (> 0) the term involving 6+! in (9) is replaced by 


@8+1 log @ o (@¥+1) 


for vy < 0 and hence this is included in the above. If A is an integer L (> 0) 
the only difference is that U contains a term p/*! log p = 0(p**+1) (A < 0) so 
that this, too, is included. 

Also we note that the form of the coefficient of 6”’+1 in (9) enables us to 
show that v, v,, and v, are odd half-integers as in reference (3). 

Finally we extend the results to the corner that occurs at the end of a 
rectangular wave-guide by treating each of the flat portions separately as 
above. There is no alteration for the free edges. For the common edge the 
only difference is that v, may be zero; when vg is not zero there is no differ- 
ence. When v, is zero th : alteration in (10) is that B,, = O(p*+ p* log @) 


and this changes the terms pf+*:+#6"+1, 6”+1 in (11) into pa+r2+2§r+1 log 8, 

6’*1 log 5 respectively. Since the new terms are o(1) our conclusions remain 

valid. 
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SUMMARY 
Using toroidal coordinates, a double space of two regions is constructed in the 
manner of Sommerfeld’s theory. The Riemann potential for a single pole is obtained 
and Jeans’s treatment of uniform finite plane current sheets is applied to the case 


if the circular disk. 


1, Introduction 

THE induction of electric currents in a uniformly conducting infinite plane 
sheet was discussed by Maxwell (1). Jeans (2) extended the theory to uni- 
form finite plane sheets, using the method of multiform potentials due to 
Sommerfeld (3). He discussed two special cases, namely the semi-infinite 
plane and the infinite strip of uniform width. 

In the present paper Jeans’s theory is applied to the case of a uniformly 
conducting circular disk. The theory is briefly described in section 2. In 
section 3 the Riemann potential for a single pole is obtained and the solution 
for the disk, together with a numerical illustration, is given in section 4. 


2. Jeans’s theory 

Following Maxwell, Jeans’s method of solution consisted in first deter- 
mining the currents induced by a sudden change in the external field and 
then dealt with the decay of these currents under the action of resistance 

nd self-induction. 

In applying Sommerfeld’s method of multiform potentials, a Riemann 
space of two regions is first constructed. The first region is to be identical 
with the space containing the magnetic system and the second will contain 
certain images of this system. The field of these images will represent the 
induced field. The current sheet itself must be the branch membrane of the 
space and its boundary the branch line. For every point (x, y,z) in the first 
region there will be a corresponding point (x’,y’,z’) in the second, where 
r=2',y=y', z= 2’, if the sign is taken noe to mean the algebraic 
equality of the magnitude of the quantities which it separates. The plane 
of the sheet will always be taken as the z-plane. 

The required Riemann potential of a single magnetic pole situated at 


[Quart. Journ. Mech. and Applied Math., Vol. V, Pt. 3 (1952)] 
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(x,y,z) must satisfy the following conditions: (i) it must satisfy Laplace’s 
equation; (ii) it must be single-valued, continuous and finite in the doubk: 
space, except in the neighbourhood of (x, y,z) in the first space, where jt 
becomes infinite like R-!, R being the distance between the pole and the 
field point; it does not become infinite at the congruent point (2’, y’,z’) in 
the second space; (iii) it must vanish at infinity. Unlike the Newtonian 
potential, the Riemann potential does not have the same values at corre- 
sponding points in the two regions of the double space. 

Sommerfeld showed that the Newtonian potential due to a single pole at 
(x,y,z) is equal to the sum of the Riemann potentials of two similar poles, 
one at (x,y,z) and the other at (2’, y’,z’). 

In electrostatic and hydrodynamic problems, natural boundaries in the 
form of finite plane sheets can be reduced to infinite plane boundaries, pro- 
vided the Riemann potential is used. It is possible, then, to apply the 
method of images. The images, however, should be in the second region of 
the double space. 

Let Q® be the magnetic potential of a simple unit pole at (ap, y. zy) and 
P(x. Yo, Zp) its Riemann potential, and let Q” be the magnetic potential of 
the currents induced in the sheet by the sudden introduction of the pole. 
Maxwell pointed out that the determination of Q” just after the sudden 
change in the external field is a hydrodynamic problem. Jeans used this 
fact to prove that just after the introduction of the pole we must have: 

QO P(x, Yo: —zZo)— P(x9. Yo: 20): (1) 


al (7) 


old P,(2%. Yo: i Zo) T P, (Xo; Yo: 20): 


oz 

where P,(%»,%.2%) is the Riemann potential due to a unit doublet at 

(Xp. Yo: 29) With axis along the z-axis. 

At any subsequent time ¢, while the currents are decaying, Jeans found 

that 2QW aan sa : 
a I 1(%p; Yo: 2+ Ut) - P, Xo; Yo: —%— et), (¥) 
Cz 

where v = o/27 and a is the surface resistance of the sheet. This can be 

expressed by saying that the image system representing the induced field 

moves away from the sheet with velocity v. 


3. The Riemann space for the disk 

In forming the Riemann space for the circular disk, toroidal} coordinates 
(sb,o,6) will be used. The surfaces % = constant are the spherical caps 
having the circle z = 0, p = a (in cylindrical coordinates) in common. The 
upper and lower surfaces of the disk whose boundary is this circle are given 


+ Bateman (4) gives the geometrical interpretations of these coordinates. 
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by u 7, —7 respectively, while the z-axis is given by o = 0. It will be 
noticed that the value of & suffers a discontinuous change of 27 in crossing 
the ‘principal disk’. 

Consider now a double space of three dimensions with two superposed 
regions as explained above. The passage from one region to the other is 
made when a point passes through the principal disk. The two regions may 
be distinguished from each other by assigning different ranges of values to #, 
namely —7 < % < z in the first region, and 7 < % < 37 in the second. 
Corresponding to a point (ys, 0,¢) in the first region, there is an associated 
point (+27, 0, ) in the second. 

We shall regard our conducting disk as the principal disk of the system 
and it will be assumed that it is situated in the first region of the double 
space. 

The ordinary Newtonian potential at (%,o,¢) due to a simple magnetic 
pole of unit strength situated at (9, o9.¢9) is given by R-1. This potential 
can be expressed (Bateman (4)) in terms of toroidal coordinates in the form 


* sinh u(cosh w—cosh a)-? 


\2r7aR 1 (cosh Oo COs us) (cosh 0o— COS yy)! —— — du. 
cosh u—cos(%—wp) 

x (4) 

where cosh a cosh o cosh o,—sinh o sinh o, cos(d6—dy). (5) 


Using the identity 


‘ { r] 
2 sinh ufeosh u COS(w Wy) 


sinh 4u| {cosh }u—cos $(xs—a)} 1+ {cosh 4u+cos }(~—y)}], (6) 

t is found that (Bateman (4)) 
i P(t. Gg: $9) + P(g +27, oy. bo); (7) 
where P (ub. Fo: bo) R-4+727-!sin-Heos }(—y)sech($a)}]. (8) 


It can be shown that P(». cy. dy) is a solution of Laplace’s equation when 


0 
considered as a function of either (thy. Fo: Do) or (u,o.¢). In fact it is sym- 
metrical in the two sets of coordinates. Moreover, it is a uniform function 
of (&,o,¢) in the double space, since its value is not altered when ¢% is in- 
creased by 47. It is also continuous in the double space except at the point 
ho, Gp, $9), Where it has an infinity of order R-! (it is finite at (y+ 27, a9. do). 
the corresponding point in the second region of the double space at which 
R-1 is infinite). It can ths be seen that P(x. 09, ¢9) has all the properties 
of the Riemann potential of a single pole at (%p, ¢.¢%9). It will also be seen 
from (7) that the Newtonian potential for a pole at a point in the first region 
is equal to the Riemann potential of two similar poles, the first at the same 
point and the second at the corresponding point in the second region, which 
is true in general for the Riemann potential as stated above. 
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4. Solution for the disk 


Jeans’s method can now be applied without change in the symbolical 


notation, except that functions of (a, y,z) will be functions of (4, 0,4). Ifat 
time t = 0a single pole is created at the point (by, 9,9) = (Zo. Po: bo), then 
just after the creation of the pole 
QO = P(—dy+ 27, a9, $9) — P(o +27, a9; $o)- (9) 
[t will be seen from (8), (9) that the current function, which is given by 
Q©/27 at the surface of the disk, vanishes at its boundary o = ©, as was to 
be expected. Similarly, 
EQ 


y 


Cz 


= —Pi(—dy +27, 09, $9) +P (Wo +27, 09.6) att = 0. (10) 


These two images will represent €Q/éz at any following time provided 
they move with velocity v in the direction of the normals away from the disk. 
This means that the expression on the right-hand side of (10) will represent 
EQ /éz at time t provided the point (yo, 79,4) is identified with the point 
(29+ vt, po, do) (Zo; Po: Po): Now 


c C ‘ 
Py (ho, 0, bo) = + 7 P(th9, 9, bo) a {RF (ho, Fp, $o)} | 
Cx Cx ). (11) 
where Fb, 09,9) = $+771sin-eos $(y—a,)sech $a} | 
CF eds CF x éa on 
Hence P, (9; %:%9) = R if 9 ) +F — . (12) 
Ciba CZ x Cay C2 025 
From the relation between toroidal and cylindrical coordinates, 
Z+tpy = acot 4(y +705), (13) 
it is found, by differentiation, that 
1 Cho = cos ys, cosh o,— 1 | 
Cx 
, . (14) 
Cod, ° ° 
a —° — —siny,sinho, | 
OZ» 
Hence 
P, (59, 9, Po) = —R-*(zy—z) Pho, 09.9) + 
+ (27a R)- | {(cosh o, cos )—1)sech $a sin $(ys—ay)+ 
+ $sin f, sinh op sech? $a cos 3(ys—y,) X 
< (sinh oy cosh o— cosh oy sinh o cos(¢—dy))} > 
x {1—cos* (ys—ay)sech? da}-!]. (15) 


It will be seen from the above that the expression for P, (io, a9, bo) is rather 
complicated. If, however, the pole is created on the z-axis at distance % 
from the disk, this expression simplifies considerably. In this case 


o, = 0, by) = 2tan-(a/z,), a= ¢, (16) 














an 


wi 


(ws, 





1bolical 
). [f at 
o)s then 


(9 
ven by 


. Was to 

(10 
‘Ovided 
1e disk, 


present 


7 point 


(14 


(15 
rather 


nce 2% 


(16) 








THE INDUCTION OF ELECTRIC CURRENTS 


and 


P,(hq: 9. bo) R-3(z,—z)[4+727-1sin—{eos $(—yYy)sech $0} |— 


- . j 9 ‘ ‘ ‘ 
, sin $(~—y,)(z§+-a*)-{cosh? $0—cos? $(s—o)} ef 
aR 


— 
~I 
— 


where R, the distance between the field point (%,0,¢) and the point 


thy. To: Po)s iS given by R? (zy—2)®-+ p?. (18) 
If R, is the distance between the field point and tk > point (—#p, a9, do), then 
R? = (z%)+2z)?+ p?. (19) 


Using (10), it is found that 
ea . ' 
[| R-3(z,—z){4—7-1 sin (cos 3(ys—yy)sech $a)}+ 


Ce 
Ry 3(zy+2){4—a7— sin-( cos $(b-+-yy)sech $0)}+ 
L.(a/2)(a®+-22)-4f — R-! sin $(4—a,)( cosh? 40— cos? 3(s—yy))-? + 
Ry ‘sin 3(%-+)( cosh? 40—cos? 4(s+-ahy))-*} | (20) 
it t Q. 
The same expression on the right-hand side of (20) will represent 6Q®/éz 
at any following time provided R, R,, and y, are changed to R, R,, and x, 


respectively, where 


R? = (z,+vt—z)*+ p? 
RF = (z)+vt+z)?+ p? }. (21) 
MT 2tan—(a/z,+7t) 


The potential of the induced field can thus be evaluated in the form 


x 


* aQM 
a | —— dz (22) 
. Cz 
and the current function as 
] Ll £ aQo 
ay Q)' :=— dz. (23) 
) c=0 ») ae 
2a 2a ez 





If the point at which we are evaluating the field lies also on the z-axis, 
then o = 0 and uJ 2tan-l(a/z). In this case using equations (20) and 
21) we get 

an . io 
' 2 3 (ys—Yq)(Z)—2) 7+ (b+ Yq)(%9+2)-*- 
a(a?+-22\-1{(Z,—z)-1—(Z)+2)“}. (24) 

A special numerical case is now considered. A unit pole is created on the z- 
axis at distance 2a from the origin at time t = 0. The normal component of 
the induced field is calculated at a point on the z-axis at distance a from the 
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origin using (24). The results are given in the following table. They sh 
that the normal component of the induced field is reduced approximately 
to 4, zs, and 4 of its initial value after times 27a/o, 9-87a/o, and 307a/@ 
seconds respectively. 
TABLE 
The normal component of the induced field due to the sudden creation 
of a unit magnetic pole at (0,0, 2a) calculated at (0,0, a) 


EQ" 





9 


00009 


2° 


222 “ 0°0554 


° 


15760 : 00500 
1601 5 0°0465 


9 9 


1382 8- 0°0427 


° 


1203 8: 0°0393 
O°1055 ‘ 0°03606 
0°09 34 0°0335 
0°05 32 0°0316 
0°0746 0°0277 








This problem shows that it is possible to obtain numerical results by thé 
method of this paper. Jeans (2) in his discussion of the semi-infinite plang 
and the infinite strip of uniform width did not give any numerical results, 

For any other inducing magnetic system the induced field can be obtained 
by superposition of the fields of its component poles. Also if we have aj 
continuously changing magnetic inducing system, then (following Maxwell 


(1) and Jeans (2)) we can calculate the induced field by integration in some 


cases by letting the change and interval during which it occurs both tend 
to zero. 


In conclusion, I wish to thank Professor A. T. Price for suggesting thif 
work and for his interest in the development. I wish also to thank Dr. A.M 
Gordon for kindly checking the mathematics. 
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